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This  report  is  primarily  concerned  with  the  analysis  of  induced  transient  current 
and  voltage  pulses  on  buried  shielded  transmission  lines,  due  to  eartli  conduction 
effects  of  nearby  lightning  d’seharges.  A numerical  method  is  presented  in  this 
report  to  determine  the  amount  of  coupling  between  a lightning  discharge  to  ground 
and  an  earth- return  transmission  line.  The  transmission  line  is  assumed  to  be  a long 
straight  horizontal  coaxial  cable  with  an  inner  shield  and  an  outer  armor,  terminated 
on  both  ends  with  typical  communication  equipment  load  impedances.  The  general  case 
is  considered  here,  in  which  the  outermost  conductor  is  not  necessarily  in  perfect 
contact  with  the  conducting  earth  but  has  a contact  impedance  with  the  earth,  as 
in  cables  with  an  outer  dielectric  covering  for  corrosion  or  water  protection, 
indirect  strikes  to  the  cable  via  conductive  coupling  mechanisms  through  the 
earth  are  considered.  Several  average  lightning  channel  conditions  and  several 
representative  buried  cable  geometries  are  examined.  The  results  arc  conveniently 
displayed  via  numerous  graphs  of  the  time  histories  and  frequency  spectra  of 
the  resulting  transient  current  and  voltage  surges. 
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CHAPTER  1 


Introduction 

The  work  described  in  this  report  is  part  of  a larger  study  program  to 
provide  protection  for  communication  electronics  equipment  against  transient 
electromagnetic  d i sturbances . Strong  electromagnetic  disturbances  may  be 
caused  by  nearby  lightning  activities  or  man-made  electromagnetic  pulses. 
Current  and  voltage  pulses  are  induced  in  nearby  cables  running  between 
buildings  or  equipment  enclosures.  These  currents  and  voltages  are  then 
coupled  into  the  terminal  equipment.  Excessive  interference  in  conventional 
exposed  metallic  communication  lines  is  indicated  unless  adequate  protection 
measures  are  provided.  This  may  require  extra  electromagnetic  shielding  of 
certain  important  communication  lines  and  associated  buildings  housing  sen- 
sitive equipment  and/or  the  installation  of  protective  devices  on  certain 
communication  equipment. 

The  larger  study,  known  as  the  FAA  Lightning  Protection  Study,  has 
been  performed  by  the  Post-Doctoral  Program  through  several  of  its  member 
universities  for  the  Federal  Aviation  Administration.  The  institutions 
include  the  Air  Force  Institute  of  Technology,  Florida  Institute  of  Technology, 
Georgia  Institute  of  Technology,  and  Purdue  University.  The  individual  par- 
ticipants in  the  FAA  Lightning  Protection  Study  are  listed  in  Appendix  C. 

1 . 1 FAA  L i ghtn i ng  Protection  Study 

The  FAA  Lightning  Protection  Study  consists  of  three  technical  tasks: 

(1)  the  determination  of  the  voltage  and  current  levels  likely  to  be  con- 
ducted to  FAA  equipment;  (2)  the  determination  of  the  susceptibility  levels 
of  FAA  equipment;  and  (3)  the  determination  of  lightning  protective  devices 
that  are  available  to  reduce  the  levels  of  (1)  to  those  permitted  by  (2). 


These  three  tasks  have  been  performed  in  parallel  with  close  interaction  and 

are  essentially  completed.^’  ^ Appendix  C lists  the  schools 

having  primary  responsibility  for  each  of  the  tasks. 

This  report  is  the  result  of  the  work  done  under  the  first  technical  task. 

This  study  is  primarily  concerned  with  the  numerical  analysis  of  the  resultant 

circuit  disturbances  caused  by  earth  conduction  effects  of  lightning  discharges. 

A numerical  model  is  presented  in  this  study  to  determine  the  amount  of  coupling 

between  a 1 i ghtn i ng  discharge  and  an  earth-return  transmiss ion  line.  Only  the 

detailed  numerical  calculations  are  to  be  reported  in  this  report  since  the 

(M 

theoretical  foundation  and  analysis  were  presented  in  an  earlier  report. 

Actually,  this  work  can  be  considered  as  an  extension  of  that  earlier  report 
since  the  earlier  report  only  contained  a description  of  the  theory  of  lightning- 
induced  transients;  whereas,  this  present  report  contains  the  applications  of 
that  theory  to  a typical  buried  shielded  transmission  line. 

1.2  Problem  Statement 

Numerous  interference  and  protection  problems  are  encountered  in  the 
development  and  operation  of  extensive  communication  and  power  systems. 

These  problems  are  caused  by  the  internal  coupling  of  such  systems  with  each 
other  and  by  the  external  presence  of  the  earth  which,  in  some  measure,  is 
involved  as  a return  conductor.  The  earth  also  serves  as  a return  conductor 
for  lightning  currents,  which  often  occasion  disturbances  in  commun icaf ion  and 
power  circuits.  Lightning  disturbances  are  largely  atmospher i c phenomena 
governed  by  the  physical  properties  of  the  air.  However,  the  behavior  and 
effects  of  the  lightning  near  the  surface  of  the  ground  in  communication  and 
power  systems  are  primarily  earth  conduction  problems  caused  by  the  finite 
conductivity  of  the  earth.  Therefore,  problems  arise  both  in  communication 
and  power  system  circuits  concerning  the  protection  of  transmission  lines  and 
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associated  equipment  against  interference  and  possible  breakdown  caused  by 
excessive  voltage  or  current  surges  induced  by  lightning  discharges. 

To  deal  adequately  with  such  problems,  it  is  necessary  to  consider 
solutions  to  the  basic  problem  in  which  the  earth,  as  well  as  conducting 
current  paths,  are  involved  in  the  lightning  discharge.  The  analysis  of 
such  problems  is  inherently  more  complicated  than  the  problem  of  completely 
metallic  circuits  embedded  in  an  insulating  medium,  since  the  great  extent  of 
the  earth  necessitates  the  use  of  electromagnetic  field  theory,  rather  than 
conventional  transmission  line  or  circuit  theory,  in  the  solution  of  most 
aspects  of  the  problem.  It  is,  therefore,  necessary  to  restrict  the  analysis 
to  fairly  simple  fundamental  cases,  in  which  simplified  models  of  the  earth, 
cable,  and  lightning  channel  geometries  are  used,  on  account  of  the  complexities 
that  would  otherwise  arise.  Therefore,  ionization  effects  caused  by  high 
induced  voltages  or  electrolytic  actions  are  not  considered.  Also,  the  heter- 
ogeneous character  of  the  earth  as  a conductor  and  an  electrolyte  are  not  con- 
sidered. Furthermore,  the  extremely  variable  nature  of  the  lightning  currents 
and  voltages  are  not  considered;  however,  typically  average  values  of  the 
lightning  channel  parameters  are  used. 

1 . 3 Overv i ew 

This  report  is  organized  as  follows. 

First,  the  geometry  of  the  cable  and  the  lightning  channel  are  stated 

in  Chapter  2.  Then,  the  basic  equations  which  govern  the  behavior  of  the 

electromagnetic  disturbances  caused  by  earth  conduction  effects  of  lightning 

discharges  are  summarized  in  Chapter  3-  These  equations  were  derived  in  the 

(M 

earlier  theoretical  analysis  and  are  to  be  studied  and  programmed  in  this 
present  numerical  analysis.  These  equations  include:  the  fields  due  to  an 

electric  dipole  in  free  space;  the  fields  of  a vertically  or  horizontally 
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oriented  dipole  above  a flat  earth;  the  mutual  coupling  impedances  between  the 

dipole-like  lightning  channel  and  a buried  wire;  the  equivalent  distributed 

transmission  line  formulas  of  the  coupling  phenomena;  the  induced  electric 

field  intensity  along  the  outer  conductor  of  the  cable  for  a lightning  stroke 

to  ground;  the  Green's  Functions  for  a distributed  voltage  or  current  source; 

character i st ic  equation  of  the  transmission  line;  the  characteristic  values  of 

the  transmission  line,  e.g.,  the  propagation  constant  and  associated  propagation 

modes ; the  induced  current  and  voltage  surges  on  the  outer  armor  of  the  cable; 

the  voltage  and  current  standing  Waves  on  the  transmission  line;  the  induced 

current  and  voltage  surges  on  the  center  conductor  of  the  coaxial  cable;  and 

the  impedance  transfer  functions  of  an  armored  and  shielded  coaxial  cable. 

The  details  of  the  above  theoretical  analysis  were  presented  in  the  earlier 

(1*) 

report.  The  numberical  techniques  used  to  compute  the  theoretical  for- 
mulas given  above  are  discussed  in  Chapter  4.  A listing  of  the  resulting 

computer  code  is  given  in  Appendix  A.  A typical  sample  of  the  input  data  is 
given  in  Appendix  B. 

Finally,  the  output  results  are  displayed  in  Chapter  5 via  numerous 
graphs  of  the  time  histories  and  frequency  spectra  of  the  resulting  transient 
current  and  voltage  surges.  Some  intermediate  results  are  also  presented. 
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CHAPTER  2 

Problem  Geometry  and  Pulse  Model 

2.1  Geometry 

The  geometry  of  the  cable  and  the  lightning  channel  are  shown  in 
F i gure  1 . 

The  cable  is  buried  at  a uniform  depth  d below  the  surface  of  the 

earth.  The  cable  is  coaxial  with  an  inner  core,  an  intermediate  shield  and 

armor,  and  an  outer  sheath  for  corrosion  and  mechanical  protection.  The 

+ 

cable  has  a length  H and  is  terminated  in  typical  load  impedances  at  each 
end,  which  are  grounded  to  the  armor  and  the  shield. 

A ground  stroke  of  lightning  carrying  a total  current  I terminates  at 
a distance  y^  from  the  cable.  The  lightning  channel  is  vertical  and  extends 
to  a height  h above  the  ground.  The  lightning  channel  carries  a typical 
"double  exponential"  return  stroke  current  pulse,  with  a rise  time  on  the 
order  of  1 is  and  a decay  time  on  the  order  of  50  us.  The  peak  current  of 
a typical  return  stroke  is  on  the  order  of  20  KA. 

2.2  Pulse  Model 

The  current  pulses  produced  in  a lightning  stroke  may  be  closely 
represented  by  a "double  exponential"  distribution 

I (t)  - l0(e~at  - e'Pt)  (1) 

The  pulse  rises  rapdity  to  a peak  value  I at  the  time  t^  and  decays  slowly 

to  a value  of  cl  at  the  time  t,,  where  c < 1,  , 

P d 


5 


(3) 


cl  = I (t  ,)  = in(e  atd  - e Ptd)  (3) 

p a 0 

A relation  between  a,  8,  and  t can  be  determined  from  the  expression  for 
the  extremum 

Sllii  - o,  (-> 


which  implies  that 
a/B  = e(“'6)tP 


-CXtn  q -t  tn 

ae  P = 8e  P 


A relation  between  a,  8,  t , and  t,  may  be  obtained  from  (2)  and  (3), 

P d 


(e’atd  . e-etd)  = c(e_atP  - e‘BtP) 


a and  0 may  be  solved  numerically  from  (5)  and  (7)  once  c,  t , and  t,  are 

p d 

specified.  A simple  and  accurate  approximation  may  be  made  by  noting  that 
t^  >>  tp.  Then  (3)  may  be  approximated  by 

cl  * I.  e‘atd  (8) 

P 0 

Combination  of  (2)  and  (3),  yields 


'P  = 'c('  - ^ 


An  expression  for  a may  be  obtained  from  (8)  and  (9), 


In  c + ln(l-^-) 


t - t . 

P d 

When  t,  >:>  t , a <<  8,  the  above  equation  nay  be  further  simplified, 
d P 

a = (In  c)/(t  - t ) (1 

P d 

Knowing  a,  a numerical  value  for  0 may  be  calculated  from  (6)  using,  for 
example,  a Newton-Raphson  technique. 
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Knowning  the  complete  wave  form  I (t)  in  the  time  -domain,  the  frequency 
spectrum  i (to)  of  the  current  in  the  frequency  - doma i n is  determined  by 


i(u)  = i0  ( — 4-  - 3-V) . 

0 \a  - ju)  B - jo)/ 

For  a typical  lightning  stroke, 

t = 1.2  ps 

r 

td  = 50  ps 

and 

I = 20  KA 

Using  the  above  parameters,  the  resulting  wave  for  l(t)  in  the  time- 
is  shown  graphically  in  Figure  2;  the  resulting  wave  spectrum  i (u)  ir. 
frequency  - doma in  is  shown  graphically  in  Figure  3. 

The  computer  program  superposes  pulse  trains  of  the  above  "double 
exponential"  form  to  model  the  prestrikes,  the  return  stroke,  and  any 
sequent  restrikes. 


(12) 


doma  i n 
the 

sub- 
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CHAPTER  3 


Discussion  of  Theoretical  Results  to  be  Programmed 

(4) 

The  theoretical  results  of  the  earlier  report  are  now  summarized 
and  discussed.  Detailed  description  of  symbols  and  definitions  of  terms 
can  be  found  there. 


3.1  Vertical  Dipole  (Free  Space) 

The  current  density  £ of  a vertically  polarized  point  dipole  with 
a total  current  I a distance  z^cnthe  z axis  in  free  space,  as  shown  in 
Figure  A, is  approximated  by 


J_  = z ld£  6(z-Zq)  . 


(13) 


Then,  the  resulting  Hertz  potential  II  is 


n = z j 


. Id£  e 


-jk0R 


4tt  u)eQ  R 


(WO 


where  the  free  space  wave  number  is  given  by 


ko  = “/  Vo 


and  the  distance  R is  given  by 


R = / p2  + (z-zQ)2  . 


Alternatively,  by  applying  a cylindrical  Hankel  transformation. 


n . Id£  I a-  T , /_  \ 

11  = z it dr  -=~  J„(t  pie 

0 0 


^0  z~zo 


[13) 


where  is  the  Bessel  function  of  order  zero  and  where  the  axial  wave 
number  £ is  defined  by 


/t2  - k2 

i K0 
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Notice  that  there  is  no  <|>  dependence  due  to  the  axial  symmetry  of  the  geometry. 


3.2  Vertical  Dipole  (Above  Flat  Earth) 

The  Hertz  potentials  of  a vertically  polarized  point  dipole  with  a total 
current  I a distance  z^  on  the  z axis  above  a flat  earth,  as  shown  in  Figure  5, 
have  the  form 


Jisi-f*  VT  p)  (t  e''a|2"Zo1*  r(,)e"V).(z>0)  07) 


-F  z , 


n =zj 


, .n  r“'  Hz 

n = z j , --  dr  J (t  p)  t (t  )e  e 
- e %toe  J Q 0 


(z<0)  (18) 


where  II  and  H are,  respectively,  the  potentials  in  the  air  and  in  the 
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earth.  Also,  r is  the  ref lect ion  coefficient  of  the  wave  reflected  from  the 
air-earth  interface  and  t is  the  transmission  coefficient  of  the  wave  trans- 
mitted through  the  air-earth  interface.  From  the  boundary  conditions  on  the 
fields  at  the  air-earth  interface. 


T / ka^e 

r(T)  =Tr  (1  - 2 a^-2 

' k F + x £ 
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where  the  axial  wave  numbers  F and  F ate  defined  by 

a e 


/ 2 ",  2 
= / T - k 


r / 2 ,2 

F = / T - k 

e e 

and  where  the  wave  numbers  k and  k are  defined  by 

a e 


ka  ■ ko  ■ “'Vo 


k = m p (e  + j — ) 
e e e w 
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Alternately,  by  applying  a cylindrical  Hankel  transformation, 


-jk  R -jk  R 

n = z j r-1— - 1 + e--,a 

' R R 
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where  the  distances  R and  R are  defined  by 


R = / P2  + (z-zQ)2 


R = / $ + (z+zQ)2 


and  the  terms  A and  A are  defined  by 
a e 


A = 2k 2f  dtx  J (Tp ) 
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k £ + k f, 
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The  first  two  terms  involving  R and  R in  (21)  for  11  represent  the 

— a 

fields  due  to  the  dipole  and  its  image  through  a perfectly  conducting  earth. 

The  terms  A and  A in  (21)  and  (22)  for  IT  and  II  are,  respectively,  the 
a e — a — e 

correction  terms  in  the  air  and  in  the  earth,  which  account  for  the  finite 
conductivity  of  the  earth. 


3 • 3 Ground  Stroke 


The  Hertz  potential  II  i n the  earth  due  to  a vertical  lightning  stroke 

to  ground  with  a total  height  h is  determined  by  integrating  the  Hertz  potential 
in  the  earth  along  the  z axis  from  0 to  h,  i.e.. 
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In  the  above  expression,  it  is  assumed  that  h » 1. 

The  corresponding  electrostatic  potential  ^ in  the  earth  is 


(t)  = (t)  = ~ j 

Hs  e J 2ttc or 


<X' 
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k f. 

dTT  Jo(Tp) 
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and  the  corresponding  field  intensity  in  the  earth  along  the  position  of 
the  buried  cable  is 


E = E 
s e 


where 
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F / dxx2  J (xp)  e e , (z<0)  (27) 

R'<;  J 0 k t + k E 


v=ys 

z=-d 


2ttwe 


0 s ' 0 


k ( + k £ 

e a a e 


/ 2 2 

= /X  + y 


3.^4  Governing  Equations 
3-^.1  Faraday  1 s Law 
When  Faraday's  Law 

fdl_  * E.  = - ju>  J dS  • B 
r 


(28) 


is  applied  to  the  contour  Cj  and  the  surface  Sj,  as  shown  in  Figure  6,  the 
rate  of  change  of  the  voltage  wave  V along  the  cable  is 


dV (x)  = _ . 


dx 


jox^(x)  - E (x) 


(29) 


where  the  voltage  V and  the  potential  (J>,  per  unit  length  in  the  x-direction, 
are  defined  by 
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3.4.2 

Amperes ' 

Law 

When 

Amperes ' Law 

<Sf*L 

; r 

• H = jto  J, 

. dS  . 
’2 

O' 

(32) 

is  applied  to  the  contour  C ^ and  the  surface  , as  shown  in  Figure  6b, 
the  rate  of  change  of  the  current  I along  the  cable  is 


d I (x) 
dx 


j co E E (x)  ( 2tt  r ) 

CO 


where  the  current  I,  per  unit  length  in  the  x direction,  is  defined  by 

I (x)  = j,  dj_  • H (34) 

3.4.3  Transmission  Line  Model 

The  above  equations  developed  from  Faraday's  Law  and  Amperes'  Law  are 
put  into  the  form  of  the  standard  Telegrapher's  Equation  on  a distributed 
t ransmi ss ion  line 


dV  (x) 
dx 


Z I (x)  + E 


= - Y V (x)  (36) 

where  the  series  impedance  Z and  the  parallel  admittance  Y are  def.ned  by 

Z " ZL  + ZS  <37) 

Y = (38) 

1 

The  surface  impedance  Z<.  and  the  longitudinal  and  transverse  impedances 
Z^  and  Zj  are  defined  by 
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(39) 
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(40) 


= . 

7 dTTxT 

dx 


'T1 


+ Z 


12 


(41) 


In  (39),  E1  denotes  the  secondary  induced  field  due  to  the  primary 
return  stroke  field  E . The  longitudinal  impedance  is  later  decomposed 
into  the  sum  of  two  terms  Z j and  Z . Similarly,  the  transverse  impedance 
Z is  later  decomposed  into  the  sum  of  two  terms  and  Zj^ . 

Strictly  speaking,  the  impedances  are  functions  of  x.  However,  only  a 
moderate  variation  with  respect  to  x over  the  major  portions  of  the  wire  is 
expected,  except  near  the  points  where  the  current  varies  rapidly.  Further- 
more, it  can  be  shown  that  when  the  fields  under  consideration  are  travelling 
waves,  these  impedances  are  truly  independent  of  x.  Thus,  a good  approximation 
to  the  impedances  can  be  obtained  by  considering  the  natural  modes  guided  by 
the  coaxial  configuration,  as  shown  in  Figure  7. 


3-5  Natural  Modes 

To  determine  the  exact  form  of  the  natural  modes,  the  potentials  in  the 
earth,  insulator,  and  outer  armor  of  the  wire  are  expanded  as 

n =(  dc  F to)  J (E  r)e'J?x  , (r<r  ) (42) 

— to  J to  0(0  (0 

n . -f  dclF*.  fc)J  AZ.r)  + F".  fc)Y  n(£,r)]e'jC*  (r  <r<r  ) (43) 

— I J I 0 I lUl  * u3  I 
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/.  4*  OD  . . 

d C,  F fe)  H lZ,(5  r)e“JC*  (r>r.)  (44) 

e 0 e ' i 

— OG 

where  F . F'..  F". , F are  the  unknown  expansion  coefficients  determined 
i ’ i e 
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from  the  dispersion  relationship.  The  transverse  wave  numbers  £ , 

to  i 

and  £ are  defined  by 
e 
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where  the  wave  numbers  k , k.,  and  k are  defined  by 
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The  impedances  , Z^j,  Z^.  Z^.j  , and  defined  in  (39)~(4l)  may  be 

expressed  in  terms  of  the  unknown  expansion  coefficients  F , F'  F".,  and  F : 

to  i i e 


I F fc)£  2J  (£  b)e'j'XdC 
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(45) 
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The  propagation  constant  £ determined  from  the  above  expressions  is 
denoted  by  £ which  is  the  eigenvalue  of  the  determinant.  Also,  the  depend- 
ence of  the  expansion  coefficients  Fj , F'.',  and  F^  on  F^  is  attained,  i.e. 
the  eigenfunctions  of  the  determinant  are 

f:  (C ) a + j io  e £ J,  Kb) 

i _ c c c I C fFTl 
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When  the  current  on  the  wire  is  a traveling  wave  of  the  form 


I (*)  = lQ  e 0 = J dc  lQ6 (c -CQ) e 
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where 
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s a constant,  the 
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v£tt'  b ( a +j  we.) 


remaining  eigenfunction 
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F is  determined  to  be 
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Therefore,  due  to  the  presence  of  the  delta  function  in  F , the  integrals  in 
(45) - (49)  for  the  impedances  are  easily  evaluated  to  give 
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3.6  Green's  Functions 

The  voltage  and  current  waves  on  the  outermost  conductor  of  a coaxial 
transmission  line  escited  by  a distributed  voltage  source  due  to  a lightning 
discharge  are  now  summarized. 

The  voltage  V(x)  and  the  current  I (x)  on  the  transmission  line  due  to  a 

distributed  voltage  source  1/ ( x * ) are  determined  by  the  superposition  integrals 

rl 

V(x)  « dx 1 G (x ; x ' ) l/(x 1 ) (64) 

? V 

I (x)  = I dx'G.  (x;x')l/(x' ) (65) 
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where  U(x')  is  the  field  due  to  the  lightning  stroke,  i.e. 


</(x* ) = Er 


and  the  voltage  Green's  function  G (x,x‘)  and  the  current  Green's  function 


G (x,x')  are  known  to  be 
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The  unknown  constants  V , which  are  the  magnitudes  of  the  forward  and 


backward  traveling  waves,  are  determined  from  the  boundary  conditions  imposed 


by  the  continuity  of  the  current  wave  and  the  discontinuity  of  the  voltage 


V = + 


-yx ' r+  +yx ' -2yZ 
e - J e e 
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where  the  reflection  coefficients  F"  are  defined  by 
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and  the  determinant  £ is  defined  by 


t.  = - 2(i-r'r+e'2Yf) 


The  propagation  constant  y and  the  characteristic  impedance  Z are 


defined  by 


I 


Y = vVZ 


(76) 
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where  the  series  impedance  Z and  the  shunt  admittance  V were  defined 
earlier  in  (37)  and  (38). 


3.7  Transfer  Functions 

Once  the  induced  current  and  voltage  surges  on  the  outer  conductor 
( p = p ) of  a coaxial  conductor  are  known  due  to  a nearby  lightning  dis- 
charge, the  induced  current  and  voltage  surges  on  the  inner  conductor 
(p  = P<)  are  determined  via  the  use  of  the  impedance  transfer  functions  for 
a coaxial  cable,  i.e. 


: = z i + z .i . 

x ! P=P>  ee  ext  e i i nt 


= Z.  I + Z. . I . 
fF  p i e ext  m i n t 


where  the  internal  and  external  impedances  are  defined  by 
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Since  the  coaxial  cable  may  be  composed  of  several  concentric  shells 
rather  than  just  one,  as  depicted  above,  the  above  formulas  may  be  used 
recursively  to  obta i n the  current  and  voltages  on  the  innermost  cylinder  due 
to  the  current  and  voltages  on  the  outermost  cylinder. 


CHAPTER  b 


The  Computer  Program 


4. 1 Introduction 

The  "LPS"  Program  has  been  written  in  FORTRAN  IV  language  for  the 
CDC  Cyber  74  computer  to  produce  solutions  to  the  electromagnetic  coupling 
problems  which  are  outlined  above. 

Several  factors  were  given  consideration  in  the  design  of  the  program. 
Efficient  coding  was  used  to  reduce  computer  run  time  as  much  as  possible. 

Full  single  word  accuracy  and,  where  necessary,  double  precision  accuracy 
was  utilized  in  the  calculation  of  special  functions.  Single  precision 
accuracy  on  the  Cyber  maintains  fifteen  digits  of  accuracy.  Maximum  use  was 
made  of  core  storage.  The  sizes  of  the  basic  solution  arrays  (up  to  100  x 100 
complex)  were  determined  so  that  secondary  storage  devices  such  as  tapes  or 
disk  packs  do  not  necessarily  have  to  be  utilized  in  running  the  program. 

4.2  Glossary  of  the  Routines 

The  program  operates  via  a MAIN  program  and  numerous  auxiliary  routines, 
e.g.,  a runstream,  program  sequences,  subroutines,  functions,  and  a library, 
which  are  listed  in  a flow  chart  contained  in  Figure  8.  Standard  I/O  and 
mathematical  routines,  e.g.,  intrinsic  functions,  are  assumed  to  be  available 
through  the  FORTRAN  operating  system. 

The  MAIN  program  controls  the  overall  runstream  and  sets  up  the  parameter, 
type,  and  dimension  statements  and  the  common  blocks. 

The  runstream  calls  the  program  sequences,  which  in  turn  call  various 
subroutines  and  functions.  A library  of  special  functions  and  operations  was 
developed  to  support  the  program. 
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A. 3 The  input  Routine 


Subroutine  PSIO  reads  the  user's  control  information,  prints  out 
headings,  and  obtains  information  for  all  operations.  The  input  cards  and 
their  formats  are  listed  in  Table  1. 


CHAPTER  5 


Numerical  Results 

Several  average  lightning  channel  conditions  and  representative  buried 
cable  geometries  are  examined.  The  results  are  conveniently  displayed  via 
numerous  graphs  of  the  time  histories  (time-domain  results)  and  the  frequency 
spectra  (frequency-domain  results)  of  the  resulting  transient  voltage  and 
current  pulses. 

In  particular,  a typical  RG-58/U  coaxial  cable  buried  one  meter  below 
the  ground  in  wet  soil  is  considered.  The  cable  is  assumed  to  be  shielded 
and  armored  and  is  coated  with  a protective  layer  of  insulation.  The  radius 
of  the  wire  is  0.813  mm,  the  fnner  and  outer  radii  of  the  shields  are,  res- 
pectively, 2.953  mm  and  3-953  mm,  the  inner  and  outer  radii  of  the  armor  are, 
respectively,  If. 953  mm  and  5-953  mm,  and  the  outer  radius  of  the  sheath  is 
6.953  mm.  The  relative  permeabilities  of  all  of  the  insulators  and  conductors 
are  assumed  to  be  unity;  whereas,  the  relative  permittivities  of  the  stabilized 
polyethylene  spacers  are  assumed  to  be  2.26.  The  absolute  conductivities 
of  the  copper  conductors  are  assumed  to  be  5-7  x 10^  mhos/m;  the  absolute 
conductivities  of  the  stabilized  polyethylene  spacers  are  assumed  to  be 
1 .0  x 10  mhos/m.  The  earth  is  assumed  to  be  wet  with  relative  permeability 
1.0,  relative  permittivity  5.0,  and  absolute  conductivity  1.5  x 10  ^ mhos/m. 

A "double  exponential"  lightning  channel  pulse  is  assumed,  with  typical  rise- 
times  and  half-life  decay  times  on  the  order  of  1.2  ps  and  50  v-s,  respectively. 

A peak  current  of  20  KA  is  also  assumed.  The  lightning  channel  is  considered 
to  be  vertical  and,  for  comparison  purposes,  is  assumed  to  terminate  at  d = 0 m, 
10m,  50m,  100m,  200m,  500m,  1km,  2km,  5km,  and  10km  from  the  cable,  at  a point 
near  the  end  of  the  cable.  The  cable  is  assumed  to  be  1 km  in  length. 
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Particular  attention  is  paid  to  the  magnitudes  of  the  voltage  and  current 
waveforms  at  both  ends  of  the  line  and  at  the  middle  of  the  line,  i.e.,  at 
x = 0 m,  500  m,  and  1 km. 

As  an  example.  Figures  9 through  16  contain  plots  of  the  computer  output 
which  resulted  from  the  input  data  contained  In  Appendix  B.  In  particular. 
Figure  9 shows  the  frequency-domain  distribution  of  the  external  excitation, 
i.e.,  axial  electric  field  on  the  outer  conductor  of  the  cable,  as  a function 
of  position  at  various  frequencies  due  to  the  lightning  stroke  directly  over 
the  load.  Figure  10  shows  the  frequency-domain  distribution  of  the  induced 
current  transient  on  the  outer  conductor  of  the  cable  as  a function  of 
position  at  various  frequencies  due  to  the  distributed  external  excitation 
of  the  lightning  stroke.  Figure  11  shows  the  frequency  spectrum  of  the  trans- 
fer ratio  for  the  current  or  voltage  between  the  inner  and  outer  conductors 
of  the  cable.  Figure  12  shows  the  frequency-domain  distribution  of  the  inter- 
nal excitation,  i.e.,  the  axial  electric  field  on  the  inner  conductor  of  the 
cable,  as  a function  of  position  at  various  frequencies  due  to  the  imperfect 
shielding  of  the  armor  and  the  shield.  Figures  13  and  1 show,  respectively, 
the  frequency  spectra  of  the  current  and  voltage  pulses  on  the  inner  conductor 
of  the  cable;  and,  Figures  15  and  16  show,  respectively,  the  time  histories 
of  the  current  and  voltage  pulses  on  the  inner  conductor  of  the  cable. 

Table  2 contains  the  values  of  the  peak  current  and  voltage  and  the 
rise  times  and  the  half-life  decay  times  of  the  pulses  on  the  inner  conductor 
of  the  cable  at  the  loads,  and  at  the  midpoint  of  the  cable  as  a function  of  the 
distance  of  the  lightning  stroke  from  the  cable  with  the  other  parameters  of  the 
model  held  constant.  The  lightning  stroke  cocurs  near  the  end  of  the  cable.  The 
program  is  capable  of  varying  all  of  the  parameters  in  the  model  simultaneously; 
but,  only  this  one  case  is  studied  here.  This  case  is  of  interest  and  is 
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treated  in  particular  since  it  is  very  similar  to  the  conditions  which 
actually  exist  during  a thunderstorm,  i.e.,  the  distance  from  the  cable  to 
the  lightning  stroke  varies  from  stroke  to  stroke;  whereas,  the  remaining 
parameters  of  the  model  are  fixed  at  the  initial  time  of  installation  of 
the  system. 

An  examination  of  Table  2 reveals  that  the  ability  of  the  lightning 
stroke  to  induce  current  and  voltage  transients  on  the  inner  conductor 
of  the  cable  at  the  load  decreases  slowly  with  increasing  distance  between 
the  lightning  stroke  and  the  cable,  as  one  would  naturally  expect.  We  also 
note  that  there  is  relatively  little  difference  between  the  values  corresponding 
to  lightning  strikes  occuring  near  the  ends  and  near  the  middle  of  the  cable, 
which  indicates  that  the  coupling  mechanism  is  relatively  insensitive  to 
axial  displacements  along  the  length  of  the  cable. 

The  maximum  current  and  voltage  peaks  observed  on  the  inner  conductor 
of  the  cable  at  the  load  are  6$2  V and  13.1  A,  and  these  maxima  both  occur 
when  the  stroke  is  directly  over  the  cable. 
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CHAPTER  6 


Conclusions 

This  report  presents  the  numerical  results  of  a lightning  protection 
study  which  was  performed  for  the  Federal  Aviation  Administration  (FAA) . 

In  particular,  this  study  considers  the  lightning  induced  transients 
on  a buried  shielded  and  armored  coaxial  transmission  line. 

A lightning  stroke  to  ground  (not  an  intra  or  inter  cloud  stroke)  near 
a buried  coaxial  cable  can  induce  current  and  voltage  surges  on  the  inner 
core  of  the  cable,  even  if  the  cable  has  an  outer  armor  and  an  inner  shield. 

The  external  coupling  mechanism,  i.e.,  from  the  terminal  ground  point 
of  the  "vertical"  lightning  channel  to  the  outer  conductor  of  the  cable, 
is  either  arcing,  i.e.,  by  dielectric  breakdown  through  the  earth  to  the 
cable  jacket,  for  a "direct"  strike,  or  conductive  enerization  of  the  cable 
jacket  through  the  conductive  earth,  for  an  "indirect"  strike.  The  interna  1 
coupling  mechanism,  i.e.,  from  the  outer  jacket  to  the  inner  core  of  the 
cable,  is  conductive  transfer  through  the  various  metal  skins  and  dielectric 
spacers,  i.e.,  skin  effects,  due  to  the  finite  conductivities  of  these  materials. 

In  particular,  these  lightning  induced  surges  tend  to  cause  excessive 
voltages  to  appear  at  the  ends  of  the  wires;  and,  therefore,  tend  to  cause 
excessive  currents  to  Plow  into  the  terminating  equipment  loads.  Obviously, 
then,  circuit  distrubances  or  failures  are  likely  to  occur  unless  adequate 
protective  measures  are  provided. 

Since  there  is  a tendency  to  move  from  "high  voltage"  electron  tubes  or 
discrete  transistor  circuitry  to  "low  voltage"  solid  state  integrated  electronics, 
the  surge  protection  of  carbon  blocks,  neon  bulbs,  etc.,  (devices  which  are 
presently  in  service)  may  no  longer  be  adequate. 
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FIG.  6a.  CONTOUR  C]  AND  SURFACE  S] 
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LSF/LNF/LTF  — ALF 


SRSC  SRCS  - — NRF 


36 


I 


T 

U. 


>-  u_ 

X <-> 
CO  CC 
1 

V 


f 


(flow  chart  cont. 


< Q 

cc  u_  ct 

co  <£  vr> 


/ / 


/feCi’HL 

/t;  oc  o:  u.  cl 

/lL  ifl  i/i  l ui 


38 


EXCITATION  IV'M) 


r 


1 


50 


I 


CURRENT  (A) 


TRANSFER  RATIO  (x  10 


POSITION  (KM) 


REAL 

IMAGINARY 

FIGURE  12.  INNER  EXCITATION  ( FREQUENCY  DOMAIN ' VS.  POSITION 


CURRENT  (mA) 


W 


I 


FIGURE  13.  INNER  CURRENT  (FREQUENCY  DOMAIN)  VS.  FREQUENCY 
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APPENDIX  A 

PROGRAM  L PS 

4-<INPJT, OUTPUT, PLOTS,  Program  Listing 

+ OLSP,  QLSG,  OCOE  , DC  If  . DCOC*OCO V,  OCIC,  OCI  V,  DFTK  , OFTIJ, 

4-  TA«*Et>=INPUT,TAPEc  = OUTPUT,TAPI7  = PLOTS, 

TAPEtl  = DuSP,  TAPE22=CLSQ, 

4-  TAPc23=OCO^ , TAPE 2h=0CIE , 

<■  TAREc5=GCOC, Ta°E26=DCOV, 

4-  TAP£27  = CCIo,  TAPE28  = OCI  V, 

4-  TAP£29=0FT<,  TAPLiJ^DFTU) 

C # It  Itltt  It  It  It  II  It  It  It  ft  tt  It  It  <t  It  It  It  It  It  H It  It  It  It  t tilt#  tilt  It  ttt  It#  It  t»t#ttt«#tlt  II  «H(t»lt#llt#»tHHt«lt#tt» 
C MAIN  PROGRAM 

c t tit  n ft  tt  it  tt  it  it  it  it  it  tit  tt#  1 1 it  it  t it  it  it  it  t ii  t it  it  it  it#  it  ttnnttt  tit  it  it  ttittttttttttttttttttit  t it  mutt 

C F A A LIGHTNING  PROTECTION  STUDY 

C CDC  CY3ER  7-,  (FORTRAN  u) 

C GEORGIA  T-CH/SCHOOL  OF  E t / A TL A NT  A , G A 3Go32 

C NORDGARO  / 12  MAY  76  (RcVISEO) 

C 

C pAk AMl  TwR  OP  = » GC 

C PARAMETER  UE-~l2*CX=Z2»OF=6-»,DT  = 128 

C PmrAMc.T^5  OA  = 1 3L  , CO=  5 

C PARAMETER  01  = b 

C PARAMETER  OW  = o*+ 

C 


INT^G^R 

L P 

INT  EGER 

OE , OX , OF  , cr 

INTEGER 

DA  ,DJ 

INTlGcR 

C M 

INTEGER 

DW 

INT  lulR 

CP,^P,°C 

INT- GlR 

JLGp,DlSQ 

INTEGER 

SCO. ,0C1F 

INTEGER 

ccoc, ocov 

INTEGER 

CCIC.CCIV 

I NT  _!  GE  R 

OFTK.DFTU 

INTl.Gl.0 

RLS,RCO,FTT,RCI , ~ Tl ,FFt 

I N T G c.  R 

M l S , M CO,  WTT  , r.Cl,WTL,WFT 

int  _ g: p 

T T , T X , T V , TH 

RE  A i_ 

MIN, MAX 

COMPL-X 

P,  C 

OOMPl-X 

< , C 

COIPu  X 

U.  0 

COMPLY  X 

i , IT 

COMPLEX 

PC.7C 

compel  x 

MR 

COM°i_l  X 

pcx, zcx 

COMPEL  X 

Y _■  x , Z j X 

CO  M PE  ^ X 

A N X , A I X 

CO M JL ^ X 

ZNX,  /Fa 

COMPLEX 

ZlM, ZuP 

COM  X 

W J , W H 

COMPl-X 

Jl  M .4 , 0 1 M M 

COMPLEX 

j i : o , j i : w , y c i / , v ; ; ^ 

COMP_c.  X 

Jl  IS , J1IS , YE Ii, Yi  IS 

COMPLEX 

S 

COMPEt.  X 

_ > W , _ 1 I , L * ; 

COMP^*-  X 

A PCM . ,JC 1 , A 3Ci 

C 0 M P L r X 

V 

Jltf.N^ION  =>(  vi  > ♦ r (bl  ) 
JIMlNCION  L(t*»o*Un 

DIM.N  ;i 

ON  <<  ;oll  ,C(  !,o„) 

J1  M_  J II 

n \ U(5,oi),(/(’,fi)  A-i 

jiMi  4jIdn  pc(  *) ,7 c(ei> 

T M N , I 

CM  --(•il' 

c 


0(65.5) 

T X ( 5 ) . T V ( 6 ) ,TH(9) 
V (61  ) 


DIMENSION  A(o3>  , 

DIMENSION  TT ( o)  , 

01  MENTION  X ( ol ) , 

COMMON/A.U/CR,lP,PQ.. 
C0MM0N/AI1/0L3P, OLSO 
C0MM0N/Aa2/JC0E. o c i e 
C0MM0N/A13/JC0C, OCO  V 
COMMON/ A Ih/OCIC.DCIV 
C0MM0N/Ai5/DFT<,0FTU 
COMMON/Bll/IMX.EMX 
COMMON/3l2/ML,RL,SL 
C0MM0N/813/MIN«MAX,M0 
C0MM0N/ai4/NDX,MDX,MYY 
C0MM0N/3l5/N8,NI 
COMMON /820/MN»NM 

COMMON/32  i./ML2»MC0,  MTT  . MCI.  FTL  * MFT 

COMMON/322/RlS,RCO.kTT,FC1,FT(.,RFT 

COMMON/823/MLS,WCO,WTT,WCI,WTL,WFT 

C0MMJN/327/MRC 

C0MM0N/33C/MQX 

COMMON/331/MM 

C0MM0N/C11/TM,WM 

C0MM0N/C12/NK,NC 

COMMON/Clo/rS, TR , TO, TP, TC 

C0MM0N/C14/AX, 8X , CX 

CO MM0N/C15/APS, BFS.CFS 

C0MM0N/C16/ARX,RRX,CRX 

COMMON /C17/ARS»8RS»CRS 

COMMON/CpC/PGPS,PGRX,PGFS 

C0MM0N/C21/PSPS. DCPS. DKPS.MPS, NPS 

C0MM0N/C22/°SKX,0CRX , dkrx.mrx ,nrx 

COMMON/C2  3/PSRS,OCRS,DKRS,MRS,NP.S 

C0MM0N/C24/T5PS, Tops, TOPS , Top s , TC FS 

COMMON/C25/TSRX , TRRX, TORX , TPRX , TCKX 

COMMON/ C2 6/ TSRS,TRRS, TORS, TPRS.TC PS 

COMMON/OiO/SI»SF*NS 

COMMON /Oli/XX ,XF ,NX 

C0MM0N/0l2/YI,YF,NY 

CO M MON/013 /W I, NF.NH 

C0MM0N/014/TI , TF ,NT 


COMMON/01 

5/MS, 

MX 

» MT 

,MW, 

MT 

COMMON/Ol 

6/IEI 

.1 

EF, 

NEI 

C0MM0N/01 

7/1 X I 

,1 

XF, 

NX  X 

COMMON/Ol 

e/iri 

» IYF  , 

NY  I 

COMM JN/Oi 

9/XHI 

♦ I 

WF, 

NMX 

C0MM0N/02 

L/XTI 

.1 

TF, 

NT  I 

C0MM0N/02 

1 /XX, 

XT 

C0MM0N/02 

E/XXI 

,X 

XF, 

NX  X , 

MX  X 

COMMON/OE 

3/XY1 

, X 

YF, 

NYX, 

MYX 

C0MM0N/D2 

4/XPI 

, X 

PY, 

NFX, 

MPX 

COMMON /Ell/RK 

C0MM0N/Ex2/RSM,RSp 

C0MM0N/E13/RAM.RA0 

C0MM0N/E1m/RC 

COMMON/E15/ZO.YO 

COMMON/Fil/UR:.,ERE,U£,E£,OE 

COMMON/Fic/URI,ERI.UI,EI,OI 

C0MM0N/F13/URA, ERA,UA,EA,OA 

C0MM0N/F14/UR2,cP2,U2,E2,02 

C0MM0N/F15/'JRS»£RS»US,ES,0S 

C0MM0N/F16/UP* ,E R1 ,U 1 ,E 1 , 01  A-2 

COMMQN/Fj.7 /URW,  ERM,l)W,EW,0W 

COMMON/Gll/ZLM, ?L° 


ooo  ooo  o o oooo 


COMMON /Hli/PCX* zcx 

COMM0N/H„£/Y3x,  ZSX 

C0MM0N/H13/RNX,  WIX 

CGMM0N/H1*./ZMX,ZPX 

COMMON/ II L /PJ, NH 

COMMON/I1I/JOWW,J1WW 

C0MM0N/I12/JJlH,JlIW,VtIW,YiIW 

C0MM0N/I*3/Ju I S < J 1 1 S * Y C I S * Y i I S 

COMMON/ Ilw/HjnS,Hi£S 

COMM ON /I. 5/tM*/,6Ml,tMF 

COMMON/Ilb/APCW,fiPCI,flPC£ 


EQUIVALENCE.  <<*C> 
EQUIVALENCE  (U,VI 


C 


DATA  CR.lP*dQ/‘.>.6.7/ 

DATA  JLS^.OLSO/cltZZ/ 

DATA  JCO_*DCIE/22.24/ 

DATA  OCOC.OCOV/25.26/ 

DATA  OCIC.DCI V/iZ.S*/ 

DATA  JFT<  ,0F  TIJ/2P.  30  / 

DATA  ) ° / 3 1 / 

DATA  DttCXfJF , OT / L 1 * 3*61*61/ 

DATA  JA,DC/o3,5/ 

DATA  J M / t>  / 

DATA  DU/o I / 

C A *_  l R S 

♦(A,0*JA*00«T r.TX,TV,TH.OM,X,Y,OW, 

♦ P,  J*JP»<*C«'J*V*£  *IT  *FC,7C*MP,D£,DX*OF*DT) 

3 T J J 
£ NO 


$ it  tt  it  ii  it  a i)  ,i  ,i  n 0 <\  >i  :i  u •!  a a ii  tt  /I  n :i  n 00  00  00  it  00000  it  tit  000  it  000000000*000  it  00000000* 


R U N 3 T s t_  A M 

it  it  n it  ii  a it  i n ‘i  n it  it  o it  ii  n ii  it  ii  n n n n ii  n ii  n ini  0 mi  n 0 it  n mi  0 nit  00  tint  it  h ii  $ * it  it  0 nnn  im  ii  unn  » 

subroutine.  R3 


♦ (A,0,3A,P0,TT.TX.TV,TM,oM,<tY,i5^, 

♦ c*D*P3*<*C*J*^*£*IT,FC*ZC*Mc,P£,=X,0FtPT> 


RUN3T<EAM 


I NT  -_Gl  R 

pc 

INTlGcR 

»E *P< , 

JF,  -T 

I NT  .Gc.  R 

P A » PJ 

I M T _ j _ R 

S‘< 

I N T G L R 

INTiG.R 

TT  , r X . 

T \l  , TH 

COMp,_l  X 

p.  J 

COM  X 

< . c 

CJM3lc.X 

IJ  , V 

CuM’L.X 

c_ . I r 

COM-’Lc  X 

= L,ZC 

COM-’LE  x 

•A 

COMJ_-  X 

Y 

DIM.N^ION 
J I Mi.  N j 1 0 N 
JIM.  JiiON 
) I M . N 3 1 0 N 
JIM. *il JN 
JIM  N j 1 ON 
JIMINjION 

DIM.  J j I 0 N 


P ( 3F  ) * C ( 2 b ( 

* (-’X  « PT  ) ,C  <FX  ,PF) 

0 ( n X . 3 T ) , V(PX,OFI 
_ ( 3. ,PFI , IT (F£ , PF) 

PJ ( 3F ) , 7C ( CF > 
ix(r>c| 

A ( JA I , 0 ( =A  , °0) 

T r ( P M » . T x ( 3 M I , T V ( P M ) 


A -3 


Trl  (PM) 


ji m: nsion  * <-■*>,  y r pw ) 

c 

COMMON/'5  ^„  / MLS, ^CC.f'TT, MCI, FTW,MFT 
C 

CILl  P$IJ 

c 

IF  ( ■'itS.Nd  . C I CALL  psLS 
MA.J.PA.FO.TT.TX.Ttf, rw,PM,P,Ot OP) 

IF  ( MCO.  NE.  .u)  CJll  °S  CO 
MA,j,PA,po,rr,rx,Tv,TH,PM,E,Pc.zc,°-:,PF) 

IF  < MFw.  •'U  . J)  CALL  PS  Tl 

»(L«JiPAiDL*tT(TXt’\/tTHtFH,/,y,^K,IT.ir,L,oC«ZC>PEiPi«pFI 
IF ( ITT  . Nr  .01  LAuL  PST 

♦ (A  .J.PA.PO.  T T . T X , TV,  TH  i HM  , , = r ) 

IF(MOI.N_.Lt  CALl  °S  C I 

♦{A, J ,PA  , PC, IT, TX  , TV, TH.PM.t, PC, ZC  , IT, MR, PE  , PFI 
IF(l-r».N£.o>  CALL  P$Tl 

♦ (A,),PA,P0.Tr,TX,Ttf,TH,OM.X.Y,JW.C*\/.i,PCtZC.Pc,PX,PF» 

TP(  IF T .NC ,u)  CALL  P3FT 

MA,0,?A,PO,rr,TX,TV,TH,PM,x,Y,Pw,<,CtU,(/.PX,PF.PT) 

c 

P.t  J JRN 
£ NO 

c*«  turn  * it  n n a a a i n a i>  'i  < > it  ii  a # a a ii  ii  n n ii  n ii  a it  n a mi  u n n n im  n intuit  tut  it  nn  it  mm  unit  it  it  n »n  n*  mu  * 
c PROGRAM  Sr  UUnUC7  S 

emit  mi  ii  ii  ii  ii  ii  it  it  n it  ii  ii  ii  it  n * ii  ii  ii  ii  » « ii  n ii  n ii  turn  it  unit  n tin  it  it » n mmuntintititittnntinu>»»mtmt»m 

SUBROUTINE  P^IC 
C 

C INPJT/JUTFUT 

C 

INTEGER  C- ,l=, °0 

INTO  GOP  PLS,RCJ»tTT  ,dCI»PTl»RFT 
INT_or_P  «i-3»t>ICO»WTT  .WC3.WTL.WFT 
KL  A L M I N • M A X 

COM=L£X  ZlM.ZL0 
C 

COMMOm/A1l,/Cx,LP,pQ 

COMMON /311/IMX.EMX 
COMMJN/B^Z/ML.PL.SL 
COMMON/ Bj.  o / MIN  , MAX,  M® 

C0MM0N/31<*/NQX,MDX,MXY 

C0MM0N/31r /•■n.NI 
COMMON/SCI /MN, NM 

COMMON/321/MLS.MCO.MTT *MCI«  R T L • MFT 

COMMON/ 3 ’Z/RlS.RCC.-'TT. PCI. -Tl.FFT 

COMMON/ dZ  £ /WLS , WCO, W TT , WCI , NTL , WFT 

COMMON/ 3 Z7/MP 

COMMON/03  0 /M  3X 

COMMJN/Cl»/TM.WM 

COMMON/Ci2/NK»NC 

COMMON /C2 l /P&P5 .PGR X , PGFS 

COMMON/CL  1 /P3P3, OCPS. D< PS, MPS.NPS 

C OM, ION/ C2 2/ PS RX, OCR X, D< * X, MF X , NR  X 

COM MUN/C’3/PSKS. OCR  3, OKFS .MRS. NRS 

C0MM0N/C2-»/TSc3,  TR°S,  TOPS,  T p°S,  TCC3 

COMMON/Cc-3/T3PX,TPftX,TCFX,TPRX,TC  = X 

COMMON/CL b/ T 3K3, TRRS, TORS, T PR S, TOPS 

COMMON/Cl t /S I , 3F , NS 

COMMON/ 01 1 /X  I , X F , NX 

COMMON/OiZ/Y I , Y F , NY 

COMMON/J13/WI.NF.NW 

COMMON/C14/T I . T F , NT  A-4 

COMM  ON/ 015/ Mo. MX , MY, MW. MT 
COMMON/Oi-.  / IE1.1ZF.N"  I 


common/  ji  ?/ : *. ; , i xf  , i 
: H.MJN/Dl  “/IYI  »I  YF»NY  I 
CUM  M;n|/Ul  V WI  , I , NWI 
COMMON/  -J2  ./ I T l , I TF  . “IT  I 
COM  M0N/_i.»/RW 
C 0.1M J.i/  *l/-  ^ 

CUMMJN/l . i/RAM,-'  £ - 

common/::*./-  c 

COMMON/  1 E/GJ,  SO 

c^mmo'./-::/j-,_*.  , ) £ , c i . re 
COMMON/Pi.2  / Jk1  ,ER  I ,UI  ,£I  . Cl 
C M 1 JN/r  - n / )r  A , r i-  A , !J  A , £ A , 0 A 
MM  . J/f  14/  l^c»c.K2,Uc»££*Ot: 
. -1"  ■ ,EF$,US.t  . 


OOMMON/FiC/J 

' . 4 L ^ * 4 ' * 4 ' . ± 4 C a. 

COMMON/Fir/ )mw p^,iu,rw,n^ 

coM.MON/r.ii/z 

lM. 2L  c 

1 

r Ok  14  1 ( 

1C  11 

) 

*) 

i. 

FORMAT  ( 

± u ' i 

. L ) 

3 

F 0 R M « T ( 

I 1 L ♦ 

u G 1 G • c 1 

*4 

FORMAT  l 

2G*. 

O.iuliO) 

5 

FORMAT  ( 

v jlw 

u.iOIiO) 

.0 

F 0 7 ( 

1 y 

If. °UT/OllT PUT" ) 

i 1 

FORMAT  ( 

"l". 

MAXIMUM  IT-t-ATIOM 

",  12. 

/ 

4 

4 

MAXIMUM  £FROk 

" 4 1 D'»  Z u • 1 ^ 

) 

12 

FORM  u T ( 

J 4 

LCj  MCOr 

",  IcG 

/ 

4 

It  t t 

4 

LOG  ratio 

“ , 1 3G2t . 1. 

/ 

4 

• • • 1 

4 

LJ'j  SCftLi 

" 4loGZii»i0) 

i 3 

fo-MmT  ( 

• l • • 

l 4 

plot  finimuf 

" , 1 JG 2 G . lo  / 

4- 

4 

plot  faximijm 

" .1=022 . 1j 

/ 

♦ 

• • it 

4 

HUT  RODS 

",  12  0 

) 

i * 

FORMAT  ( 

POLYNOMIAL  DEGRtF 

',  1 20 

/ 

4 

f 

CALCULATE  FtRIVATIYE 

1 20 

/ 

4 

* 

table  order 

",  120 

) 

- 5 

FORMAT  ( 

It  It 

J 4 

lUAORflTURt  GF  I 0 

",  I2u 

/ 

f 

4 

1 NT£  GCA  T ION  GF 10 

",  120 

) 

- 6 

FORMAT  ( 

• « *« 

4 

TITLE  NUMBER 

",  120 

I 

oi 

F ’ R M A T ( 

J 4 

iODc  ST  9 OKi 

",  I2G 

/ 

4 

• 

MO  Jl  COL'PLl  (OUTER) 

",  I2u 

/ 

4 

4 

MOTE  TRANSFER 

",  120 

/ 

4 

4 

Mo  J_  CO'IPLt.  (INNER! 

",  120 

/ 

4 

««  •« 

4 

TOOL  T a A NS M I S3 ION 

",  120 

/ 

4 

II  M 

4 

Mo Jl  INVEPSE 

120 

I 

2 2 

FORMAT  ( 

"0", 

Rl  AO  ST  A OK r 

",  120 

/ 

4 

4 

RlAO  COUPLl  (OUTER) 

",  120 

/ 

4 

• 4 «• 

4 

R^AC  TRANSFER 

",  120 

/ 

4 

4 

^ A G OOGPLF  (INNER) 

",  120 

/ 

4 

4 

-lAO  Transmission 

",  120 

/ 

4 

It  It 

4 

■<  l 3 IN  VLRSl 

",  120 

1 

23 

FOmMA  T ( 

J 4 

•^^iTr  STROKE 

",  120 

/ 

4 

tt  it 

4 

write  couple  <out-_p> 

",  120 

/ 

4 

tt  tt 

4 

•fPITt.  TP  ft  NS  Ft  F 

",  120 

/ 

4 

• t it 

4 

write  couple  (inner) 

",  120 

/ 

4 

4 

WFITl  transmission 

",  120 

/ 

4 

tt  it 

4 

w k i r f INVERSE 

",  120 

) 

27 

FORMAT  ( 

It  tt 

Li  4 

MO0c  ROOT 

",  1 2 o 

1 

3 0 

i >R  M A T ( 

4 

Mu JE  clO  T 

",  120 

) 

jj 

FORMAT ( 

•*r\ 

MULT  I PL  HR 

".1PG20.1J 

, 1 X * "T I ME 

4 

4 

MUlTIcLI-R 

" , 1 PG2  0 . 10 

,1  X,  "FREQUENCY 

3c 

c Q P M A T ( 

• t It 

U 4 

NIJMRlR  of  PULSES 

",  120 

. IX,  "THE 

4 

tt  It 

4 

(UM9FP  OF  DIJLSES 

",  120 

,1 X « "Fkc DUE  NCY 

*4  C 

format  t 

Li  4 

pulse  gu.ss  prestrike 

",  1 PS  2 0 • 10/ 

GUESS  i»E  TURN  ".1PG2C.1C/ 


COM4  1 N" / 
DOMA IN") 
DOMAIN"/ 
00 M A I N” I 


f 

* 

PUL  St  GUESS 

Pf  STRIKE 

1PG20. 1C  ) 

41  format  ( 

U 9 * 

PULSE  P.AK 

PR.ST RI<£ 

1PG2C . 1.  / 

♦ 

«•  « 
♦ 

JlCAY  constant 

STRIKE 

1PG20. 1 J/ 

4- 

9 

DECAY  CONSTANT 

PRESTR  IKE 

1 PG2  0 . 12/ 

¥ 

• • • 
* 

MOO.  O.CAY 

PRE STRIKE 

125  / 

4- 

«•  • 
* 

MODE  N U M 3 2 P 

pRESTRIKE 

I2C  ) 

4 2 FORMAT { 

irv 

p'JLSE  PEAK 

RETURN 

1PS2G  . 10/ 

•f 

» 

3-CAY  CONSTANT 

RE  TURN 

1 °S  2 0 . 1 j / 

4- 

* 

DtCAY  CONSTANT 

Pt  TURN 

1 PS  2 0 . 10/ 

4- 

t 

MO  JE  DECAY 

RETURN 

120  / 

4- 

* 

MODE  NUM 3t  R 

qc  TURN 

120  ) 

.,3  FORMAT ( 

j",1 

p'JLSE  PEAK 

RE  I TR IKE 

lPG2u . 10/ 

4- 

••  • 
» 

DECAY  CONSTANT 

RESTRIKE 

IPS  2 0. 13/ 

4- 

* 

Dt_  C A Y CONSTANT 

^E.oTkIK£ 

1 PS  2 0 . 1.  / 

4- 

» 

MODE  CECAY 

RE  STRIKE 

120  / 

¥ 

* 

MJOc  N J M 0 l ^ 

RE STR IKE 

120  ) 

*♦4  FORMAT* 

Ti ME  START 

£k£5T£I<«_ 

1PS2C.10/ 

¥ 

«•  < 

* 

T I Me.  CISE 

PR.  STR  Kl 

1 PS 2 C . 10/ 

¥ 

* 

TIME  CECAY 

PPtSTRIKc. 

1BG20.10/ 

¥ 

••  • 

• 

TIM.  PULSE 

PRE S TR  IKE 

1 PS2  0 • 1.  / 

¥ 

« 

T I Me.  CONSTANT 

pRlSTRKE 

1 P 32  0 • 10  ) 

FORMAT* 

■ • i 

D v 

TIME  STmRT 

r\F  T Uk N 

1 DSE  0 . 10/ 

¥ 

• • • 

9 

TIM.  RISE 

RE  TURN 

1 nS2  o . 1 j / 

b2  FORMA?  (‘ 


4- 

‘•TIME  CECAY  BIT 

JRN 

", 

4- 

••TIM-  PULSE  Re  T URN 

• « 

4- 

‘•TIM.  CONSTANT  --.TURN 

• « 

tb  FORMAT* 

0 

"TIME  ^TA.tT 

T PIKE 

• • 

4- 

“TIME  RISE  *E3 

trike 

it 

¥ 

"TIME  C^CAY  PlS 

TPIK. 

"TIME  PULSE  ' L S 

TRIKE 

• • 

4- 

"TIME  CONSTANT  RES 

T P I Kt 

M 

50  FORMAT* 

U 

"INITIAL  EXCITATION 

",1=020  . 

A 

0/ 

¥ 

"FINAL  lXCITATICN 

*M°G2J. 

1 

j/ 

¥ 

"NUMUfR  OF  SANDS 

"i  led 

) 

51  FORMAT* 

0 

" i N I T * 3 L ^GSxTlON 

“ . 1 3G  2 J . 

A 

o / 

4- 

"PINAL  POSITION 

•• . I p ; 2 j . 

1 

J / 

4- 

** NUMBER  OF  CAN,  » 

” . i : j 

) 

52  fur  ia  t ( 

0 

"Initial  Fc-EQUEN^y 

" « i " 3 4.  0 • 

A 

0/ 

♦ 

"PINAL  FREOU1NCY 

".1=020  . 

1 

J / 

♦ 

“NUM-jt  R OF  PA  NOS 

”,  120 

) 

53  FORMAT  ( 

wl 

"INITIAL  TIME 

" . 1 =S2u  . 

A 

o/ 

♦ 

"final  TIME 

" , 1 PS  2 4 . 

1 

. 4 / 

4- 

"NUTsr'R  Or  DANGS 

120 

I 

54  FORMAT  < 

0 

"moDl  excitation 

« leu 

/ 

"MUOe  POSITION 

I?u 

/ 

¥ 

"TOD"  POSITION 

",  12  0 

/ 

¥ 

"MOD.  PREOUENCY 

",  120 

/ 

¥ 

“,M  )_  TIM! 

",  12  0 

) 

55  FORMAT ( 

J 

"INITIAL  - X c I t a t ion 

1 2 J 

/ 

♦ 

"FINAL  excitation 

I 20 

/ 

¥ 

"NUM PEP  3 c SAMPl.S 

",  Ii0 

) 

5 o FORMAT  ( 

w 

"INITIAL  POSITION 

120 

/ 

¥ 

"FINAL  POSITION 

I 2 u 

/ 

4 

" NU  MPE  R of  SAMPLES 

120 

) 

5 7 forma r< 

u 

"I  ,ITIAL  F - t C U ■ NC» 

",  120 

/ 

¥ 

"FINAL  Ft.C  . NCY 

'•  r T 

« i _ J 

/ 

¥ 

"NUMOLR  of  SAMt-LeS 

",  120 

1 

5 M FORMAT  * 

"initial  time 

",  1 20 

/ 

♦ 

"FINAL  TIME 

",  T20 

/ 

¥ 

" Nt  Mf>c  P Uc  SAMPLES 

”,  12  0 

) 

bl  FORMAT ( 

t 

"R». OIOS  Hl-e 

",*»r,20 . 

4 

J ) 

RADIUS 

V.  ) l'l~ 


• : • i > 
^ M i _ L "* 


'I  I MU  3 
•>4.1)3 

irv"3 


. ' s _ u 


. 13/ 

. l<i) 

. : a / 


1DS2  u • 10  / 
1PG20. 13/ 
1PS2C .13) 
1PG2C. 13/ 
i.P320  . 13  / 
1PG2C . 1.  / 
1 PG2  u . i:  / 
1 cr>2  0 . 10  > 


A -6 


o3  F 0 1 ii 7 ( ** 


"y«ai!»s  apmck  plug  ", 

1 U • 

6-,  eormat  c*  ", 

"RADIUS  SHEATH 

1 JG20  . 

o 5 format  <*v*» 

"Ol  PTH 

• • 

* 

1PG2Q . 

♦ ••  •*, 

"SEPARATION 

1RG20  . 

/i  format  c*u", 

*RELAT IVE 

PERMEABILITY 

l A X TH 

^ M «• 

•re  lat IVE 

PEPMiTTi vi r r 

lARTH 

^ »•  »• 

•absolute 

°EFME ABILITY 

E AR  TH 

♦ " ", 

’AbSOLL  Ji 

BtRMITTIVlT Y 

l A R T H 

^ il  «l 

•ABSOLUTl 

CONDUCTIVITY 

lARTM 

72  FORMAT  ("u  ", 

•ptuATH/c 

PlBMlAEIlITY 

SH.ATrt 

•PLLATIVE 

PERMITTIVITY 

SHE  ATM 

♦ ••  " , 

"AcjJLl/Tt 

PERMEABILITY 

SHl ATH 

♦ " ", 

•absolute 

PERMITTIVITY 

SHEATH 

v " 

•ABSOLUTE 

CONDUCTI VI TY 

SHAATH 

73  FOkMAT(’*v", 

* <c  L A T I VE 

D£  °M£  A SIL I T Y 

A RM  JR 

• y^L A T IVl 

PEkMITTIVITY 

A RM  JR 

•AHSClUTl 

PEP  ME  ABILITY 

ARM  JR 

^ ••  M 

'absolute 

PERMITTIVITY 

A RM  JR 

^ M II 

•ABSOLUTl 

CONDUCTIVITY 

ARMOR 

7>,  FOR1AT  <"o  ", 

V *l AT I Vl 

PERMt ABILITY 

2 

♦ " " f 

• y_*.A  T I VE 

permittivity 

2 

♦ ••  " , 

•m  JSOLUTE 

PERMEABILITY 

2 

^ «•  M 

•a  SSOLUTl 

PERMITTIVITY 

*> 

L 

‘AdSOLUTt 

CJNCt CTIVITY 

2 

75  FORMAT  CO", 

U.lLI  T 1 V/f 

°FK ME  ABILITY 

SHIELD 

♦ “ ", 

•RELAT  IVE 

PERMITTI VITY 

ShIElO 

♦ " ", 

•ABSOLUTE 

PcP  ME  ABILITY 

SHIELD 

♦ '•  ", 

•A  USJLUTl 

PERMITTIVITY 

SHIELD 

•a  ISOLUTl 

CONDUCTIVITY 

SHIELD 

7b  FOwMATCO". 

VlLATIVE 

p£RMl ABILITY 

▲ 

♦ 

*R _ -AT IVE 

permittivity 

1 

* A P SOL  L T E 

PE RM -ABILITY 

i 

♦ , 

•Absolut-' 

Pfr-MITTIVITY 

1 

’absolute 

CONDUCT  I VI  T Y 

X. 

77  FORMAT 

*,  L A T I V* 

PR* Ml  ABILITY 

WIRE 

t " ", 

’ : u « I 1 V l 

PL  *•  T I T T I V I T Y 

W I RE 

♦ " 

* * • _>  0 L U T • 

P1  - Mt  A P I L 1 T Y 

WIRE 

♦ " ", 

•A  , SOLUTE 

P:.  RMITTIVITY 

WIRE 

♦ **  " f 

’ll  . L ) ' T l 

CONDUCTIVITY 

nf  I ^ r 

a:  formatcv, 

"lJAO  IM^EDANCc  MINUS'*. 

l ( i JG2 

♦ ••  •*, 

"LOAl)  IIP 

edancl  plus  ", 

2 < 13G2 

r t a d c c y , l i 

1 MX, ■ “X 

r-aji:-1.  ■> 

*1L,4L,SL 

R ; A 0 ( t •>  , -) 

M IN,  M A X , 

MO 

•k  : U j { ^ -J  % « ) 

N-X ,M*X  , 

MXX 

P.ATicy,  i) 

•n , n x 

( y » 4 ) 

IN  , N M 

yc  a oi  l ) 

UlS.MCO, 

-TT.MCI.MTL.MPT 

Rladic-.  i> 

y L > * O • 

RTI.RCI.RTL.RFT 

y f a ) it p . i • 

W i_  > * w C O t 

wTf**)(CIfWTL«'yF"r 

yt a jicp,  : i 

*1- 

r„a ni.p,  ii 

If  < 

yt-AO  ('  ■>,  w i 

r m , wm 

t-  A J ( Z'i  • 1) 

. * ,‘,r 

yf.  ao  toy.  z ) 

^ j t Pj* 

x « 3G*<  S 

a ) i ;y . 5 » 

’-PS. oc° 

S.LKaS.MP;,N  PS 

v A)I.O,  5 1 

X ,OT- X , M*  X, NR  X 

1 A O ( C y , 5) 

° >ys, dor 

i.OKRS.MR  o.Nys 

c) 

r ,J‘. ryps. tips, teps, tops 

•y  r A J ( .y  » 

x , t„rx , tpr  x , rc 

*<  X 

A tl  ( o ^ * I 

r jy  s . tp.- 

J.TDRS, TBRS.TCyS 

y t a 3 ( c y , •>) 

) A • ->  F * N J 

yE A ) ( w < , -) 

A i , X E , N X 

lu) 

10  > 

10/ 

101 

",  IPG 20. 10/  ! 

", IPG  20. 10/  I 

“,1PG20.1C/  i 

", 1PG20 . 10 / ' 

•MPG20.10) 

", 1PG20. 10/ 

", 1PG20.10/ 

", IPG 20. 10/ 

", 1PG20 . 10  / 

",  1 PG 2 L . 1 0 > j 

", 1°G20.1*/ 

",  IPG2 0 . 10  / 

", 1PG2U.10/ 

”, 1PG20.13/ 

", IPG 20. 10) 

", 1PG20.1C/ 

", 1PG2J . 10/ 

",  1 PG  £ u . 1 0 / 

",  1.PG20.10/ 

1P02C . 10 l 
", 1PG2C.  10/ 

", 1PG2C . Il  / 

",  IPG 20. 10/ 

", lpG2u .10/ 

", 1PG20 • 10  ) 

", IPGOO.Ij/ 

", 1PG20. 10/ 

", 1PG20 . 10  / 

", 1PG20 . 10/ 

", lPG2u  . 10  ) 

", 1PG2G . 10/ 

", 1PG2  0.1  it 
",  lPG2ii • 1 * / 

",  , p',20 .10/ 

1 DG  2 1 . 1 * ) 

. 10)  / 

• 4b  ) ) 


A-7 


...  ..  .. 


r 


■X  . 1 ( . •> , 

n.n  ,-iv 

a ; ( l ■>. . 

w i » .*  r , n w 

• ( •-  < , 

T c , I p . JT 

1 1 

!-> . , my  , Mw  , mt 

■ ^ a j t x . 

i J 

1.  I.l-P.V  I 

> ' - 1 L < • 

1 ) 

IXi.ix17**.*! 

X i ) (l  \ < 

lYI.lYF.MYI 

x -Me-. 

; ) 

1 A I . I 4 F , N W T 

A .)  < , 

_ ) 

c Ti , i rr, nti 

AH  x , 

» j 

AW 

HI  A, 

£ ) 

A _>  M , *■  S B 

7 A j (C\» 

) 

•X  A M , » A = 

.AH  H, 

L t 

?S 

. « J ( H, 

10 

X . A ')( 

L ) 

5 0 

ft  ( !■?* 

1 H £ . r.  A £ . 1 ' 

■<i  A3  ( 3 A , 

. J 

u x : , x i , o : 

..  )(H, 

J*"  A » c .*t  i 3 A 

r . A 3 ( . a , 

. 1 

U a 

X w A .H  ~\ . 

- 5.  OS 

A ( 0-X . 

l-l.Cxi.Oi 

: A lI'Hl 

) 

. J w . t 3 W • 0 W 

'■  - ( . 7 , 

(_  ) 

/c  1 

A ) 1 . - , 

- 1 

7 L 3 

S-?C 

A ) (0^» 

_ ) 

HO  ?IF<MI0.l  .,0)  3 £ T U A.  N 

A-i r ( L " 

l lu) 

I r _ ( lF 

, 11  I 

I 1 X • T M X 

A*  1 ,Il3 

♦ let 

ml . pl , sl 

w - i r i l 3 

.Hi 

MIN, MAX, M 3 

(L~ 

. , o 

MJX 

W - I T £ < l ° 

• A 5 • 

N0.NI 

W F I I c ( L 3 

, lot 

MM 

W - I T 3 ( L 3 

i is) 

111 

w- 1 r_ cl 3 

.cl) 

f'LS.f'CO.MTT.FCI.t'TLtf'FT 

A-  I I.  ( L-= 

,cct 

' cT  i - Cl  t Ti-11c.Tl(cTT 

A I f _ < L ~ 

, £3) 

mSi’A  CO  * W T T inCIiNTl.AFT 

ii  . '■  i ( L 3 

. £7) 

IP 

W~»  I T £ CLP 

. J C t 

1 3 < 

wr i r £ (lp 

.-!> 

T M , W x 

- I Ts.  (4.  o 

. c ) 

i<  , \C 

A'  I Ti  < L ° 

.-cl 

3oPS , P0F  X , FQL 5 

W-  1 1 . t LP 

. - 1 1 

FiFS.OCPbtOKFs.MCi.NFS 

a rtitP 

. - £ 1 

3>PX.0C?>,0k;FX,MxX.n3x 

wpi  i i l 3 

, ^-t 

• )jr)<  5 ♦ WPS • NPS 

w ■ I r t ( l p 

T jF5. TaFl. TOF3, IPPS»TC”S 

W-IT-  (_•> 

.Lb) 

rJ--,r^'x,Tui-x.T33x,Tc^x 

£ ( L P 

. LCt 

TSpS*Tr\,PS*TOKS«  T°-?S*TC^S 

w -r>  i r £ < L 3 

oi.1 

bl.SF.NS 

W ■ 1 1 £ ( L 3 

.bit 

x r , x c , nj  x 

A x I T c.  ( L a 

• 1) 

Y I , YF , MY 

At  1 1 £ ( L* 

, 5>  2 > 

A I . W c , N W 

I-  i I - ( L ° 

.tit 

r i , r f , n r 

a~ i r _ c l 3 

, y*.t 

MS  * MX, MY ,MW, MT 

A-1 Ti (lP 

. 15  > 

1 £ I , I £ F , M£  I 

( . ■ 

• 3b) 

IXI . I X F , MX  I 

(Rll  f L ’ 

, 561 

1 Y1  . I Y c , N Y I 

• ■ I Tt£  f L.P 

« b 7 1 

IWI , I A F , MW  I 

W-  ITu  ( L ° 

. 50 1 

it:,itf,mt  i 

■ IT  ( t 

t 6 1 ) 

rr  A 

w - 1 r £ ( L p 

« ©21 

-3M.A3P 

I 


c 


c 

c 

c 

c 


c 


c 


c 

c 


c 

c 


c 


c 

c 

c 

c 


c 


WRI TE  <Lp, 6 51 
WR  I TE  ( L ° , c h ) 
WRITE ( L P » 6 5 ) 
WR I T £ ( L ° , 7 1 ) 
W R I T£  (L3i  7cl 
W R I r £ ( L 3 , 7 i ) 
WRITE(L°» 7* I 
WRI  T c (L P,  75) 
WPITc(LP,  7 c ) 
WRIT£ (L °,77) 
WRI  Tl  lia. 61) 


* A M , R A D 

KS 

□ 0,  SO 

URE , ERE  » UE  , Et  * OE 
URI  ,tRI ,UI ,EI,OI 
URA,£RA,UA,£A,OA 
UR2 ,ER2,U2.E2.02 
URS,ERS,US,ES,0S 
0R1 .ER1 ,U1 , E 1 , 01 
URR,ERW,UW,EW.0W 
ZLM  , ?LP 


RETURN 

£N0 


3U9R0UT iNt  p$L3<£.0,PA,P0,Tr,TX,TV,TH,PM,p,Q,PS) 
LIGHTSING  STRGKc 


I N T c i»c  R 
INT  EGER 
INTEGcR 
INT  EGt  R 
I N T c G_  R 
CO  M°  L c X 


°s 

PA  ,30 

DM 

C"\ » L P,  PQ 

TT,Tx, TV,  TH 

p.Q 


DIMENSION  P< P£) , 0 (PS > 

0IM_NSI0N  A< PA) , 0 (°A , PO ) 

DIMENSION  T T ( P'*)  f T X ( c M ) , T V(PM»  , TH ( pm  t 

COMMON/Alu /C  R » L P » c Q 
C0MM0N/C11/TM,WM 
C0MM0N/C1 2 / JK,NC 

COMMON/OIS/Mc,MX,RY,MW,MT 

10  FORMAT ("1", "LIGHTNING  STOKE") 

CATA  Ni/i/ 

OATA  TO.WO/u.u.i.O/ 

WRIT£(L°,I0(  ICALL  SRSC  t C A L L SPLS<p,T,pS> 

CALc  TPr,X(A,0,Pf,P0,TT,TX,TV,TH,Pr.P,PS,T0,TM,NK,MT,Nl,NK,Nl> 
CALL  TPGX(A,O,PA,PO.TT,TX,TV.TH,PR,0,3S,WJ,WM,NC,MW,Nl,NC.Nl) 

RETURN 
c NO 


SUBROUTINE  P SCO ( A ,0, PA .Pu,TT,TX,TV,TH,°M,E,PC,7C,P£.°F> 
COUPLING 


INTEGcR  =c,PF 
INTEGER  P A * P 0 
I N T c Gc  R D M 
INTEGcR  CPfL3,0! 

INTEGcR  OCOE * OCIE 
INT-.G.R  R'LS,rCO,PTT,RCI  .FTL.RFT 
INTEGER  WLS,WCJ,WTT,WCI,WTL,WFT 
INT  c Gc  R rT,TX,TV,TM 
COMPLcX  c. 

COMPLcX  PC, 2C 


A-9 


OlMcNilON  c(Pc,pF) 


I 


ooo  ooooo  o o ooo 


r 


■ 


DIMENSION  PC ( °F ) I?C(PFI 
DlN.NsION  A<  PA  > . 0 (Pi  , F0> 

DIMENSION  IT  < PR)  ,TX<CM.TV(PM»»TH(PM> 

c 

COMPON/flit/CP.LF.c-i 
COMMON/ Ax  2/0 COc  . DC  IE 
COMMON/B^E/RlS,  RCO  * HTT  , PCI  • R T L » RF T 
C 0 M M 0 N / 3 1 / ^JLbi  WCOt  rt  TT  ,WCl,WTi_,WFT 
COMMON/B’X/MM 
COMMON/  01  o /l  I ,£F  ,N-1 
C0MM0N/D1a/RI,RF,NW 
COMMON/ DX5/ME  , MX • MY  . MW, M T 
COMMON/OXo/Iti,IfcF,N£I 
COMMON /Oi P/INI, IWF, NW  I 
C 

11  FORMAT  ( “ PRO  P AG  A T I ON  P £P  A M£  T C R S " ) 

12  FORMAT  C’X'VCOUPL  ING  (OUTc.R)*’) 

C 

MM  = 1 
C 

HRITE(LP.ll)  *CAuL  SPSft .PC,/C,P£.PFI 
RRIT£(lP»1c> 

CALL  tpgk 

♦ (A,OtPA,PO,rT,TX,Ttf,TH,PM,-tPt.,PF, 

+ E I » EF,  ME,  ME  , I£1  , IEF,  Nt  I , W X , WF  , NW  , MW  , I W X , I WF  , NWI  ) 
C 

IF(«ICJ.Nt_.U>  CALL  FTrW(DCO£,L,o£.PF) 

C 

RETURN 

iNO 

SUBROUTINE  p^,T(A,C,PA,P0,TT,TX,T\?,TH,pm,MR,PF) 
TRANSFER  (OUTER/INNER) 


I N T cGlR 
INTEGER 
INTEGER 
INTEGlR 
INTEGER 
COMPLEX 


° F 

PA,  PO 
PM 

CR.LP.PQ 
TT , Tx , T V, TH 
MR 


DIMENSION  M / ( PF ) 

D I Me.  N j I ON  A ( PA ) ,0 ( PA , FO) 

DIMENSION  TT(PM>,TX(PM>,TV(°M),TH(PM) 


COMMON/Axu/Ck,LF,FQ 
COMMON/D13/RI,RF ,NW 
C0MM0N/Di5/M£,MX,MY,MW,KT 

10  FORMAT  (’WTRAKSFcR  ( OUT  t R/ I NNE  R ) " ) 


OATA  N 1 / 1 / 

HR  I T£(LP.10t  {CALL  SRTfMP.oF) 

CALL  TPGX ( A,0,  PA , OO, TT . TX , TV, TM.PR, MR,PP, WI, MF.NW, MN,N1 , NW.Nl » 


RETURN 
£ NO 

SUBROUTINE  PSCI ( A,0, FA  tP0,TT , TX , TV, TH,PM,£ ,PC,2C» IT . MR,PE,PF» 
COUPLING 


A-10 


o o o 


c 


INTEGER  Pc.,PF 

INTEGER  a6i D0 

INTEGER  ° M 

INTEGER  ck.lp.pq 

INTEGER  0C0l,0CIE 

INTEGER  RLSt RCO,RTT,fiCItRTL,RFT 

INTEGER  WLS, WCO, W TT, WCI , WTL , WFT 

INTEG.R  TT,TX,TV,TH 

COMP<_eX  E,IT 

COMPL-X  PC.ZC 

COMPLEX  M~ 

C 

DIMENSION  e(Pe,°F)  ,IT(F£  ,FF) 

JIMZNSION  PC(PF) ,7C(°FI 
QIMENSION  MR(PF) 

01  MI  NS  ION  A ( PA ) i 0 (PA i D0 ) 

JIMlNjION  TT(PMI,TX(PM),TV(PM»tTH(PM) 

C 

COMM0N/AlJ/CR,LPt°Q 
COMMON/A^/DCOc  , DC  IE 

COM.MON/32E/RLS»RCO,PTTtPCI.RTL»RFT 
C0MM0N/12?/WLS« WCO.WTT,WCI»HTL. WFT 
COMMON/931/MM 
COMMON/Dli /XI , XF  ,NX 
COMMON/Di i/A I ,WF .NW 
COM  MON/ DIE /Me ,MX,MY,MW,MT 
COM MON/ 01 7/1 XI , I XF, NX  I 
C0MM0N/01 9/1 WI, IWF.NWI 
C 

11  FORMAT  ("I'VppOPAGATION  PAR AMe T£ RS" ) 

*2  FORMA  T ("i'V'CJJPL  ING  (INNER)") 

C 

mm=o 

c 

W-vI  T - (L°,  11)  TO  A L E SRCClf  ,PC,2C,IT,MR,Pt,PX,PFI 
W*  ITE  (LP,  U) 

GAEL  T P G N 

♦(A,J,PA,PO,TT,TX,TV,TH,PM,E,PE,PF, 

♦ Xl,XF,NX,Mx,:Xl,IXF,NXl,WI,WF,NW,MW,IWI,IWF,NWI> 

C 

IFJWCi.N-.  ,U>  CALL  FTCW(DClL.t  ,PF  ,PF) 

C 

RE  T JRN 
eND 

C*  

SU3R0UT INL  °STL 

♦ (A,J,PA,POtTrtrx,TV,TH,PM,X,Y,°W,C.V,i,PC,ZC,P!r,PX,PF) 

TRANSMISSION  LINE 

INTEGER  P;.Jx,PF 
INTEG.R  PA  * °0 
InT.GlR  pm 
INTEGER  cW 
INT  l6. R C" »e  p. PQ 
INTlGlR  OCOr.  • JO  IF 
INTeG.P  OCOC.OCOV 
X N T l Gc R DC  I C « OC I V 

I N T e G e R RLS»RCO*FTr ,ACI*PTL,PFT 
IMT.GeR  Wl  Jt  WCJ,  WTT  , WCI  , r.  TL  , WFT 
iNT.GeR  TT.TX.TV.TH  A-ll 

COMPL- X C.V 
COMPL.  X * 


I 


o o o 


c 


c 


COM^L^X  PC.ZC 
com^l_x  y 


J1  Ml.  Mi  ION 
DIMENS  ION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


C (PX . > , V ( FX , °F| 

Lt°:,an 

PC(  pf>  ,7C(PF) 

a ( pa » , o (Pa , »o) 

Tr  (-‘>'1)  , TX  (“Ml  , TV<  aMI  ,TH  (OM> 
A ( J«  ) , X ( 3W  > 


COMMO«l/fixC/Cr<,<.P.rO 
COMMON/AlE/OCO;.  , DC  IE 
COMMON/AiZ/DCOC, DCOV 
COM, ION/ A 1 */ DCI C, DC  IV 

COMMJN/3E£/RLS,fiCO,RTT,FCI.RTL,FFT 

COMMON/TSc/MLS.WCO.WTT.WCI.NTu.WFT 

COMMON/921/MM 

COMMON/DIi/Xl.XF.NX 

COMMJN/Oit/Y I , YF ,NY 

CUMMO  J/Oli/WI.NF.NW 

COMMON/l/x5/M_  i MX  , - Y«  MW,  MI 

COMMON/DI  7/ IX  I , IXF.NX  I 

C0MM0N/3*M/IYI,IVF,MYI 

COMMON/JH/INI.IWF.NWI 

c 

i.0  format "transmission  line'** 

c 

IF  ( (RTL.Nl.C)  .AND.  (Mm.fC.1)  ) CALL  FT  OP  ( OCOc. , E , PE  » pF ) 
IF( (PTL.Nr ,U > .AND. (MM.Eu.C M CALL  F T DP  { G C I E , £ ♦ °F  , PF  > 
C 

WPIT^(L3.Ij>  JCALL  ScTL(X«Y«PW*C.V.t-.PC.7CfD£.0X.pF( 
C 

IF  (MM.  El),  i)  CALL  TPGK 
♦(A.O.PA.PO.rr.IX.TV.TH.PM.C.PX.PF, 

*■  XI.XF.NX.MX.IXI.IXF.NXl.FI.WF.NW.MW.ir/i.lWF.NWI) 

if( MM.sa.t)  call  r pgk 
♦(A.o.^A.Do.rr.rx.fv.TH.PM.c.ox.PF, 

♦ YI.YF.NY.MV.IYI.irF.NVI.WI.WF.NW.MW.IWI.IWF.NWI) 

IF (MM.tO.C)  CALL  T PG  < 

♦(A,0,DA,PJ,Tr,TX, TV, TH.PM.V.FX.PF, 

♦ YI.YF.NY.MY.IYI.IYF.NYI.FI.WF.NW.MP.IWI.INF.NWI) 

C 

IF ( (WTL.Nt.O) .ANO. (MM.EO.IM  CALL  FT  OW  ( OCOC. C.°X ,pF) 
IF  ( ( wr*..NE  .0  ) . AND.  (MM.tU.  C > ) CALL  FTDW(DCIC,  C,°X,PF» 
IF ( (ATL.NE.L) .AND. (MM.EG.C) ) CALL  FT DW ( OCI V , V , PX , °F > 
C 

RETURN 

END 

SU3R0U T I Nl  PjFT 

♦(AfO,PA,PO,TT,TX,TV»TH,PM,X,Y,PF,X,C,U*V»PX,0F,PT) 


FOURIlR 

TP AN5FOFM 

INT£Gt < 

p» ,PF.°T 

INTcGuR 

= A,JO 

INTEGER 

= M 

I fj  r L G c.  p 

PW 

I N T l5lX 

CP ,L  P, PC 

INTEGER 

DCIC.DCI V 

INTEGER 

OFTK.  OFTIJ 

I N TE Gt  R 

PLS.RCO.PTT 

INTEG.R 

MLS.WCO. WTT 

INTeGl  < 

rr.rx. rv, tm 
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o o o 


complex  k,c 

COMPl-X  U,Y 
COM°Lt  X Y 

OIM£NSION  K < Px , FT ) ,C ( PX . PF ) 

OIM-NSION  U(PX,PT),\/(PX,PF* 

DIMENSION  A < dA ) ,0 (PA , OQ) 

DIMENSION  TT(pM),TX<PM),TV(PM),TH(pm» 
DI  Me.  No  I ON  X(3HI  , Y MM) 

COmmon/aiu/cr.lf.fq 

coMMON/4iw,/ncie,nciv 

COMMON/Ai 5/OFT <*DFTU 

COMMON/  -3  2e./RLo,fiC0,PTT,RCI,RTL,RFT 
COMMON / f3£3/MLS*WCO»WTT«WCI»MTL»WFT 
COMMON/Oi£/YI,YF,NY 
COMMON/QIh/T I , T F , NT 
C0MM0N/D15/M£ , MX ,MY,  MW,  MT 
COMMON/Oie/IYI,IYFf NY  I 
COMMON /Q£C/i T I , I TFfNT I 


10  F ORM  A T ( ” 1 1 


i F T " I 


IF  { RF  T • N£ . U ) CALL  FTCMCCI  C,C,°X,CF> 

IF(RFT.N£.0>  CALL  F T Ck  ( DC  10  , V , °X  , CF  ) 

C 

MR  IT._(L°,iG) 

CALL  SRCr  1 (X,Y,PW,K,C,'J,V,PX,PF,FT) 

C 

CALi.  T °GK 

♦(A,0,PA,DO,rr,TX,TV,TH,PM,K,PX,°T, 

♦ YI,YF,NY,MY,IYI,IYF,NYI,TI,Tf,NT,MT,ITI,ITF,NTI) 

CALL  T P r>  K 

♦(A,J,PA,PO,TT,TX,TV,TH,PM,U,Px,PT, 

♦ Yi, YF, NY, MY, lYI.IYF, NY  I, TI.TF, NT, MT.ITI.ITF, NT I» 

C 

I F M F I . NF . t I CALL  F T 0 W ( OF TK , K , 3X , P T ( 

IFMFT.NL. 1)  CALL  FTCW(PFTU,l),  px.FT) 

c 

Rt TORN 
END 

C hhiihh  * it » <i  mu  h ii  n a n n n n ii  hh  n n ii « n n n n n tin  n it  im  ii  n h h h h im  n ii  n *«  * # a hh  n urn  # » h n an  h « * * # n 
C oUBRJUTINiS 

c H UK  tit  Hit  It  H ft  ll  till  n It  II  it  i)  II  II  II  II  ‘I  Iffl  il  I' ll  H II  ll  II  ll  II  II  /III  IIH  H Him  II  It  II  *11  * H III!  ll  IIH*  II  0 H » HII*  * HD  »«  0 

oUOPOU  T 1 Nc.  SRC 

CONS  T A N T S 


COMMON /F;./Om  , L.  F - , 'If  , c t . OE 
Com  io  n/f.  c./j-’i  i z i- 1, mi  • 1. 1 . o I 
COMMON/Fi  '/'.I pa  .ERA  .IJA  ,t.A.OA 
C0MM0N/Fli./JRt,£Rw,U^,.  3 , C2 
COMMON/FI^/ JRSm.  PS,  US,  r'S,  OS 


COMMON /Flo/  IR: 


.»■  l ,ni,E  i,oi 


COMMON /Fl 7 

/ J W , ***,'!  N , 

Un  = DF ( J=>_ ) 

‘■.t  ::  F !£►  L 1 

JI  = JFC.1-I  ) 

»•"  I = r F { l S.  I ) 

JA  = JF  ( KA  ) 

5 ; A * £ F ( : A ) 

U2aUF «UP2) 

*■-.  L=tF  (r.p  ?i 

US  = 'JF('JRS) 

• • C-LMr-J 

!J„ = JF (U- . ) 

< 3 1 =L  F (£P 1 ) 

'JW=  J F ( ' J R M I 

* £ W = £ F ( t.P  W ) 
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ooooooo  O O O DO 


RtTURN 

_N0 

SUB  -<  J JT  I N j-vSC 
C 

CONj TANT S 
INT.GlR  CF.LP.pQ 


comjti/J^/^.L^ , 3 1 
COMMON/Ci 7/APS, ?PS, CPS 
C(JMMJi,/C*o/Af\x,RAX,Ct-X 
C0NM0N/Ci7/ARS»pRS*CFS 
C0MM0n/C2L/pGPj,PSRx  • PG  r S 
COMMON/C.1.,/  Jo PS  , DC 3S,  i<  PS  , N PS 

CJM  1 CM  / C 2 0 / J S -x  X , D C A X , D K >•  x , M A x . N >.  x 
COMMON/Cc  ; / °i»5  , DCRS.  OK  PS,  MAS  , NR  3 
COMMON /C2A/T  S°S«  T ~PS « 1 OPS.TPPS,  TOPS 
COM-10U/C2t/T5’xX,TARX,TOi*X,TPI-X.TCEX 
COMMOfJ/Ccr  / T jp^j,  T - R Li  • TrcS  * T G » T C R S 


1 FORMAT  (Auli 

j i 

ii  format 

' AL  ° H A 

p3E  STRIKE 

", 1PG20. i:/ 

♦ '•  '* , 

' JlTA 

PR E S T ",  IK_ 

", IPSlC. 10/ 

* G A MM A 

P **  L S T R IKE 

", 1PS2G . 1 J / 

♦ ,Bj"  , 

•pLlSL 

PlAK 

p- c strike 

",  1 pG2l .10/ 

VULSf 

MAXIMUM 

p-L STRIKE 

", 1PS2C . 10  / 

♦ "j". 

*PULS£ 

G'Jt  SS 

prestrik- 

, xPGluiIC ) 

12  format (M:Mf 

*AL°HA 

RlTURN 

", 1PG2C . I3  / 

♦ " ", 

*Bc  TA 

RE  TURN 

", 1DS2C . 10/ 

t " ", 

•gamma 

RE  TURN 

",  lPr.2u.13/ 

♦ "j"* 

*°ULSE 

p _ A K 

Rl  TURN 

",  1 DG2 ( . 1J  / 

♦ " ", 

1 JULSL 

MA  X 1MUM 

Rt  TURN 

", 1PS2Q.1G/ 

* "G4\ 

•3ulSl 

GU*  SS 

RETURN 

", 1PS2G. 10) 

13  FORMAT  ("G". 

* A L pm  A 

RE  STRIKE 

",  1 PG2  u • i.0  / 

f " ", 

*dt  TA 

RE  S TR I KF, 

", 1PS20.1C/ 

4-  M " , 

'Gi  MUS 

R.STRIKl 

", 1PG2C . 1 C / 

♦ “ J '*  * 

’PUL  SE 

°CAK 

RlSTPIKE 

", 1PG2G. 11/ 

♦ ••  ••  ^ 

•pGLSr 

MAXIMUM 

R - STRIKE 

",  IPGrG . 10/ 

♦ M3M'. 

• pul  s _ 

GU£  SS 

Rl  STR.IKl 

", 1PG20. 1C ) 

IF ( NPJ( N: . ^1  CAlL  S *■  C S 

♦ (TRPS.TCPS, T PPG, T CPS, DC  PS, OKFS.MPS, PSPS. PM PS, PSPS, APS, BPS, CPS) 
I F ( NRx  . N;_  . J ) CALL  SRCG 

M TRRX  , TDRX,  TPfxX  , 1 CRX  . OC*xX  , OKR.  X,  MR  X,  PSRX  , PMRX  , PGRX,  ARX  , BRX,  CRX) 
IMNRS.NE.G)  CAlL  SPCS 

♦ C TRRS. TORS. TPRS, TCRS.GCRS , DKPS, MR S , p SR’S . pmR j , PGRS. A *S , 8RS, CRS > 
RE  A 0 ( C R , i)  MlC  tlF(MLC.tO.G)  RETURN 

I P ( NPS • Ni • o)  WRIT' ( L P , 1 X ) AcS,BFS,CPS,PoPS,PMPS,PGPS 
IF  (NRX.iVL  #0)  WR1TE(LC,12)  A R X , OR  X , C° X , PS R X , P MR  X , PGR  X 
IF ( NRS, N_ . L ) WRI TE (LP. 12)  A 3 S , UR S , CRS , P SR S , P MRS , PGR S 

RETURN 

lNU 


SUBROJTIN.  SRCS(TR.TO,TP,TC,PC,DK,M.  , CU , CC , C G , A , 0 , C ) 
CONSTANTS 

JOUBLE  AT, .MAT 

OOUJll  p-x-CISION  p T , f MR  T A-14 


oooo  o oo  oo  oooo  o oo  o no 


w F 


R£A,_  NRF 

c 

C0MMCN/C13/T SX, TFX, TOXtTPX, TCX 
CGMMON/Ci^/AO, 30, *0 
C 

_ X T £ RnA  L COF 
£XT£RNAl  C03F 
C 

r r<x - r r rrox  = ro  f t - to-tp 

c 

40  = - AlOG ( DC) / T K fl  = A o 
30=AJ*CG  «3sNRF(60,C0F,G00F» 

C = 0 JIF(Mi  » C = - A L 00 ( OK  ) / T C 

C 

AT=A*TR  JlMAT-JIXP (-AT ) 

3 T = 3 * T R iz M 3T  = 0- X 3 ( -9T > 

C 

cc=co*  (£MAT--:M3n 

c 

re  r jrn 

£ NO 

0¥¥»*¥***¥¥¥¥¥*, 

SUBROUTINE  3RlS(P,Q,PS) 

PUlSc 

INr.GcR  °S 
COMPLEX  p , Q 

DIMENSION  P<  Pit  , 0 (°S  > 

CALL  SRL<(P,PSI 
CAll  3 R L C ( 0 , P 3 ) 

it  1 J <i  J 
1 NO 

311301  T I v G-.<  (x  ,PS1 
3UL  3=: 

INT. j-X  P _> 

Rfc  A L < r 

CCMPllX  < 

OINl  J3IUN  k<  p^» 


Clm ion/c- i/r  u, wf 
COMMON/CIc  XNK , NC 
CCM'"IOi</  "I  5 /i  , 1 < , v Y , FW,  MS 

TI  = ^ f.TF  = TM  INT=N<  tyT=Mr 

00  tUu  XT-*,NT,i 

r=AF(ri,TF,,Nr.:r,'-n  ?« d t i =kf ( t» 

i Co  cofjrxNU- 

Rl  TURN 
t N 0 


SUS”  Jl/T  IN.  ^PLLir.o-) 
P U L ->  l 
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L 


0 0 M P l . X CtCF 

c 

JIMcWSION  C(JS> 

C 

COM1 JN/Cil/T  1,  WN 
COMMON/O-C/'K.NIC 
COM  ■UN/0i5/'1£  , IX, vy.  If, MT 
C 

><1  = 0 5><F=hM  S\M  = NC  5MW  = MF 

C 

)0  mi  I«  - 1 > Mw  , * 

4=AF (WI ,WF,  <W,  Iw,  M»<>  TC  ( I W)  =CF  { W» 
lOu  CONTINU- 
C 

r jxn 
£ NO 

Q*********************#**#**»»**»*X¥****** 

SUB-XOUUNt  S~?^<  £ ,3C,  ZC,  PE  ,PF) 

C 

C S T fXO  < w 

C 


i NT  ij _ 

u_ 

*1 

a 

COM^L-X 

COM  3lc  X 

PC,ZC 

COM  D X 

oc 

c 

01  “1C  Nil  ON  c ( ar  , 3 F 1 

01 M£  Ni ION  PC(PF) ,ZC(cF> 

C 

COMMON/G-l  /SI , 3F,KS 
COMMON/Oii /XI , XF , NX 
COMMON/Jlc/.-II  , 4 F , N H 
COMMON  / 015/"tj«  MX, MV,  M w . H T 
COMMON/0££/L)I,JF,NU,mU 
COMMJN/0££/VI,VF,NV,MV 

commjn/o£i./x:,zf,nz,mz 

c 

ui=x i tvi-si  tzi=wr 

UF  = XF  «VF=oF  {ZF=WF 
NU=  NX  SN»/  = N_>  JNZ  = NW 
1U  = MX  XMV=Mj  JMZ=«>< 

c 

00  x C L IW— i,NW,i 

W=AF(rtI,WF,NW,IW,MW) 

CALL  iRy(W,JC(IW)  ,ZC  (IWI  * 

JO  LOb  15  = 1* NS,  1 

C = AF<sl,SF,iJS.  IS,  M3) 

1 (1S*I<0=CF(  3 , V4  > 
iOO  CONTIN'JL 

C 

•XiT-J-xN 

£N0 

SUfliXOU  T I Nt  5>XY(W,3C»ZCI 

c 

C °R0  ° AO A f I ON  CONiTANT/CHARACriFISTIC 

c 

c p a r a m£  t = p p-e=iO 

c 

I NTc.G£  R 

INT^G£R  CP.LP.PO  .. 

COMPLtX  WN.WNF  A*lt> 

COMPLEX  NI.WIF 


L 


Is 


COMPlc-X  PKG,  ZKG 
COMPLEX  RTS 
COMPLEX  T 3C 
COMPLEX  ZX.ZL.ZT 
COMPLEX  YS.ZS 
COMPLEX  KI.KIF 
COMPLEX  < C t <CF 
COMPLEX  PC.ZC 
COMPLEX  OR 

logical  test 

c 

DI ME  NS  ION  RTS(lCi) 

C 

COMMON/Ai0/CR,LD»C0 
COMMON /Bl./IMX,CMX 
COMMON /927/MR 
COMMON/B3G/MX 
COMMON/Dti/XX.WX 
COMMON /El3/RAM*KAP 
COMMON/cl  -t/.RS 

COMMON/Fil/iJR£,  EPE.UE  ,EE,  OE 
COMMON/Hli/<C,KI 
C 

EXTcRNAL  0 KF 
C 

DATA  PR/lO/ 

OATA  k.  2/1.*»a>.c:135  6/ 

DATA  N°R/C/ 

DA  7 a nr  /!/ 

OATA  TEST/. FALSE ./ 


c / 

10  FORMAT  ("u",”FRtOUENCY  ",  1PG20.1U  I / 

11  FORMAT  < ”u  " t “FRE  QOt  NC  Y ",  lPG2w.ll  /( 

♦ " "."WAVE  NUMBER  ", 2(1PG20.10) / 

♦ " ”, "waVE  IMPEDANCE  ", 2(1°G2u.1C»  / 

♦ " " , " Pf  0 c A G AT  I ON  CONSTANT  ",2(*PG2G.ll)/ 

♦ " "."PROPAGATION  IMP. DANCE  "»2(lpG2o.lL) / 

♦ " "."PtOPAGAT ICF  CCNcTANT  " , 2 ( IPG  2 C . .. 0 » / 

♦ " ","AOMI T TI VI r T ",  2 ( iPG  20 . 1 C ) / 

♦ ..  "."impeDIVI  TV  ", 2(1  pG  20.11)  / 

♦ " "."CHARACTERISTIC  CONSTANT  ", t ( 1°G 20 . 1 C I / 

♦ " "."CHAPACT.RISTIC  IMPaDANCE  " , 2 ( 1PG  2 0 . 1 w ) ) 

12  FORMAT  {"  "."PROPAGATION  CON_.TANT  " . 2 ( 1 “G  2 u . 1 0 ) ) 

13  FORMAT ("  ","LMAE ACTtRISTC  IMPEDANCE  ",  2 ( 1 °G 2 0 . 10 ) ) 

C 

JC-Jt  F.C=2_  iOC=0£ 

C 


WN=W  <F  (w.UC.EC,  30  S F<G  = Wf. /P2 
Wl=WlF(W,UC._C.OCI  F ZkG  = w I 

c 

IFIMR.ED.UI  GO  TO  100  TIF(w.EO.O)  GO  TO  lwC 
C 

wx  = w 

RTS(l) =C 

CAL.  S?M(np-:,NR,cT3,im»,lMX,.mx,PnF,TEST) 

TPC=RTS(x) 

CALC  jRZ(W,TPC,Zx.Zl.ZT) 

Y 3s  i/ZT 
ZS*ZX»Zl 

<Cs<CE CYj.Zj*  A-17 

< I s< I F (YS.ZS) 

c 

I F ( i x • N ,1.  ) W-1T:  ( L ' , 1 1 ) /..WN.-s1 1,c<j,7<G,T  PC  ,YS,Z-5*XC,XI 


/ 

/ 


I 


c 


PC  = <C 
2C  = < I 


t,'=CAuS(  "K  > ilf  (CC.?-  (f-KO)  ,N_  ,C  1 £ 
IF  ( _-R.LT  ,LS>  f_  TM?N 


F=CAJS  (QP/DKb> 


CONTiNU- 

°C=a<0  Ew-IT_(l.c.1c)  FC 
/C=/<o  fW-ir_ (LF,i3>  7C 

•?lT  JPN 
ENO 


• *¥**♦*** 

SUBPOuTI 

N£  P<F( TPC.DF Tt 

P-GPA3A I 

ion  Constant 

lNT£Gc.P 

os  0 

C 0 M p l _ X 

CI 

COM  = Ll  X 

- MW,  ^ MI  . L M£ 

CO  M °L  c X 

WNH,WNI«WN> 

COMPCt-  < 

A D CN  « ~ °C  I , APCt 

C 0 M P L _ X 

W.F, A ?CF 

C 0 M p c.  X 

2 I W , x I r 

COMPll  X 

r °c 

COMpL.X 

A - G 

COMPLEX 

OF  ,Nc,HFM,  Mc  p 

COMPLEX 

J L W , 0 1 >,  W 

COM  3l  ^ X 

Jb iw , ji i h , rG I w. 

COMPLY  X 

JbI^«JiIS«YC  I , 

COMPLEX 

HitCfHiES 

COMPLEX 

NJfOJ, WJ 

COMPLlX 

N M , 0 M , rf  H 

COMPLt  X 

T 

COMM JN/ J2./X  , * 

COMMON/-.  * - /■>■  X . PH 
COMMON /£!*«/<$ 

COMMON/‘ri*/K£.E.*r.,l:..LS,Oc 

COMMON/Fi^/'jKi  , CM  I ,IJI  ,L  I * Cl 

COMMON/Fi3/0'XW,:.PW,l|W,EP,CW 

COMMON/  UC  /WJ,  WH 

COMMON/IIa/JOWH.JIWW 

COMMON / li2/Jl)IW,JIIWfYGIW,YiIW 

COMMON/ II 3/ JO IS*JiIS.rOIS,YlI3 

COMMON/ I^H/rtUt  3, HIES 

COMMON/Ii5/_MW,£MI,tMF 

COMMON/ II 6/A PC H, A °CI ,ACC£ 

CI-CMPL<((u«Oti.C) 

emw  = ow*-ci * w* :w 
a.mi  = oi»ci*w*-:i 
-Mt  = Ot.»CI*H*  Li. 

WNW  = WNF(W,UW,:_W,OV>)  ?APCW=ACCF(  T CC , WNW ) 
WNI  = HNF{W,llI,c.I,Oll  *r.DCI  = APCr<TPC.WNII 
WN£  = WNFM,U  ,.CfO-)  t A PC.  - A PCF  (T°C,WNEI 


PIw=  (l  1W/t  MI ) * ( APCI/APCW) 


f 


RIt  = (tMi/t  MI  ) * ( APCI/AFCt  » 


ORJ  = j SARG=APCI»F  W 

CALL  FNHF  ( CRO,  AR  G , BF,  NF,HFM,HFP> 

JGIW^OF  *rGlW=NF 

ORJ  = i SARG=APCI*RW 

CALL  3NHF  ( 0^0,  ARG  , SF  , NF.HFM.Hr  P> 
JiX  W = BF  |YiI W = NF 

JRD  = 0 *ARG=A3CW*RW 

C A L l.  BNHF (ORO, ARG, BF , NF,HFM,HFO) 
JG  WW  = OF 

JR  D = 1 J ARG  = A PC  W*  RW 

CALL  BNHF  ( OR  J , At  G ,BF  , NF  ,HFM,  HF°» 

J1 WW  = BF 

JRO  = J £ AR  G = A JC I * - S 

call  BNHF  <0RJ,ARG,-3F,NF,HFM,HFP> 

JCIS  = BF  t Y C I S =NC 

JRD - 1 SARG-A  DCI*F  G 

CALL  BNHF(ORJ.ARG.BF.NF,Hrn.HFP) 

J 1 1 S = 3 F 5Yi.IC  = MF 

OP.O  = 0 £ARC=ApCl*cS 

CALL  dNMF <OR0,ARG,9F,fF,HFM,HCP> 
BG£  S = HFP 


JPD  = 1 J A R G = A ° C , ' * l 5 

CALL  3 N M F ( 0 • },A=G,3F,NF,HFM,WFP) 

-Uc.S  = HF° 

NJ  = R IW*  J1  WW*  jo  I W-JG  WW*  J*  X W JNH=LI‘_*H1  i • J „ I S -H  ? L a*  J 1 1 S 
oj=riw*jj.lr*yoIr-j:mw»yiik  soHr.-i.-.  5*r:x  -hcf c*yii<5 
WJ=NJ/ X J tWH=NH/ JU 

JET  = W J - R H 


RETURN 

ENG 

SUtJP  JUT  XNl  jR/(^,TPC./S,/L,/T) 
XPP^CANC; S 


COHJllX 

Cl 

C JM  Jl_n  < 

T PC 

COHPL_  X 

cMW,  - MI,  -_M-. 

COM^Lt  X 

A F cw . A PC  I . A»Ci 

C 0 M J 1.  l < 

JC WW, JiWW 

COM1 

Ju iw * Ji xw * yC i w , y i : w 

COMPLEX 

JlIS.  JilS.U  1 S , Y 1 1 -j 

COMPLl X 

N L S , H 1 t S 

C 0 m p L --  x 

WJ  ,WH 

COMPLl  X 

°c 

COM°Ll X 

CJ.OY, JU 

COMPLl  X 

0- , JS» JW 

CGMjLl X 

Rf.-'T,  J. 

COM  °L i X 

cc 

COMPl.  X 

cs»clp.cl%ct 

COMPLl  X 

C JY , CM 

COMPL  . x 

' 

4 •> 

o o o 


c 


c 

c 

c 

c 


c 


c 


c 

c 


c 

c 


c 

c 


c 


COMJ|_x  7LP.Zl1.ZL 

com->l_x  zij./n,/T 

C0M‘1JN/^„/vh,.-k 

COMM JN/rii/J  »eP£.UC .Cl  »0t 
COM'ION/F’.i/JP. » « C k I . 1 1 1 . r i.OI 
COM  '10,l/I»i/WJ,  WH 
COM-IOil/I**/Jt'WW,  JlrttJ 
COMMJN/I_£/JoIW,J*I.V,V0IW,riI'W 
COM-IOM/IiL/JuIi.JilS.YOIS.XIlS 

common/ i:  *-/hoi i.w;LC 

CUM.UN/  I*  S/t  MW  , LMI  ,£  Mf 

COMMON /I»c/APCWfAFCI. AFC f 
data  =>i /i  . 5 42  tc./ 

CI  = L -IPLX  (U.  1 . ^ . j 1 


PC  = WJ 

0J=(JuI3-JljIW)/JlWW 
Dv  = < y u is- y Cl w > / j; lw 
OH " H J t S/J1*N 

JP- Jjnw/Jiww 

os= jl i s-uc*Yu :s 

OW=JlIW-^C*Yi  lw 

R D = ♦ (lML/lMI  ) * (iPCW/Jt>cl)  * (J1WW/JM 
RQ=-<cMW/£MI»*(AFCW/APCI)*(J1«W/OW)»PC 

Ri.  = «’(£.Mk./_MI)*(flPCI/floCL)*(HPCW/aFC.)*(JiWW/Mdi>*(as/0W» 


CC=L*°I*RW*LMH*aPCW 


cs=<apc'w»apcw>/cc 

CLP=  (CI*w*'JI*eMI  ) /CC 
CLl}=(CI*W*U£*£ME)/CC 
CT=;/CC 

CJY=RF*OJ*RQ*OY 
CH  =R(.*OH 

ZLP=-OLF*CJY  fZTF=-CT»CJY  S7b=CS*r!R 
ZLQ-*CL1*CH  /7T0-*CI*CH 
ZL  = ZL D+  Z L .3  f ZT  = ’T=>*ZTO 

RLT'JRN 
t.  N 0 


SUORJJTINi.  S^CC ([  *PC.ZC.IT,MP,P£  , p P 1 

COAXIAL 

C A 9L  e 

1 N T £ G £ R 

°£  , °F 

INI- GLR 

CF.wP.PQ 

COM°L£X 

W NX . WNF 

COMPl.X 

WIX, WIF 

COMPLEX 

£,IT 

COMPLt.  X 

PC.ZC 

COMPLEX 

MR 

OIMeNSlON  c(p£ , PF) . IT <P£ , PF>  A.20 

Cl  Me.  N$  ION  PC(PF)  ,ZC<PF) 

DIMENSION  MP(DF) 

o o o 


c 


C0MM0N/AiG/CR,L°,°3 
COHMON/03C/MC 
COMMON/Oli/XI,XF,NX 
COMMON/Ol2/YI,YF,NY 
C0MM0N/0.,3/WI,wF,NW 
C0MM0N/015/WS.  MX,  MY,  MW.MT 
C0MM0N/D22/U1 ,UF ,NU,  MU 
C0MM0N/D23/VI , VF ,NV,MV 
CO MMON / 02 h/7 I, ZF.NZ.MZ 
C0MM0N/F16/UR.,,EP1,U1,E1,01 
C 

10  FORMAT  ( "0  " * "FR£  QE<l  NC  Y •*,  1PG20.10  ) 

12  FORMAT  (“  “ , "PROP  A G A T ION  CONSTANT  " , 2 < IPG  2 G . 1 C > > 

13  FORMATC  ","CHARACTERISTC  I M F£  DANCE  " , 2 ( 1 PG  2 l . 1 ti ) > 


C 


c 


102 

101 


UC  = Ui  SEC=ti  EOC=01 

U I - Y I f V I =X  X tZI=Wl 
UF  = YF  *VF=XF  £ZF=WF 
NU  = NY  i N V - NX  £NZ  = NW 
MU  = MY  JM  V = MX  £M/=MW 


OO  101  IW=1,NW,1 

W=AF(Wl,WF,NW,IW,MWJ  t I F ( MC • NE , 0 • 

WNX=WNF< W , UC. £ C, OC)  JPC(IWI=WNX  JIF(MC.NE.G) 
WIX=WIF  (W,UC,lC,  OC)  «ZC(IW>=NIX  ?IF  ( MC.  NIL.  0) 
00  lOZ  I X =1 , NX  , 1 
S<IX,IW>=IT<IX,IN>*MR(IW> 

CONTINUE 
CONT 1 NUE 


WRI  T£  (L  °i 1 0 ) 
WRITE (LP,12) 
wpi rr  (lp,i3) 


RETURN 
£ NO 


C 


SUBROUTINE  SRTl(X,Y,ow,C,V,..,cC,ZC,Pt,PX,°F> 


W 

WNX 

WIX 


TRA  N SMISSION  LINE 


INTEGER 
INTEGER 
COMPLEX 
COM^Lc.  X 
COMPLEX 
COMPLtX 
COMPLc  X 
COMPLE  X 
COMpLlX 
COMPLEX 
COM°ll X 
LOGICAL 


Pt , °X, PF 
PW 

C,  V 

L 

PC,ZC 
WNX, WIX 
ZMX, ZPX 
zlm, zl° 

CNF , YNF 
Y 

OF 

test 


DIMENSION  C ( Px , PF ) , V ( PX , PF) 
0 1 Mt_  NS  I ON  E(Pl.°F) 

3 1 Ml  NS  I ON  PC ( PF  ) , ZC < =F ) 
OIMlNSION  X(PW>,Y(PW) 

COMMON/R15/N3,NI 
COMMON /P3 I /MM 
COMMON/ 02i/XX,ZX 
COMMON/D2l/X 1 ,XF, NX, MX 
C0MM0N/02E/YI , YF ,NY, MY 
C0MM0N/02-,/ZI,ZF,NZ.M7 
SOM  ION/GII/’.m, ZCc 


A-21 


o o o 


COMMJN/Mi  7/RNX , W I X 
COMMON/Hi-//  lx* 7PX 

-XTIRNAl  G*F,V<F 

DATA  1M/U/ 

ca  Ll  ^raf  o,  p* , y i . yf, r.r . my) 
do  iou  if=*,nz,i 

ZX=AF(ZI,  7F»  NZ  » 1°  » MZ I 
RNX=PC(i°> 

RIX=ZC (IF) 

ZMX=ZlH  SIF(MM.Nt.O)  ZMX  = W I X 
ZPX-ZLP  IiF( MM. NE. J > ZCX=WIX 

call  sro< im.ip.e ,p^,pf,y,pw> 

00  100  IX  =1  * MX  * 1 
X X = A F (XI,  XF.NX,  IX, MX) 

TEST  = (MR.Eil.C).OF.(MM.EQ.i) 

IF  ( TEST)  C(IX,IF)=QF(YI,YF,NY,  TY,NO,CKF,X,V,ow) 
TESTMMR.t  Q.u)  .OR. (MM. £0.21 

IF ( TEST)  V (IX  , IP)  =TF ( YI  ,YF ,NY,MY ,NQ, V<F,X, Y, PW) 
100  CONTINJl 

s 

RETURN 
£ NO 

SU9R0U  T I NL  CRT ( MR , o F ) 

TRANSFER 

INTlOcR  pf 

INTEGER  CR*LP.PQ 

RE  A c MSI,MTA,MTO,MTS.MTT 

COMPLEX  SI.SIF 

COM  °Lt  X TA,TD,rS*TF 

COMPLEX  TT 

CO  M °Li.  X MC 

DIMENSION  MP(PF) 

C0MM0N/A1u/CR,lP, 3T 
COMMON/Oi 3/R  I * WF  * NX 
COMMON/ 01 5 /MS* MX .MY.MW.MT 
COMMON/E 1Z/RSM,RS° 

C0MM0N/E1 S/RAM, « A ° 

COMMON/F1J/URA, EPA.UA.EA.OA 
C0MM0N/F*.-./JR<l,EPZ,UZ,EE,02 
COMMON/F 15 /URS.t RS.US.ES,  OS 


FORMAT  ("0",  “ INDEX”,  IX,’*  FRE QUE NCY” , IX 

► " ARMOR", IX,"  SP A CER” , 1 X »” 

► " TOTAL" 

FORMAT  (”  ",I10,1X, 6(1  PC. It  . 5,1  X)  ) 

WRITE(L°,lO) 

DO  100  I W S1 , NW , 1 

R=AF(RI,HF,NR, IW.MR)  SIF(W.EQ.O)  GO  TO  iJO 
SI=SIF(RAP,R,UA,FA,OA)  iMS I =CA eS ( S I I 

TA=TF(RAM,RAP,W,UA,EA,OA)  tMTA=CAES(TA) 
TO=TF(RSP,P.AM,W,U2,E2,OE)  *MT D=CA Fa ( TO) 
TS=TF(RSM,RSP,W,US,ES,OS)  t*T S=CAPS ( TS)  . .. 
TT=SI*TA»TO»TJ  fMC(IW»=TT  tMTT=CA6S( TT) 

RPITl (LP.il)  IW.W.MSI.MTA.MTO.MTa.MTT 


SURFACE", IX , 
SHIELD", IX, 


I 


ooooo  o o o o ooo 


0 N 1 I NU  ; 


ICO 

c 

•<£  TURN 
£ NO 

C* »*••*»*»*»*.*.****»♦***» 

SUBROUTING  SkO(IX,IY,W,PX,PY»7.PZ> 

DATA 

INTGGER  PZ 

CGM°l_ x w , z 

DIM.N3I0N  HlJX,oy| 

OHM-NSION  Z(JZ) 

IF(IX.tO.o)  r.o  TO  1C1 
IF ( I Y . £ T. l > ,o  TO  l02 

1C1  CO  NT  I NU  - 
JO  liu 

Z<IS)  =W ( I'  , I Y) 

i i u co  n r i nu  _ 

R c.  TURN 

102  CONTINUE 

JO  12G  IS=>1.PY,1 

Z(lS>=rf<lxtI5l 
12C  CONTINUE 
Re.  TURN 

£ NO 


SU3R0UT  INC  GPCFT(X,Y,PW,K,C,Ut\/,PX,PF,PT) 
FOUJIc.P  TpANSFJCM  (CONTINUOUS) 

INT_  3£R 

I n r g g _ r 

COM  Ji.:  X 
COMPL- x 
Complex 
CO M 3L  c X 
COMPLEX 

c 

JIM  Nil  ON  K ( 3<  , c T > , C ( FX , PF  ) 

MM  N _>  I 0 N U(=x,rT),Z<FX,PF) 

JI  Me.  NjION  X ( PJ>  , Y (OWI 

c 

jommon/  ?ic'/  n . mi 

COMMON/Di 2/YI , Y F , NY 

I J M -1  j M / iZ/rt  1 • X F i * X 

COMMON/  .‘./TI  i Tr  ,KT 
C G M M J N / .5 1 • / MS.MX,FY,M«,MT 
; o m m j t / ; 1 1 / r x,xx 

c 

£XT_n!NAi_  IN’- 

c 

)A  T A il/i/ 

C 

CACu  j P A * (XfPWfWItJFtNWfMW) 

JJ  iC-  IY=1,NY,1 

<X=«F ( YI , YF, NY , I V,MV ) 


F*«PF,PT 

K » C 

U.  Y 

IKF 

Y 

IF 


c 


A-23 


CALL  s’0(IY.If'1tCtFX«PF«Y.PW) 

30  ill.  IT  = .,ljr,I 
TX=4F<  TI , TF, NT , IT ,MT ) 

<(IY«iT)=IF(0I.wF*Nr'(,RN,Ni»lK’F,X»Y*oW) 
iCi  CONT IN J£ 

CALc  jF'{  1Y,  IM, V,FX,OF,V,f«| 

00  * u — I T = 1 * N T • I 

rx  = «F ( T I, TF.NT , IT ,MT  » 

J(IY  * I T>  = I F ( R I » R F iNW|Kln,Nl,  IKF  < X « Y , P R J 
lu<L  CO  NT  I N'J_ 
lOu  CONTINU. 

C 

RETURN 

ENO 

£********»**********»******♦*♦♦»*****♦*****#*******#**¥**¥****¥***♦**#** 

SU'IROJTIN-  GRAF  (A  . P F f/ItXF*NX*MX) 

C 

C ARRAY 

C ' 

INTtGc.R  PF 
C 

01 M_  No  ION  A ( Jr  ) 

C 

30  10l>  IX=i.,NX,l 

X0=AF(XI, XF, NX, IX, MX)  |A(IX»=XO 
1C J CO NT I NUb 
C 

R_  T JRN 
ENO 

C inti>  it  »*»  v mi  mnnHi  » >nt  mum  v ini  n it  <t  tin « it  oh  it » »»  # /« # # o « # tv  it  n » # * » n » if  1 1 it  I »*  * # # » # # 

C FUNCTIONS 

REAL  FUNCTION  UFIURI 
C 

C PERMEAOIlITY 

c 

DATA  UO  / 1 l.  • 5 bo  4E  ” C 7 / 

C 

U=UO*UR  I U F = U 
C 

RETURN 

C.NQ 

c..*. «•••••*• 

RE  Ac  FUNCTION  £F(ER) 

C 

C PERMITTIVITY 

C 

DATA  1 0/ 3 • 85  *♦  £ E • 1 2 / 

C 

L = EU*tK  tE  F = t 
C 

RE  TURN 
c.  NO 

RE  A L FUNCTION  <F(T) 

C 

C CURRtNT 

C 

INTEGlR  PI 

c 

COMMON /ClD/APb*  9 F i * C PS 

C0MM0N/Ci6/ARX,PRX,CRX  A-24 

COMMON /Cl7/ARG»PPS«CSS 

COMMON /C£l/PiPS*OCPStO<P5>*MPS»NPS 


c 


c 


c 

c 


C0M:M0N/Ct;,_/-:,3RX.0CR>  , OKI-  X .MRX  «NRX 
COMMON/C23/pSmS,  OCRS,  U<F3,MMS,  NR3 
COMIOr-l/CiH/TiPS.T-JPS.TOPS.rFPStTCPS 
COM.iGN/C2  5/T3-'X  , T RRX  , TO^X  , TPF  A , I CSX 
C0MM0N/C~6/  TSRS  , TRRi,,  TDRSf  T PR3,  TC^S 

data 

DAT  a NS/  1/ 

o a r a ps/*i/ 

PPS=  ST  F (PSPS,  T 3 °S,  T°PS,  T , A°_>,  3 Pj,  0 = 5 , PS  , NP3I 
PRX=STF<P3RX,TSRXfTPRX,T,AXX,TRX,rRx,NS,NPX> 
PR  s = 3TF(P3RS,TSRS,T°PS,T.ARS,0RS,CR;>, MS, NkS) 


PT :PP_ 


|KF:pT 


Rt  TJmN 
£ NO 


COMPil X FUNCTION  CF(W) 

C 

C CURRENT 

c 

INT^GiR  Fj 
COMPLLX  pps.ppx.pfs 
COMPLEX  S f F 

complsx  pt 

c 

COMMON/C15/A°3»cF  :,C°£ 

COMMON/C1 C/A  RX, PAx , CP* 

COMMON /Cl 7/AFS, Sf~ , C»S 

COMMON/Ct »/PpPS, OCPS, CK  Pj ,MFS, Nc^ 

C0MM0N/C22/O 3 A a , CCRX , D*>  X ,MRX , NA  X 
COMMOfo/C2  :•/-■>-  S , £ C R >*C  K f I . MR  r , NF  ; 
COMMO  l/Cn»/TjFC*T“J;j(7FFSiIocS>TCcj 
COMMON/CLp/Tj-X,  t - - x , T x , T x,  r c - < 
C IHMQH/C2  • / r , j , i Cf  J.TPRS. TCI 

c 

OA  T 5 M3  /- . / 

OAT  A N 3 / i / 

JATA  OS / F 1 / 


° P S = j F f ( p . r> , i ;fs.  Top  ),W,AoS.'JP_,,rP3,FS,N°S 
p-  X - 3F  F t p ■ R < , T J x . T°F  * ♦ W . Am  X , P-<x  , C'  * , NS  . UP. X 
P"S=,jFF  (P  )“  )i  T »-'3*r-Pr  (I-i/I-'ji1-;  3»C*P,'13,NRS 


c 


r--=T3*I3*TP 

ttstu-ts  ?iu:J3*,*r*rr  m-i 

S UM=3U'1+  r.nl 
loo  COMUN.I 
o T F = iU’M 

C 

■ c TJin 
- NO 


oOM^ol  X 

F L'N  J T i 3 N CFF(o5,Xi, 

3 T n 0 < . 

< 300  )L:  . > £ONL  NT  III  ) 

I N T r.  3 _ 3 

3 ON 

COMJl.X 

C 1 

COM°_l X 

U£ iC  U vC 

COM3 !_  •_  X 

A S , A T 

CCMP^i X 

Cut.  T 

CO  M Jc : X 

Tc.*> 1 ♦ : JM 

C;=CivP L X ( L . 0 , * . O » 

c 

oFF:  J 3IF  < N->.  •_  J.  G > ~ ET'JFN 

CA=AtCI«n 
0 l!  s *3  ♦ o I * W 
CC=-SGN*C*CI‘^ 

V 

C S - 13  j * < 1/CA-  1 /C-?  > 

«i  = -Ci*W*To  £Cr  = C£XP>fA<3» 

C 

SUM  = j 

30  1 G . I " x • ‘4  ^ tl  t IS  - 1 * “1 
AT-'Io'CJ’T’3  CT  . F^CiX?  (A  Tl 
suh=  ju  m f r f f "1 
1 u u CONTINUE 

3FF=CS*CT* OJM 
C 

f>-  r jfn 

ENO 

C»  * 

•XEAL  FUNCTION  (Pj,  rs,  r ,4, 

c 

C CIH^ENT 

0 

30  0 HLt  PRECISION  A T , £ MA  T 
OOUdL;.  PPFCIOIJN  IT  i £M1T 
C 

«E  AL  < 

C 

<S  = o TIc-cr.LT.T3>  3L  TURN 

C 

TT-T -f  , 

A T - A * 1 T f . NAT  = 3£XF (- 1 T) 

•)T-  i r t:  MOT  = 3*  Xp(-9T> 

* = Pj»  < £ 1A  T-  MOT)  i<S  = kt 

c 

Pt  rux  j 

ENO 

■Cr  Au  FUNC.1  100  OUT  (X> 

C 
C 
C 


,rl 


»#***¥•*#»*»+♦#*.**.*».  + .»«««*«  *»»#»..» 

TF|Fi4tliCt3  ">  N t NF  t 


CU'X-XeNT 


A-26 


oooo  oo  oo  oooo  oo  o o ooo 


DOUBLE  PRECISION  AT  * E MA  T 
00U  3Ll  PRl  C I 3 ION  RT.E.MRT 
C 

C0MM0N/C17/T  3, IR, TO,  TP,  TC 
C0MM0N/C1L/A , D, C 
C 

AT=A*TR  idMAr=OiXP(-AT) 
3T=X*TR  JlMBT=QEXF <-QT> 

C 

F=BT*zHBT-aT*£ia T tCQF=F 


RETURN 
£ NO 


REAL  FUNCTION  C00F(X> 


CURRENT  DERIVATIVE 

D0U3l£  PRlCISIOS  AT,^HAT 
DOUBLE  3RE  Cl  3 1 ON  0 T » E MR  T 

COMMON/Cl3/T3,TR,TO,TP,TC 
COMNON/Cl L/A , B,  C 

AT=A*TR  JlMAT=OEXc(-AT) 
3T=X*TR  Ic.M3T  = OlXd(~0T) 


FD=TR* (1-3T) *lMRT  SCOOF=FD 


RE  TURN 
ENO 


COMPLEX  FUNCTION  WNF(K,U,E,0) 


WAVE  NUMBER 

COMPLEX  Cl 
COM3Ll  X UC,c.C 
COMPLlX  wne,wn 

CI=CMPLX ( u .0 ,1. C ) 


UC  = U 

E C = E 1 1 F < W . NE  • (J  ) lC=E-C1* (O/W ) 
WN£=-W*W*UC*£C  IWN=CSOpT (WNE)  IwNF=WN 


RE  TURN 

ENO 


CON  3L  l X FUNCTION  WIF(W,U,F,0» 

WAVE  INPlCANCE 

COMPLEX  Cl 
COMPLEX  UC.lC 
COMPLlX  WI£,WI 
C 

CI=CMPlX (t . J I ) 

C 

UC  = 'J 

cC-£  UF(W.NE.O)  lC=F “Cl* <0/W ) 

c 

wi£  = *L'C/.c  5w:  = csnRT  <wii  ) ;wif  = wi 


oooo  oo  oooo  oo  oooo  oo  ooo 


c 


RF TURN 
END 

COM  aL^  X FUNCTION  KIF(Y,Zt 

CHARACTERISTIC  If't  DANCE 

COMPLEX  Y,Z 
COM^L:.X  <12, KI 

<I2=Z/Y  B<  I = Cb  Jp  T (KI Z ) {<IF  = <; 

RETURN 

END 

¥*¥¥*¥**¥***¥**¥¥¥**¥¥¥*¥*¥¥*¥***¥****¥****■¥¥¥***¥***¥¥******¥*¥***¥¥** 

COMPLEX  FUNCTION  <CF(Y,Z» 

PRO°AOA  T ION  CuNSTANT 

COMPLEX  Y,Z 
COMPLEX  KC2.KC 

KC2=Y*Z  B<C=CSQRT  (<C2 ) $KCF=<C 

RETURN 

END 

¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥*¥¥¥¥¥¥¥*¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥•¥¥¥¥¥¥¥¥¥ 

COMPLEX  FUNCTION  APCF(TPC,WNI 

AXIAL  PROPAGATION  CONSTANT 

COMPLEX  T PC » W N 
COM°LEX  A PC 

APC  = CSQRT  (TpC*TFC~WN*WN)  -,APCF  = AFC 

RETURN 

END 

¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥*¥¥**¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥¥ 

COMPLEX  FUNCTION  CF(X,W) 

DIPJLl. 


c 

c 

c 

c 

c 


REAL  <0 
COMPEL  X CO.CF 
COMPLEX  EO 
COMPLEX  E 

C0MM0N/E1 5/ZD, YT 

COM  MON/ FI  1/UREtEP.E  ,UE.  EE,  OE 

DATA  PI /3 • 1 *1592  65/ 

UC=U£  *EC=EE  *0C=0E 

<0=l/(2*PI*JCt  *C0=CF(W)  I£O=K0*C0 
R2  = X*X«-YT*YT«-Z0»Z0  *P=S0RT(P2) 
E=£0MX/{R*»3>  » BOF  = E 


RETURN 


C 


COMPLEX  FUNCTION  CIF (R.W.UtE.Ot 


A-28  , 


c 

C SURFACE  IMPEDANCE 

C 

I NTE  i,E  P 0 = 0 
COMPLEX  WU,WNF 
COM°LcX  Wl.NIF 
COMPLEX  AR& 

COMPLEX  MBF,MI,F,MMFM,MHFP 

complex  u, ii 
COM  PL  t.  X D.-t 
COMPLEX  31 

c 

3 A T A PI/3.iL„59263/ 

C 

rlN  = WNF(W,U»_«0) 

WI  = WIF{W,'J,£,0> 

C 

ORO  = J SAKG=NN*R 

CALL  M3NHF< JRO. ARO.MCF, MN^, MHFM, MMFP> 

IG=M3F 

c 

ORQ  = I SAPC-=WN*  R 

CALL  M3NHF(JP0, ARG.MPF, MNF, MHPM,MHFP) 

Ii=M3F 

C 

OR=iG/Il  S3I=(PI/(2*rl*F» >»0k  iiIF=3I 
C 

Xi TORN 
END 

q,,,,,,,,*,,,,,,,,,,,,,,,,,,,,,,*4*,»^*******^**^^****^^******^********* 

CO  M = cc  X FUNCTION  TFfPI,RO,P,U,E,0» 

C 

C TRANSFE5  = A T 10  3 

C 

I N T E G£  P 0 = 0 
COM°i__X  WN.PNF 
COMJccX  A G 

COMPLEX  M OF . MnF , MHFrt, MMF P 
C0M°Lt.X  UO.XlO 
COMPLEX  IuI.XlI 
C0MdL_X  U0,<-0 
COMPl.X  lil.Kil 

com  =l c x o^ r , a: r, k jr 

COMPL-X  ip.jp 

c 

DATA  JI/.-  . i ai>  >2  65/ 

c 

HN  = = ( X ,1,  , ^ , Ol 

c 

100=0  F < l 0 = 0 
I L I = j 5<ul=u 
I l'J  = J ?<*D=U 


ODD  o o o oo  ooo 


IoI=mqf  *<ci  = ;-mfm 

c 

JAPr-=WN*RI 

CALL  MJNHF  (0-l),  OF  G * M p F , MT E . FHF  M,  MHF  3 ) 
IiI=MJF  tKiI  = MHF  M 

c 

Jt;  T = I j 0 riF<  vl.NL.O)  OET=IUO*KiI*IlI*<CO 
TiFfRI.N  cut  wTK=IAJ*Ml-Ui*<l  J 
RJK  = 1 f.lF<  RI.NF  .0)  KJk  = io  I*C  1 + Ii  1*<C  X 

c 

I?=({<.*PI*KO»/WN)*(WIP/'DfcT» 

J\  = WJF/C-T  -<TF  = JP 

C 

RE TURN 
pNO 

r*. ♦****♦*#♦.*******.*.*.****»***♦*. *..*♦*♦•** 

COMPllX  FUNCTION  C<F (XS.XP, X.Y ,PW) 

KlRNAl  (CURRENT) 

INTcGlR  = W 
COMPLEX  Y 
COMPLEX  E X « _ OX 
COM °L  _ X OF , CGF 
COM  PL -X  <F 

OIMENGION  X(PR),Y(PW) 

COMMJN/P1 L/NO , MO , MC 
COMMJN/3._5/Xl,XF,NX»MX 

CALL  TLU(X,Y,PR,NX,NO,XP,LX.EiJX.MO.MC> 
GF=CGF(XS ,XP)  SKF=LX*GF  iC<F=<F 

RET JRN 
END 


COMPLlX  FUNCTION  V<F ( XS , X P, X . Y . °W ) 

KERNAl  (voltage) 

INTcocR  =R 

complex  y 

COMPLEX  EX,;.OX 
COMPLEX  GF.VGF 

COMPLEX  kf 
C 

OIMcNilON  X(PR»,Y(PWt 

C 

COMMON/Oiv/NL),  MQ.MC 
COMMON/OEE/XI t XF.NX.MX 
C 

CALL  TLU<X,Y,PW,^IX,NO.XP.tX,ECX,MO,MC» 
GF  = VGF  ( X S » XP ) IKF=tX»GF  <SV<F=<F 
C 

RETURN 

iNJ 

......................... 

COMPLtX  FUNCTION  CGF(X,XP> 

c 

C GFt_N.;  FUNCTION  (CURRENT) 

c A -30 

INTLGtR  C^lP.pO 
RL  Al  L 


o o 


COMPLEX  AML  , AY-1,  AYP,  AZM,  AZP 

COMPLcX  tML.dYH, EYP,£2N,E2.B 

COMPLEX  OC.OET 

COMPLtX  FC.ZC 

COMPUlX  zlm.zlp 

COMPLEX  PCM.RCP.FCF 

COMPLEX  QCEX  P 

COMPLtX  GF 

LOGICAL  LGT,LLT,LTT 

COM'ION/Ai  t,  / CR,  LP  , PQ 
COMHON/D22/XI ,XF ,NX,  MX 
COMMON/322/YI,YF,NY,MY 
COMMON/M*  J/PC.ZC 
COMMON ZHIh/ZcM,  Z L r- 


.0  FOWtftT ("u'VEF-XOI-  (CGF)  » PARAMETER  Lj.MIT;>  NOT  SATISFIED", 
♦ (lFGEo.  iC) //) 


LGT=(X.GT.xF).OR.(xP.GT.YF) 

LLT=(X.LT.aI)  .J;.  (XP.LT.YI) 

L T T - L G T .O-.LLT 

IF  <L  TT  > WMTE  e_F,  id  ) XI,X,XF,YI,XP,YF  ?IF(LTT)  STOP 

L=YF-YI  {SI  = l JSF=L 

Y =SXY ( Yi , YF, SI, SF, X ) 

YP=SXY<YI.YF,SI.SF,XP> 

A D - Y - Y P ?YO  = ABS(AC)  '»  Z D =2  *L  - V D 
AS=Y  HP  tYS-AHS(AS)  iZS  = ?*L-YS 

A ML  = -L  * PC*  L S_ML  = CClXF (AML) 

AYM=-PC*YC  «EYM=DCEXD (AYM > 

A/P=-PC*YS  i- Y P = OCE X = ( A Y ° ) 

AZM=-PC*ZO  F;ZM=OC^XP(AZM) 

AZP=-PL*ZS  5^ZD=rC_XP(AZP) 

PCM=PCF(/lM, ZC) 

RCP  = RCF(Zi.P,  ZC) 

IF ( X .L  T .A  F)  OC  = c Y M*FCM»PCP*tZM-ACM*L  Yp-kCP*E  ZP 
IF (X.GT.Xr)  JC  = l YM*PCM*FCF*t  ZM-RCM*r  YP-<CF*E  ZP 

OlT=Z*(1-MCM*FCF*EMLI 

5F= ( DC/D_  T ) / ZC  ICGF  = SF 

RE  T JRN 
:no 


COMPLEX  FUNCTION  VGF(X,*P) 


jRE-.NS  FUNCTION  (dOLTAGt) 

INT^G.P  Cl-  ,LP,30 
P:  A L L 

COMP_cX  aml.aym.ayp.azm.azp 

JOMPLu  X -.Ml,.  YM.l  YP,£ZM,E7P 
CO  M 3 L _ X OVf  J-.T 
COMPl^X  JC,ZC 
COMPl.X  ZlM.ZlP 
0 0 M P u l X JCM.PCP.PCF 


I 


o o o 


c 


c 

c 


c 

c 

c 


0 0 < ur-  < P 

C0.1oL;  < '.F 

„ L U i.  y M L L ' f * k. 


T T 


JJMYJ.«/A..i./C'.  »lF  , - ) 
c o h -t o i/oze/aI.xf.nx.mx 

COM-iOJ/Tw  '/YI  , YF  , • y , my 


COM 

i £•  / P C 

, ZC 

c J N 

HU  UP 

' 

./7l 

■ U / L 

c 

Fur ' 

MAI  < " 

Lt 

* K ' . T- 

( \T 

; f ) : paf 

A M E T " P LIMIT  > NOT 

SAT  IGF 

" 

• i 

« *.  0 

( ;pc 

L C » //  ) 

lgt 

I . 

, XF  ) 

. j - . 

( , 

, C-  T . Y r ) 

ui.T 

= U .LT  , 

.XI) 

• J •*  . 

UP  , 

.lT  . Y I » 

L T T 

- Lo  T • 

x 

IF  ( 

w r r ) 

W r 

■ ir_ 

(L'i 

i J > 

* I t x ,n  « 

y i . <<■ , yf  5IF  (l  r r i 

STOP 

L - Y 

F-YI 

5. 

u-j 

$ 3 F 

r L 

EO", 


Y =SXY{YI.YF,SI,CF,X  » 
YP=iXY(Yi  . YF i j I > Sr, A XI 

Al)  = Y-YP  fY  0 = Ai5j  < A 1 ) »ZO=:*L-YO 
A 5 - Y ♦ Y J S Y S=  A • 3 ( A ; ) SZS  = 2*L-YS 


c 

Af'L  = - J*  PC*U  fr  Y L = C C E * F ( /.  Y L ) 

C 

AYM  = -F>C»YO  J..  Y 1 = DC- X°  ( A Y M ) 

AYP  = -f>C*YI  Y 3 = DCc.XP  (A  YP  ) 

AZM=-°C*ZO  $l7M=CC£xp<AZM) 

AZP=-PC*Zi  J-ZP=OCExP(AZPt 

c 

PCM-PCF ( ZLM, ZC) 

XCP=P-F ( ZLP. ZC) 

c 

IF(x.lT.X°)  OV  = i y v-~ CF'*pC  r*E  Z “♦•JO"'*  Yp--vC  = *EZp 
if  t x . ,t.xf)  Liy  = iY»'-PCM*^CP*-Z,'i-F-cM*  yjfpcp*ezp 
c 

ot  r=z*  ti-A.CM*pc=*rML) 

c 

jF  = jy/3iT  I y G F = G F 

c 

i^lTIRN 

END 

p, ».».»».. *»*,,,»»*»»»»»♦»•#•*♦**»»**********♦*»*»»»»»*******»**»******• 

PEAL  FUNCTION  SX  Y (XI  , XF  . Y i , YF  . X ) 


SC  A c I IG 

3 = (YF-YI) X(XF-XI)  JY  = SMX-Xl)+YI  ?GxY=Y 
C 

■it  T JPN 
-NO 

C*** 

COM  Dtt X FUNCTION  FCFIZL.ZD 

c 

C -XEFlElTIjN  coefficient 

c 

COM^l-.X  /i,ZC 
COM  JL:  X PC 

C A 

PC= ( Zu-ZC ) / (/L UC>  F?CF=PC  A'32 

c 


#*• 


RE  TURN 
E NO 

£*****»¥*¥¥¥**♦¥*¥¥¥¥¥¥*¥¥¥¥¥¥¥¥¥*¥¥¥¥*¥¥¥¥¥¥¥♦¥¥¥*¥¥¥¥¥¥¥*¥¥¥¥¥¥*¥¥¥*¥'• 

COMPLEX  FUNCTION  ]KF(XO*Xl,XfY*PR) 

C 

C INVERSE  KEPNAl  FUNCTION 

C 

INTEGER  °W 
REAL  KP,<I 
COM°L£X  y 
COMRLlX  KX.OX 

complex  gf,igf 

COMPLEX  < F 

c 

DIMENSION  X(R*)  ,Y (PW) 

c 

COMMON/3l.*/NO,  MD.MX 
COMMON/ 01 3/WI , WF ,NW 
C 

OATA  3I/3 . 141B9E 65/ 

C 

CALL  TlU(X,Y  ,Prt,NW,NOtXI,KX,K:OX,MO,MX»  i KR  = RE  Al  ( K X ) tK  I = A I M A S ( < X ) 
GF=IGF(X0,XI»  SGk=REAL ( GF ) tGI=AIMAG(GF> 

C 

<F  = C*(<R*GR-<I*GI)/(2*PI)  X I K F = <F 
C 

RETURN 

END 

£***♦»****»***»***»*****»**»»*»»»*»*¥***»»*•*»»***»»»*«*»*«»*»**»»*•¥**» 

COM^LcX  FUNCTION  IGF(X,xT) 

C 

C INVERSE  GP.Lc.NS  FUNCTION 

C 

INTEGER  SGN 
COMPLEX  Cl 
COMPLlX  CA 

complex  ocexp 

COMPl_X  OF 

c 

OATA  SGN/ + 1/ 

C 

C I = C M PL  X ( 1 i J i . i LI  J 

c 

CA=SGN*CI*X*XT  IGF -9CL XP (CA ) IIGFiGF 
C 

RE  TURN 
ENO 

C LIBRARY 

II  It  * t *#  ft  « it  ##**/»##  « <M  X M « M # T A ll  1 II  A # IM  M 0 ##*### 

REAL  FUNCTION  AF  ( X I , X F , NX , I X , M X I 
C 

C ARGUMlNTS 

C 

INTEGER  CR.lP,R(0 

Rl  Al  l3F.lNF,itF 
C 

COMMON/AiE/CR.LF , = 3 

C 

iu  FORMAT  ("Q’VENROP  (AF)  : PA-<AMET-P  lIMIT  j NOT  JA  T I SF  IE  0" // » 

C 

IFUX.LT.l  ) hRITE  (^R.lt  ) A-33 

IF(IX.ST.NX)  rtRIT  ( L 2 • 1 w I 
C 


oooo  o o ooo  ooo 


AF=XI  5IF(NX.i3.1)  RETURN 

c 

IF(MX.lT.l)  X = lSF  (XI ,XF ,Na, IX  t 
IF(MX.EQ.C)  X = LNF  (XI  , XF.NX,  IX  > 

IFtlX.OT.J  X=lTF(X1, XF.NX, IX) 

c 

AF  = X 
C 

RE  T'J  RN 
cNO 

R E A L FUNCTION  LaFTAI, XF.NX, lx) 

ARGUMENTS 

INTEGER  CF,,_P,P3 
REAL  liiO 

COMMON/A,.  u/CR.LP.pQ 

10  FORMAT  ("G**, "ERROR  (LSFI  : PA-AMETER  wIMITS  NOT  S A T I 3 F 1 ED  "/ / ) 

1FCX.LT.*  I WRITE(LD,U> 

IF(IX.GT.NX)  WRITE(L3,lC> 

N:lX-i 
0=NX-* 

OIF  = XF-X  I 

X=XI*  (N/0C3IF  ?LSF  = X 

RETURN 
£N  J 


REAL  FUNCTION  lNF  (X  I , X F . N X , I X ) 

ARGUMENTS 
INTEGlR  CR.lP.pQ 

real  n.d 

c 

COMMON/AlC/CR,LP,DQ 

c 

DATA  RASt/2.  7AStdl'}/ 

C 

10  FORMAT  ( " 0 " , " RROR  (LNF)  : PARAMETER  LIMITS  NOT  S A T I SF  I EO"/ / > 
C 

IFIIX.LT.1  ) WRITEUP, 10) 

IFCX.GT.NX)  WRITEUP. 1C) 

C 

LNF=0  XIF ( (IX. EQ. 1) . ANO. (XI .lU.G ) ) RETURN 
C 

X L I = 0 ilF(xI.Nc.O)  XLI  = ALOG(XI) 

XL  F = 0 IIF(XF.Nt.C)  XLF  = ALOG(XF) 

C 

N=IX-1 

0=NX-1 

OCX  =XLF-XL I 

XL  = XLlM  (N/O) * OL  X 

X=ALF(XL,3ASE)  tLNF=X 

c 

RETURN 


REAL  FUNCTION  LTF(XI. XF.NX, IX) 


.Mb*l 


noon  ooooo  o o no  noon  non 


ARGUMlNTS 

INTEGER  CR  » L ° . PO 
REAL  N,0 

COMMO  T/AiL/CR.LP, FO 

OATA  lASE/lu.C/ 

10  FORMAT ("Q**,"_RROR  (LTF)  : PARAMETER  LIMITS  NOT  SATISFIED"//) 

IFJIX.LT..  I WRITE  (lp, ID 
IFJIX.D.NX)  rt  RI  T E (L  D , 1 0 ) 

lTF  = 0 tIFJ  (IX. £0.1)  .ANC.(XI.EQ.G))  RETURN 

XL1=0  UF(XI.Nt.O)  XLI=ALUGlO (XI) 

XLF  = J JIF(XF.NE.O)  XLF=ALOGlC (XF) 

N=IX-1 

J=NX 

DL  X PXLF-XLI 

XL  = XL  I ♦ (N/O)  *OLX 

X=ALF(XL,BA3E>  ILTF  = X 

RETURN 

ENG 


RIAL  FUNCTION  AlF ( ARG , 0 A Sf ) 
ANT  I-l JG 

AL  = JASE**  a*-  J i A L F - A L 

Rl  turn 
ENO 


SUbROUTlNt  i Ob (0.°A, PO.N A, NC, ML, RL.SL I 

SCALE  <lH) 

JOUbLE  PR:  Cl S I ON  A 
C 

lNT-Gt.R  PA.PJ 
INT^GlR  C".LP.pQ 
C 

1)1  ME  NS  I UN  O(PA.PO) 

C 

COMMON /A* C/CR.lF.pQ 
C 

i.0  FORMAT  ("L'V't.R.ROF  (SOP)  I PARAMETER  LIMITS  NOT  SATISFIED"//) 
C 

IF(NA.GT.PA)  WRITE (LP.lC) 

IF(NO.GT.PO)  WRITr(La*lu) 

c 

IF(,ML.E1.0)  Rc  tu°n 

c 

00  i uu  I A =1,  NA  1 1 
JO  1 1) c NO*  1 

A»=0(i A.IOI/RL  f A - A B S ( A * ) 

J(IA,IOI=u  TIMA.  .0.3)  GO  TJ  UC 
IF(Ml.LT.j)  XIA.IO)  =SL*OLOG  (A)  A'35 

IF  ( Mu.  GT  . t,  ) O(IA.IO)  =SL*OLOGlu(A) 


o o 


n j co  *j  r i --rj_ 
c 

Rc  T Jv.M 

1NO 

f»M*»»M«»*M»r*»»M»M»*MMMf»»t»**f»»>»»*»*«*»*****«*«M*«*******«* 

iUJ^o jr  i inmx  (o  t pa  . pj  * na * no,  xmn.xmx  .iin..^  x . 10 » 

c 

c -ii.-i/  i»x 
c 

INTEGER  PA,P0 
I N T _ ->:_  - C'.l3.^ 

REAL  TIN, MAX 

c 

DIMENSION  0 ( P A , c 0 I 

c 

G 0 M M 0 *) / A 1 X / C x . L F » " ) 

c 

*0  FORM*!  <"l "."ERROR  (MNMX)  ! oApAM£  T£*  LIMITS  NOT  SA  T I 5F  IF  C"  //  ) 

c 

IF(NA.GT.FA)  rt-elTL  (LMl  1 
IFCNO.GT.PQ)  WRITE <lP.1G ) 

C 

4MN=0(l.i > 

AMXr0(l*») 

00  * 0 J II  -ii  NA  i * 

00  1 U j IO-I1NO14 

AMN  = A.MIN1<A1N»3(IA,I0)  » 

A MX  =AMAXi < A -lx  , J ( I A , 10)  ) 

1 1»  u CONTINUE 
C 

1 F ( M 0 ) lul^jcubl 

1C1  CONTINJl 

XM.-js  AMO 
XMX  =AMX 
GO  TO  1 1 u 
1J2  CONTINUE 
XMN=MIN 
XMX  =MmX 
GO  TO  liil 
105  CON TINU_ 

X MN=MIN  JIF (XMN.GT. AMN)  XMN=AMN 
XMX=MAX  JIF (X  1X.LT.AMX)  XMX  = AMX 
110  CO  NT  I N'  J _ 

IF (XMN.GT . 0 ) XMN-0 
IF(XMX.LT.u)  XMX=0 
C 

RET JRN 
2ND 



SUORJ'JT  iNt  TOi_(A*C*pAtPO.NA.NO) 

C 

TABuL 

INT EGER  c A , P0 
INT^GcX  CPtLP.PO 

JIM-NGION  A(pAI ,0(PA,P0* 

COMMON /Alu/CR»LP«°Q 
C 

10  FORMAT  ("u"."c.kROR  (Ti3L>  « PARAMETER  LIMITS  NOT  SATISFIED"//) 

11  FORMAT!"  "«IlC*iO(l°G2o.lC)> 

C A -36 

IF  ( NA.GT.  F-A)  WRITE(L°tlC) 


c 


IP  ( NO.GT  . POT  ^IT-:(L3,;f  I 


00  a o u IA  = 1,*(A,1 

WRI Ti ( L o . : 1)  IA»A(IAT  . (P(IA.I0»  , IG=i ,N0, It 
lCo  CO  NT  IN  Jl. 

C 

RETURN 

lN[) 

* 

SUoROUT I ALT<».  .PA.PO.NA.f  0. MIN, MAXI 
C 

C 3lO  T 

C 

C ° A9  A Me  T X Ni  = 'i 

C PAR  AM_T  -1  - PP  = . J I 

C 

INTtGf<  PA.J0 
INTe&EP  °~ 

INT_Gl»  CP.^P.°C 
INle&eP  * J 

1 N T G _ R axIS.OLNK.SMCL 
i<l!l  MIN.  MAX 

c 

OIMENjIJN  A < pa ) , 0 (PA ,F0) 

01  M_  NS  I ON 

UIM_NJICN  lIN:.  (1111 
C 

COMMON/Ai  j /L-: . Lc.  rC 

c 


OATA  I MN , IMX/i, lCl/ 

OATA  2amo/:x 
J A T A AXlaV"*"/ 

OATA  0LN</"  "/ 

OATA  5 M3i  /"l  ” , ” Z “ • " 5 " , M a"  . "■}”  , ,,feM  , " 7"  , •*  * ",  "9"/ 
OATA  NS,3-  / (,.ul/ 


C 

11  P OP  M A T ( j .IoX.PX.I^jIu.  T , ' t u X , 1 p G 1 G . A,2x) 

.2  FORMAT  lux, 10X..01AU 

13  POP.  TAT  ("♦",I5»5X, 1-09,2.1X1 
l.  c OR  MA  T C*  ••) 

C 

^L0P_=G  * IF  ( MIN.  Nt.  .MAX  ) S L 0 3t  = ( I M X- I T \ ) / ( M AX  - M I N ) 
C 

WRI Te (L- , il»  MIN, MAX 
C 

00  1 j.  IP  = l,3t*,l 
LINE ( I3>  =AXIG 
lu.  CO  NT INU_ 

WRiTe(LP.l-) 


WRlT_(LP,»t)  (LIN  (I°I,IP=1,PP,.) 

C 

4 0 0 1a*  Im— 1,N»*,1 

WMi  T_ (LJt let  I a , A ( I A ) 

00  lie  I0=.,Nu»1 
00  ili  IP=1,P-,1 
L I N . ( I 3 1 = o L W < 

113  CONT INJ_ 

I2=IlOPl*  (7.PO  -MIM41  li. I N_ ( I n =AX I 3 


o n 


! 


i 


i 

£ NO 

C****»*****»*****»*»****»*»****»»***»*****»**********************»****** 

SUORUOTINc.  o \F 

♦ (A,0,PA,PG,IA,\,C,YIN.FAX,TT,TX,T\/,TMtPM,TS,FN,P,KF,''K,NM,MCtMC»MG> 


c 

c 

c 

c 

c 


c 


c 

c 


CAcCO  iP 

PAR  A M-_  T£ 

",  PG  = 5*2 

P A R A M _ T £ 

P 1 = X J £ 

PAR AM_ T £ 

« P j -5 

INTcGi* 

P A , p J 

INTGG. r 

Pl-I.p’l 

INT-^G.R 

C--  ,lP,pQ 

I NT  G Gl  R 

P r , J J , p 3 

INTtGe  «! 

T T , r X , T V 

INTcGlR 

T o . 5 I 

INTcGcR 

OF" 

INTu&lR 

PF 

INTlGl  R 

ON.  GO 

INTtGGR 

PH  » 3 V 

INTEGER 

3 P C i S PS 

I TA  T G c R 

j • C S • x c 

I N T l G _ R 

101,  c. OI, 

Ri  A L 

MIN, MAX 

RlAl 

MOD 

LOGICAL 

L I (J , L Gvj 

JI M£N » X ON 
OIMJNji JN 
DIMENSION 
0 I H_  NS  I ON 
JIM_  IS10N 
DIMENSION 
DIMENSION 


- ( = , 0 (PA . PO> 

TT (PM) , TX ( PM) ,IV(PM) , TH I’M ) 
T , ( PN> 

13FP  (r,U  ) 

On  T A ( j.  u £ ) 

(5) 

Xrl  <5  > , Y“ (5  ) 


COMMOI/Aiu/CR,LF,cQ 

DATA  PS/5  !£/ 

DATA  PO/1 L £/ 

DATA  P5/5/ 

DATA  X0,Y0/Yx.u,+3.5/ 

J A T A XI/,  Y //*(j  . 0 ,*0 . 0/ 

DATA  XH , YH/t J . U . ♦ 3.G  / 

J A T A X 
J A T A 

J A T m XT  ,YT/  + *..JG,*2.CG/ 

DATA  X^,  YP/«-x  . 00  ,-L.ib/ 

DATA  x '),YN/*2,00,-2,25/ 

OATA  XC, YC/ flu . 0 , +0 . 0/ 

DATA  XlfYl/YU  u.+O.C/ 

DATA  XF,YF/Fj.G  , + Cl.O/ 

DATA  po/Ci,/ 

DATA  0N*00/-3, *3/ 

DATA  PL/12. 0/ 

OATA  H l • V L / *6 • 0 , ♦ 5 • 0 / 

JATA  HA,VA/  + i).0,+T3.0/ 

DATA  HS,Ai/»0.1U5,«-0.C/ 

DATA  *F, AF/*u . 1G5. H . C / 

OATA  HN,AN/*u.iG5,»0.0/ 

JATA  JP, A F/+u . 105 , ♦ J . 0/  . 3ft 

JATA  Ol/f’.j/ 

DATA  IW/H/ 


OATA  LC/00/ 

OATA  NT/8/ 

OATA  SS/*i/ 

OATA  CS,SC/-1, *10/ 

OATA  SO/+1.J/ 

OATA  NO/+2/ 

J A T A NS/ *5/ 

OATA  S 3 S / c . 1 ic  » 3 i 4/ 

OATA  SOI ,£ OX , EOF/1 , -999. *999/ 

C 

IF(MIN.tQ.MAX)  R£  TURN 
C 

LIN=. FALSE.  ii F ( (MA . EQ. C » . AND . (MO.cQ. 0)  ) LlN=.TRU£. 
i_OG  = . FALSE.  BIF(  (MA.NE.Ci)  . OP  . ( MO.  Nc.  . 0 ) ) LOG  = .T-?Ul. 

C 

I F ( (MA.NS.G)  .ANO.  (MO.tQ.C))  LT  = -l 
IF((MA.N£.u) . ANO. (MO.NE.G ) I LT=  0 
IF<  ( MA.EQ.G)  .ANO. (MO.NE .0 > ) L T = ♦ 1 

C 

ne=MG  *IF (MG.tQ.EOI)  M0=8OI 

IF ( M8.EQ.bOI ) CALL  PLOTS(bFR,P0, PQ, «) 

IF(MO.tQ.BOI)  CALL  PL OT < X 0 , V 0 , ON ) 

C 

CALL  PLOTMX(PL) 

IF (MA. EQ. U)  CALL  SCALE(A,HL,NA,SS) 

IF(MA.NE.O)  CAwL  LGSCAL ( A .HL.NA , SS) 

NC=NM*LW  ilF(NC.GT.LC)  NC=LC 
PH  = -NC  tP'/  = +NC 

AMN=ABS(MIN)  *AMX=A3S (MAX ) t MOO  = A M A X 1 ( A MN , A M X ) 

VT  M = -MOD  irtTM  = A (NA*SS  + 1) 

voy=+moo/ui  jhqx  = a (na*ss  + 2) 

IF  (MO.  El.  o ) CALL  AXIS(XV»YV.TVtcV«VL.VAftfTM.VOY) 

IF(MO.NE.l)  CALL  LGAXIS  (XVt  YV,  TV,  FV,  VL.  V A,  VT  M,  VDY  ) 

IF  (MA.tO.b)  CALL  A X I S (XH, YH , TH , PH , HL. M A , HT M, MOX ) 

IF (MA.NE.CI  CALL  LG A X I S ( XH , YH , TH , CH . ML . M A . H T M , MDX ) 

CALL  PLOT ( XI , Y I , 00) 

30  1 G 1 I 0=1 » NO  t 1 
00  1J2  IA=1,NA,1 
OATA(IA)=C(IA,IO) 

102  CONTINUE 

OA  T A ( NA  ♦!  ) = VTM 
DATA  ( N A ♦ 2 ) = V JY 
oFC=SPS(IO) 

IF (LIN)  CALL  L INt  ( A »DATA«NA.SS*NS. S°C) 

IF (LOG)  CALu  LGLIN£(AtOATAfNA,iS»NS«S°C,LTJ 
XS=XM(IO)  iXO=XS*SO 

ys=ym(io)  jyo=ys 
ST=TS< 10) 

NK=MN*LW  JIF (NK.GT.lC)  N<=lC 
CALL  SYMBOL ( XSfYS.HS, ScCtAS»CS) 

CALL  SYMBOL < XO.YO,  MS, FT  »AS»SCI 
lGi  continuc 

CALL  SYMBOL  ( XT  ,YT  .MF  • TT  ,AF  « NK. ) 

IF  (MO.. Net)  CALL  NUMBER  ( XN,  YN,hN  , c , A N * NO  ( 

IF(MP.NE.C)  ca  uL  symbol ( xf, YP.HF, Tx , ap, no 
CALL  dlOT  (XC.YCtOM 
C 

Mt=M3  BIF(MG._Q.tOI)  Fc=£CF 

IF ( M£. £0. EOF)  OALw  °uCT (XF.YF«£OF) 

C 

Rt  TURN 
t NO 

£*#**»**#.*****.*»****#¥*****♦***♦♦#**  a - 39  **************************** 

s')8RouriNr  r-’Gx 


I 


o o 


< A , ) . r 


- . t r , r a , r v. th , pm 


J-4  T - 


P i ' , a a « A F , N A , M X , X a I • I X F , N X I) 


c 

c 


c 


c 


c 


c 


c 


tails/ ?uT/b  -a  ta 


- A < A M l.  T ■ 

R P J = l 

I M T _ ol  R 

F r 

INT_G-  R 

P M , P N 

. r,  | _G,_R 

D A ,->0 

I NT  i 

CF ,LP. PJ 

InT  j.R 

TT.TX.T7 

INTEGER 

T G 

REAL 

MIN,  'll  X 

IJM JL. X 

DATA 

dimension  data(po) 

JIM  .NSICN  A ( ■»«  ) , 0 ( 3A  , °0) 

JI M_  No  ION  Tr<PM>,TX(FM),rv<PM),TM(PM> 
0 I M N o I 0 N T X f ■ ) 


) ,M  M 0 • X A * l / J R * L P * c Q 
G0M1J  T/H12/.ML.RL,  S<_ 
COM  M ON/ R MIN, MAX, MC 

COMMON/  ) 2 U / 1 X,  NP 


1 FORMAT ( M A 1 u ) 

2 FORMAT ( " 0 “ , 1Aij) 

J FORMAT (iu ilJ ) 

- FORMAT  ( 2oiu. J , lit  t 
*1  FORMAT ("0 


INDEX" 

A F tj  U M t.  N T 


2 


♦ 


REAL" 
IMAGINARY" 
MAGNITUDE" 
ANGLE  (RADIANS*" 
ANGLE.  (O-GREES*"//) 
->  =£AL  "/ 

->  IMAGINARY  "/ 

->  MAGNI TU3L  "/ 

->  ANGLE  (RADIANS)”/ 

->  ANGLE  (DEGREES)"//) 


J A T A RN/o/ 

DATA  To 

♦ /’’R  c.  Ai_  "."IMAGINARY  "."MAGNITUDE  " . "ANGLE. 

)A  T A R.MJ/L. 0 , U / 


A 0 ( C-x » 

i ) 

(TT  (IT)  , IT  = 1,PM,1) 

RlA0(CR, 

1 1 

( T X ( I T 1 , I T = 1 , PM, 1 ) 

• . < 

1 ) 

( T /(  IT)  , IT  = 1,  PM.l ) 

1 A D(CR» 

1 ) 

( TH ( IT) , IT=1 , PM, 1 ) 

R A0(CR, 

2 ) 

MA , MO, MT , MP,  MG 

A )(CR« 

4 ) 

MIN, MAX, MD 

"."ANGLE 


IFKMT.EQ.OI  .AND.  (MP.tG.U*  «ANO.  < MG  . c.  Q . 0 ) * Pt  TU3N 

IF ( ( MT . N_ . u ) . OR. ( IP. NE. 0 t ) WPITE(LP,  2)  ( T T ( I T) . I T = 1 , PM, i ) 


C 


-ALL  A0Z(A*0,PA,P0, DATA, PC', ND»YI*XF,NX*MX»IXI*IXF,N<U 


■ < A = N J (Nu  = 5 
I f ( m J.NE.O ) 
N A s N J (‘10  = 5 
i - (MT.Nt.fc* 

i F ( M r . J : ,U) 


CALL  sod (O.FA .po.na ,n 

R p I T t (IP, 11) 

call  tml  (A.C.PA.PO.NA 


. ML, PL, SO 


NO) 


A -40 


o o 


.NA  = ND  JNJ  = 3 

CALL  IN  MX  { 3,  °A  , PC,  NA  , NO  , XhN  . XMX  , MlN,  MAX,  10  » 

(A  = ND  SNJ=3 

IF  ( 1°. Nl . C)  PFITE  (lP.12> 

IF(MP.NE.o)  CA..L  CLT ( A ,0, PA , D0, NA.NO, XMN. XMX ) 

N A = N D i N 0 = 3 
IF(MS.NE.C)  CALL  G«f 

♦ (A,  0 , 3A  , PC,  NA,NC,XMN,XMY,TT.TX,Ttf,TH,3M,TS,PN,°,MC.MN,NM,MA,MC,  MG) 

pe  r jpn 

£ NO 


SUdPJUTINC  FPGK 

+ (A,o,°A,pj,rr,rx.T\f,Th,PM,oATA.Dx,or, 

♦ XI,XF,NX,MX,IXI,IXF,NXI.YI,YF,NY,MY,IYi,IYF,NYI) 


GF  A 3H  ( T-PL£/DLOT)  (FFIQUt.NCY  DOMAIN) 

c 

I N F £.  p 3 X , 3 Y 
INT.oIX  3A.nJ 
XNTlGiP  3 M 
XNr^.,_x  TT,rx,rv,rH 
COM°llX  DATA 

c 

OIM.NilON  A (PA > , G ( PA , P0> 

DIMENSION  TF  (PM)  , TX  (PM)  , T \/  ( o N ) , T H (PM) 

DIMENSION  DA  T A ( °X  , DY ) 

C 

DATA  NX/2/ 

C 

DO  ajl  X ' : * i 'O  * * 

XF(XP.^T.l)  CALl  \z=>xY 

♦ (VI,\/F,NV.MV,XVI,IVF,NVI<C5I«PF,NP,M3fX=I,iPF,NaX, 
i XI,XF,NX,MX,IXI,IXF,NX1,YI,YF,NY,MY,IYI,IYF,NYI) 

XF ( X 3. cX. I ) CALL  V°x  y 

♦(VX,OF,NV,MV,XVX,:VF,NVI,fI,CF,NP,MP,IPI,XPF,N°I, 

♦ YI, YF, N Y, K Y, IYI, IYF, MY  I, Xl,XF, NX, MX, IX  I, IXF, NX  I) 

CALL  T3GC 

♦ <A,J,Pd,o0,TI,TX,T\/,TM,PM,JATA,PX,PY,IP, 

♦ VX,WF,NV,Mi/,iVX,IVF,Ny/Itf'I,FF,.NP,'iP,XPI,XPF  , N°I  ) 

1UC  CONTlNJc 

C 

PlTUPN 

-NO 

p****»,***»», *»****»», ************************************************** 

iuJ-XJlJTIN;  TPGC 

♦ (A,O.PA,30,Tr,rx,TV,TH,rM,DATA.pX,c>Y,MZ, 

x XI,XF,NA,MX,XXX,XxF,NXX,Yl,YF,NY,MV,XYI,IYF,NYI) 

C 

C T Adc-/aL DT/GPA °h  FATA  VS  > OF  Y 

C 

C pA*M£T-:f  pN  = -> 

c 

XNTcGl?  3x,3Y 

xn r.cr?  pm 

Xf.T--  IlX  3 , 3 J 
I N T . ,,l  3 C - , L P , 3 0 
INT-o_P  XT  ,Tx,TY,TH 
xntlG.p  :s 

I n T w G : 3 )Ji,CCl,t  01,  ”CF 
PE.  A L N IN,  MAX 

C0M3ucX  JATm 

A -41 

DIM  N _>  I C N nArA(P».PV| 


r\ 


c 


c 


c 


c 


c 

c 


JIMuiNilON  A ( PA  ) , 0 ( 3A  , PO  > 

UlM:  Na  ION  TT  ( -Ml  , TX  < PR)  , T V(pM)  , TH (PM) 

dimension  r.is) 

COMMON/ A it  /CR.L<=  . "Q 
COMMO^/  U'/  T l , XL  iSl 
COilMON/Blc/TIN.MAX.MO 
CCMMJN/XeL/MN,  NM 


DATA  PN/B/ 

DATA  IS 

♦/"RcAl  " , I M AG  I N AF  Y "."MAGNITUDE  ".“ANGLl 

DATA  H/./ 

DATA  BO  I , COI.lOI . EOF/ 1,2 . -999, *999/ 


1 FOk  TAT ( 6 A 1 0 ) 

FORMAT  ('•u",  BAiu  ) 
FORMAT  (*GIiJ) 


4 FORMAT (2G10. 0 
1 1 FORM A T ( "u " t " 


•f 


♦ 


0 ) 

I MOP  x" 

ARGUMENT” 
K _ A L " 
IMAGINARY" 
MAGNITUDE" 
angle  (RADIANS)" 


1 2 F GRM  A T ( 

+ 

♦ 

13  FORMAT  ( 


ANGLE  < DEGREE  5 >"// ) 

1 - > c E A L ”/ 

2 ->  IMAGINARY  "/ 

i ->  MAGNITUDE  M/ 

h ->  ANGLE  f P A D I A NS ) " / 

5 ->  AnGL  r (DFGRFCS)"//) 
PARAMETER  :",iPG2C.lO> 


READ(CR,  1 
RcAO(LR,  j. 
READ (CR,  1 
READ(CR,  1 
REAO(CR,  2 
RE  A D ( C R » A 


(TT(IT) , IT=1,PM,1> 
(Tx  (IT)  , IT  = 1,F,M,1) 
( TV ( IT) , IT  = 1, PM,1 ) 
(TH( IT) , IT=i,PM,i) 
MA , MO, MT , MP, MG 
•(IN,  MAX  . MO 


IF ( ( MT .£0 . 0)  . AND.  (MF.lO.O ) . AND.  (MG. i 0.0 ) ) RE  TURN 
IF  ( ((MT.  Nt.  0 > .OR.  (MP.  Nfc  . ( ) ) WRITE(L°.  2)  (TT(IT),IT  = 1 


DO  10U  IP=IYI,IYF,NYI 
MC  = C 01 

I F ( ( ML-  . £Q  . 8J  I ) . A NC  . (IF  .EC.IYI))  MC  = DJI 
IFUMG.dQ.tOFI.ANO.  (IP  .dO.IYF))  MC=tOF 
IFUMG.tO.tOI  > . AND.  (IF  .EO.IVI))  MC=0JI 
IF  ( ( MG  • E 0 • 1 0 1 > .AMO.  (IP  .ED.  I YF)  ) MC  = EOF 
IF(  ( MG .£Q .dO I ) .ANC.  (IYI.EO.IYF)  ) MC=EOI 
CALL  aoxy 

♦(A,0,PA,P0,0ATA,PX,PY,NZ,P.IP.M2. 

♦ XI,Xf,NX,MX,IXI,IXF,NXl,YI,YF,NY,MY,IYI,IYF,NYI) 
NA  = NZ  «N0=5 

IF(MO.NE.O)  CALL  SDO (0, PA , PO, NA , NC, ML.RL . SL) 

IF ( (MT.Nt .0) .OR. (MP. Nt .C» ) WRI TE ( L P , 1 A > ° 

N A = N Z SNO  = E 

IF(MT.ME.G)  WRITE (L°, 111 
IF(MT.NC.G)  CALL  TBL(A,O.PA,co,NA,NO) 

NA=NZ  iNO=3 

CALL  MNMX (0, PA , PO, NA .NO t XMN, Xhx ,MIN, MAX, MO) 

NA=NZ  SNO=3 

IF ( M°. NE . L ) WRITE  (L°, Id) 

IF(M°,NE.wI  CALL  plT  ( A , 0,  PA  , pO,  NA  ,f)0  . XMN,  XMX  ) 


ANGLE  "/ 


PM.l) 


on  o o oooooo  non 


Nfl=NZ  f N0=3 

CALL  GRF 

MA,OtPA,PO»NA,NC.XMN»XMX,TT,TX,TV.TH,°M,TS.PN*P,MQ,MN,NM,MAtMOtMC) 
lOli  CONTINUE 

c 

r c r u r n 

END 

SUBROUTINE  AOZ(AtO,PA,°GtOATA,PU,NOtXIfXF,NX,MX,lAl,IXF,NXI» 

AESICCA/ORDINATE 

INTEGER  PC 
INTEGER  P A » a0 
INTEGER  CR,LP,°0 
COMPLEX  OATA 

DIMENSION  UATA(PD) 

0 1 Ml.  ns  ION  A ( PA  ) , 0 ( Pfl  , PO) 

COMMON/ Alu/C*,LF,FQ 

10  FGKHAT  C’L", "ERROR  (AOZ)  » PAPAMcTER  LIMITS  NOT  SATISFIED"//) 

IF  ( (Nx.GT  .PO)  .OP.  (NX  .GT.PA  ) ) WRIT-<LP,ilj> 

10=  IXF-IX  I 
ND=I D/NXI ♦! 

IA-* 

DC  iCl  IX  = IXI,IXF,NXI 
A(IA)=AF(XI,XF,NX.IX,MX) 

0( IA ,*) = REAL (DATA  CIXI ) 

J ( I A tt  ) = A I M AG ( DA  T A (I  X ) ) 

J ( I A , o ) = CA3S(0ATA(IXI) 

An>  = cAnG(DATA  ( IX)  ( $AO  = ROF(AJ) 

J ( I A , H ) =AA 
0(IA,5»  =AC 
I A = I A *1 
ACC  CONTlNUt. 


RETURN 

END 

SUSP  OU  TINE  AJXY 

♦(A,0,PA,PO,DATA,FX,PT,NZ,P,IP.MZ, 

♦ VI,VF,'JV.MV,IVI,IVF,NVI,PI,PF,NP,MDfIPI,IFF,NPI) 

AESICCA/OPDINATE 

C 

INTEGER  PX,JT 
INTEGER  PA.^O 
INTcGtR  CPfLp,3^ 

COMPuc.  X DATA 
LOGICAL  LVX,L°Y 
LOGICAL  lVY,ljX 
LOGICAL  Ll/ A 
C 

DIMENSION  UATA(PX,PY) 

0 1 Mi.  N j I ON  A ( =A  ) f 0 ( PA  , PO) 

C 

COMMON /Aii//CR*LF»FQ 

C A -43 

AJ  Fj^UT(",",".<<C-  ( T PGC  ) s «A=  AM-f-^R  _ I M I T S NOT  SATISFIED"//) 


oooooo  o oooo  o ooo 


c 


i_  V/  X = ( M .£0..  ) ,H0.  (NV.GT . PX  > 
uPY  = MZ.:_Q.i>  .AND.  <N=.r,T,  ry  ) 
i_VY=(MZ.;_0._l  . ANU. <N°.OT. PX) 
lFXs  ) ,A\ 10.  (MV.r.T.Pf) 

L / A = (Nrf.GT  • 3 A ) ,OP.(NP.GI  .PA) 

C 

IPtU/X.  J-?.L  »n  W-ITE(LP,Il) 

IFCLy/Y’.C-'.LPX)  IU  (Lc»  li  ) 

1 F ( u V *, ) /ir-  IT*  (lS5!  iu  ) 

C 

ID-Iy/F-IVI 
NZ=I  J/NVi  ♦ * 

C 

° = 6 F ( P I , °F , NP  * i=  , FO) 

C 

IA  = i 

JO  iju  I/=I/I  , I/F.NVI 
A(IA)=AF<VI,VF,NV,XV,M/> 

CALL  XYV  = <IX,IY,I\/,I  = ,MZ> 

0(14,*)=  (■  £-,_  < OATA  (I  X , I Y > ) 

0 ( I A , t ) =AIMA&(DATA(IX,IY)  ) 

0(IA,i)=  CA3j(0ATA (IX.IY)  ) 

AR=CANG ( DATA ( IX , I Y ) ) JA0  = FJF(AF) 

0<XA,L)  = A *• 

0 ( X A , -j  ) = A 0 

1 A = I A 1 
lOu  CONTINUE 

c 

RET  URN 
END 



SUdROUTlNc.  /PxY 

♦ <V1,V/F,NV.M\/,IVI,I\/F,NVI,FI,FF,Np,MP,IPj.,IPF,N3I, 

♦ Ai,AF,NA,MA,XAI,IAF,NAItOI,OF,NO,NO,IOI,XOF,NOI) 


7ARIA8LE/c«-<6Mc  TlR 


i/I  = 6I  51/1=1 
V F = AF  5H/F  = Nd 
NV=  N A 5N/I  = * 

■N  V = M A 


$PI=OI  JIPI=IOl 
|Pf=0F  fIFF=I0F 
JNP=N0  iNPI=N0I 
i M P = MO 


RcTJRN 

END 


SUBROUTINE.  XYV3(IX, IY, IV, I3,MP> 


1N0*X 


IF(MP.EJ.l)  IX=XV  KIF (MP.tl.Z)  IX=IP 
IFIMP.tT.i)  i Y = I P $IF(MP.tQ,2)  IY  = IV 

RETURN 

END 


SUBROUTINE  T L U ( T X , TY ,PT »NT»NO»X»Y,YO»NO,MX) 

TABLE  LOOK  UP  < I N TERPOL A T ION/* X T R A P J LA T I ON/O  I FF£ RE  NT  I A T I CN ) 
PARAMETER  OP=lO 

A -44 

INTEGtR  dp 
INTcGER  3 T 


J 


1 


INTlGlR  CR.LP.PQ 
COMPLEX  T Y 
COMDt_t  X QKMJ 
COMPLEX  U,QO 
COM°Ll  X Y,YD 

c 

DIMENSION  TX  (PT)  ,TY(PT> 

DIMENSION  P(HI  ,0(10  .QD(10) 

C 

COMMON/AIu/CR.LP.pQ 

c 

OATA  L>P/iu/ 

C 

10  FORMAT  ("u'VlPROP  (TLU)  I ° AR  AM£  TF  R LIMITS  NOT  SATISFIED"//! 
C 

IF(NT.GT.PT)  WRITE  (L  c , 1 C ) 

IF(NO.GT.NT)  WRITE ( L D , 1C ) 

IF(NO.GT.lG)  WRITEClP.IU 
IFCNO.LE.  1)  WRITE<LP,1C> 

c 

YD  = Q 
NTL  =1 
NTtJ  = NT 
NDH=ND/2 

100  CONTINUE 

NTO=NTU-NTL  $IF(nTD.LE.1)  GO  TO  110 

NT  A = ( NTl  + NT'J)  /2 

IFTMX.LT. G)  ASSIGN  lOl  TO  IX 

IF(MX.GT.G)  ASSIGN  1G2  TO  IX 

GO  TO  IX,  (101,102* 

101  CONTINUE 
IF(X-TX(NTAM  115.1*1.112 

1G2  CONTINUE 

IF (X-TX (NTA) ) 112,111,113 

113  CONTINUE 

NTL=NTA  {GO  TO  100 
112  CONTINUE 

ntu=nta  j Go  to  iCO 
111  CONTINUE 

Y = T Y (NTA)  £IF(MO,'Q.3)  RETURN 
NTL=NTA 

110  CONTINUE. 

NOL=NTL-NOH+l  BIF (NDU.GT. NT)  NOL=NT-NOii 
NDU=NTL*NOH*l  SIF(NUL.LT.  1)  NDL=1 
DO  131  J- * » N D SK  = NOL*J-l 
P ( J)=TX(K) 
a C J ) = T Y ( K ) 

ao<  j>  =c 

131  CONTINUE 
I=NO-l 

DO  lwl  J=i,I  <L=J*1 
DO  lw*  < = l , NO 
°<MJ-C(<) -P( J) 

1KMJ=U(K) -Q ( J ) 

XMPJ=X  -P(J) 

XMP<=X  -P(Kl 

•1  (<)  - ( < M F J * ] (K)  - X M3<  * 1(J)  * / P<WJ 

lQ(<)R{D<f'J*XM=>J*OD(<>-XMP<»aO(J>)/°<,1J 
1L1  CONTINUE 

Y =1  <ND)  tic(X.tO.TX(NTA)l  Y=TY (NTS ) 

YO=DU(ND> 

C A-45 

RETURN 

END 


I 


ooooo  oooooo  ono 


C**y************<M,***#**j,»*#»********i,*¥, ♦****♦**♦***♦♦#♦•#**********♦*# 

COM^LcX  FUNCTION  IF ( XI , XF ,'U  , MX . np, F )CT  4, X ,Y , PWI 

NlWTON/COTlS  DARAPOLIC  RULl 

INTEGER  P W 
INTlGlR  Ck.lp,'1'! 
ccm-’uc. < y 
COMPLEX  L VIUJJD 
complex  yi.yf 
com  jllx  y , ro 

COM jLl  X FNCTN 
COMPLEX  I 

logical  r;  st 

JIM-.NjIJN  X ( * , ¥ ( ^ ) 

COMMON/AIl/CR.LF, = 0 
COMMON/DEi/X  X, XXX 

10  FORMAT  (■•.•’,  "-xXOP  (IF)  : fA^AMCT'F  LIMITS  NOT  SATISFIED'*//) 
call  PARirr(  vp,  u.  st; 

IF(N°.LT.5>  ^ITl  |lP,U)  flF(TEST)  RkITw (tn. li > 

1=0 
ML  = -1 
NX  Mi = NX  - i 
NPM1 =N3- l 
NPM2  = M°-2 

30  1 0 u IX=1,NXM1,1 

IM=IX 

I°  = IX+i 

XM=dF(Xl,XF,NX,IM,MX) 

XP=AF(Xl» X F, Nx , IP, MX ) 

X 0= (XP-X  M >/NPMl 

YI  = FNCTMXX,XM,X,yfPW) 

YF  = FNCTMAX,XP,X,f,PW> 
cVt.  N = U 
000  =0 

00  lOi  IE  = 2 « NPMi , 2 
XE=AF<XM,XP,NP.IL.ML) 

YE=FNLTN(Xx,X:l,X,Y,PR) 

£VtN=.;LN»r-; 

101  CONTINUl 

JO  102  I0=i«NPM2,2 
XO=AF(XM,XP,NP,IO,ML) 

Y 0 = F NC  TNfXX«XO,X,Y*°W) 

ooo=ooo»yo 

102  CONTI  N'J- 

1 = IMXO/3)*(YI*4*E\»lN*?*OD0*YF> 

lOi)  CONTINUt 
I F = I 

RETURN 

END 


SUOROOTINc  PARITVTIT.TFST) 

OOU/lYEN 

LOGICAL  TEST 


A«46 


IC=  < 1 T/2  > *2 

IFlIC.t'J.IT)  T_pT=.TRUc. 

IFCIC.Nc.ITI  T£ST  = . FALSE. 

Rc.  Tl)n  N 
ENO 

COMPLEX  FUNCTION  QF ( X J , XF , NX , MX , NO , F NCT N . X ,Y , PR > 
GAUii/Lc!o-.NJ^£  JUAOrATUfE 
PARAMETER  6u=5 

integer  r,L. 

I N T E G L 3W 
INTEGER  CR*LP,°Q 
COMPLEX  it 
COMJc_X  I 

COMPLEX  T IkH) SUM 
compl.x  fnctn 

OIM^NjION  X ( PR)  . Y ( PW ) 
i9IM  ENSIGN  RL(5*5) 

DIMENSION  W._  (5. 5) 

COMMON/AiO/CR.LFt = 9 
C0MM0N/021/XX.XXX 


1 SL/5/ 
i Rl. 

lUUJuuCoG 
iCJjOjCuu 
lOCjjollu 
. JoGoOCGQ 
lUUOGuoo'J 
1 RL 

> UUuOoOLJ 
lOJUUJUU 

.QGGGwGCu 

•UUUuJuOJ 


, - G • 5 7 7 3 5 G 2 7 « 
, *G.67735G£7, 
,*u.CCCouoGC, 
.♦J.uOGuOOCO, 
.♦u.ClGCQGjC, 

.♦i.OGooCGGC, 
* * I * uCCCOuuo  » 
.♦O.GOOCJGOO, 
.♦u.GOOGGGCG, 

Uj.t'jjiiui/3!  v 


-G  • 7 7m  5 Job  7 ♦ -G  • 8oll3o3i  » - 0 < 9G617984, 

♦ u,CJ000uo0»-Q.3J9991Ji»»~Q.53fi46931f 
«U . 77459667 i»0 .33999l3o»+6. OCGOOGOO, 
♦0,COOOOGbO,*0. 8611 3631, ♦0.53846931, 
♦0.0GD06CG0«+u«GoCG0630,*0,5C61798m/ 

♦C,55p55556»»o.3m785m3m,+Q«236J2689» 

♦ 0.  8 8889889,  ♦,:  .6521 -,51m,  ♦3.47862869, 

♦ 0. 55555556, ♦0.6521*»51h, ♦0.56888889, 

♦ u.C0CG0ti00,^G.3‘«785H94,»d.‘»786u669, 

♦ C . Co  jOOOoO , ,0 • G00G  u u 3 3 , +0. 23692689/ 


FORMAT  ( "o  “ * "r.  RROc  (OF)  t PARAMETER  LIMITS  HOT  S A T I SF  IE  O"  / / » 
I F ( r 4 ,J  . G T . t,  o ) RR I T . ( L c ♦ 1 G > 


x = c 

NXM»=NX-1 

JO  1 0 1 1X=1,NXM1,1 

1M=IX 
I F — X X ♦ 1 

XM  = AF(XI, X F , N X , IM,MX» 

XF=AF(Xi, XF, NX , IF, MX  ) 

X A= (XP-XM ) /2 
XC  = ( XP+XM I /£ 

SUM  = 0 

90  * J 2 19-1, NC*^ 

X S = X C tRo.  ( HG.  iO)  *x  A 

TEm  M = RL(N  U,n»*F  NOTH  (XX. XS,X,T,PW  I 

SUM  - SUM  ♦ T i;  RM 

CONTINUE 

I = I ♦ X « • SU  M 

CONT  X N 'Jo. 

)F=i 


o o o 


r j R N 

EN  j 


REAL  FUNCTION  NR  F ( X , F F , FO  F ) 

NCwT 0 J/RA FHSON 

1 NT  E Gl  R C-<,La,PQ 

COMMON/-.ll/CR,LP.cQ 
COMMON/ 311/ I MX ,t MX 

iO  FORMAT  ("j",":_RROF  ( N R F ) « MAXIMUM  ITERATION  LIMIT  E XC£  EO£  O’V/ ) 

X M = X 

JO  i J i II=i,IMX.l 
F = F F ( X M ) 

FD  = FOF ( X M ) 

XP=XM/2  IIFfFJ.M.ul  XD  = XM-F/FO 
0 1 F = A B S ( X c - X M ) 

ERK=JIF  ( XM.NE . 0 > ERF=GIF/xM 

IF  URR.LT.EMX)  GO  TO  IOC 
X M = X P 

1 Cl  CONTINUE 

HRIU  (L  = , 11  > 
lOu  CONTINUE 
NPF=<P 

RETURN 

END 

SUBROUTINE  SRM < KN , N, RTS , M AX  I T , E PI , E »2 , FN , FNR£ AL) 

; MULLER  METHOD 

; He  Na ICI , "e  LEMc  NTS  OF  NU MF  R I C A L A N A L Y S I S" , PP.  US-2 

: MULLER. "A  M.THOO  CF  SOLVING  ALGEBRAIC  EQUATIONS'* 

: MATH  TAtlLeS  AND  AIDS  TO  COMPUTATION  iO  (1956),  °P . 208-215 


COMPLEX  P T S ( 1 ) 

logical  fnreal 

COMPLtX  XC,X1,X2,FXC,FX1,FX2,RT,FRT,FPPT,0EN 
COMPLEX  A2,  U2,  Cl,  02, E2.F2.H2, 02, Q2l 
E°51=lP1 
£PS2=_P2 

IF(iPsl.Lt.U.u)  cPSl=i.CE-i2 
IF<£PS2.Li.o.C)  E°S2=l.CE-20 
<NL=<N 
NL=KN*N 
10  CONTINUl 

IF(KNL.LT.Nu)  GO  TO  2C 
RETURN 

20  <NL=KNL+1 
< 0 U N T = u 

iO  X L =R Tb ( < N l • - u . 5 
ASSIGN  4 J TO  NN 
RT  = X0 
GO  TO  SO 
40  FX0=FPRT 

X1=RTS(<NL) *0 .5 
ASSIGN  50  TO  NN 

RT  = X1  A -48 

GO  TO  SO 
50  FXi=FPRT 


o o o o o 


r 


X2  = ’?TS  (<NL) 

ASSIGN  feu  TJ  NN 
nU=xc 
GO  TJ  10 
oO  FX2=FPRT 
H2=X2-X1 

32= (X2-X1) /(Xi-XC  ) 

65  CONTINUE 
Q21=a£+i. C 

A2=a2*(FX2-a2:*FXl*Q2*FX0) 

32=  < JZl+Ut  >*fx  2 -Q2l*021  *Fx i ♦ Q2*02*FX 0 

C2=02l*Fx2 

32=32*32~l.U*A2*C2 

IF  ( < FNkc.  A L ) . ANO.  (<?£AL  <021  .LT.O  . 0 > ) 02  = 0. C 
02  = C S IP  T ( 02 ) 

^2=32+02 
F2  = 32 -32 

IF<CABS(F2).GT.CAPS<£2) ) E2=Ft 
IF(CA3S(Ll)  .GT.u.O)  GO  TO  7 Cl 
£ 2 = 1 • u 
70  CONTINUE 

Q2=-2.0*C2/i2 
FX  0 = FX I 
F X i = FX  2 
H2=H2* 12 
?T=^T+H2 
ASSIGN  75  TJ  NN 
GO  TO  1C 
75  CONTI  N'Jt 
FX2=FP*T 

IF  ( CA3S<Hu  ) . LT  ,t.P3i*  C A3S  (FT)  > GO  TO  i„C 

IF < < CA3S< F» T) .LT .F0S2I . ANO. <CA3S< FPFT) ,cT. £=>32) ) GO  TO  ICO 
IF < KJUNT . GT. MA X1T ) GO  TO  ltiC 
IFICABS<FX2> .LT.lt. u*CABS(FXl>)  GO  TO  ou 
H 2 = 0 . ?*“2 
92=J .5*32 
Tr  = T-?-H2 
10  CONTIN'.Jl 

<ount=<oijnt*i 

C A c L F N ( R T , F •?  T ) 

FPF,T  = F’?T 

IF(<Nu.lT .21  GO  TO  OX 
00  30  I = 2 • X N L 
JtN  = -?T— ?T5<I-1) 

IF<CA3S<3tN) .LT.FFS2)  GO  TO  02 
90  FPP.TsFPPT/OEN 

31  CONT INU_ 

GO  TO  NN»(<*U»5ii*6"*75> 
loO  CJNTINUt. 

xTS(<NL) = -<t 
GO  TO  . 0 

32  XT  = 3T*i).Cul 
GO  TO  10 
END 


SJJ  X03TI  Nt  V JM  c ( CFO.  AFC,  , 3F,  l.F  , HFM.HT  J) 


OlSSuw  functions 

NEU3AN3  FUNCTIONS 
3SN<tt  FUNCTIONS 


INTEGf?  0-0. OA 
30D 

C3M»_..<  Cl 


A -49 


I 

J 


oooo  o o oooo  ooo 


V 


c 


c 

c 

c 


COMPLEX  A G 
G 0 M P l.  L X jF»'IF,HFI*f‘-iF  = 

DATA  POtt/x.u  *2*./ 
DATA  DPL /L.oJ./ 

DATA  c.o/ 

0A=IA&3  (D  O) 

100=CA JS (AkG) 


DF  = »j  5Ip(OA.-5.C)  PC=F* 

NF=-PJLi 

HFrt=3F-CI*N- 

H F 3= 1F«-CI*NF 


c 


c 

c 


c 


I r ( IOJ.^T-0  » 

IF  I'ICQ.lT  • "3  > 

IF  ( ( 'IUD.Gl  ..Ptl 

IF 


A A3  . ( fOC  . L-  . I FUI  I 

r- : o . r,  t . b 3 j » 


p_  rjPN 

C«  LL  SRAS(0A,ARG,dF,NF,HF*1,FIFP) 
CAUw  5 RRS  ( J A * A P G t E F y NF  f H F M « HF  P ) 
^AAw  3 P A L < 0 A • A p G « FF  * N F ,F<rM,HFP) 


IF  < JRH.  AT  . u > A.H'PN 

3 F = ( ( -1 ) * *u A I * 3F 
NF=  ( (-*  ) **l)A  ) *NF 

HFM  = JF'-CI*NF 
MFP=3F+CI*NF 


■?l  r jF;-j 
END 


£**********♦****♦*********«*#****#♦♦**♦*****#***♦******¥****•**♦**•**♦•* 

SUJPjjriNc  SFUS  <0*0, AFG.BF.NF. HFM.HFo* 


ASYMPTOTIC  EXPANSION  (SFAlw  ARGUMENTS) 


INTlj- < JPO 

0 0 M p L I X Cl 

COMPutK  ARG.AOii 

CC M °ll  < BF,  4F,HF^,HFP 

DATA  Pl/3.  Ui5  DtfcJ/ 

CI=CM3lx(C.J«I.G) 


A D£ - A*  G/I 

FP=FF ( OP  0 ) 

FM  = FF(jp  D-il 

3F  = + (AOE**  ( + JitlD)  ) /FO 

MF  = -(AOI**(-ORD))*(FM/ai)  «IF(JfiD.FD.J|  NF  = ( 2 /PI  ) • CL  OG  ( A RC,  I 
HFM=DF-CI*NF 

HFP=8FfCI*NF 


PETU9N 

CNO 

**********»*********¥¥**********•***•« 

SUa^OUTlNl  J->To(OFfl,lFG,0F,NF,HFM(Mf3| 

ASCENDING  SEPI-JS 

INT c.Gi  P JP0.0Mi.0  = i A*5 

COIPR:  < Cl 


A 


o o o 


c 

c 

c 

c 

c 

c 


c 


c 


c 

c 


c 

c 

c 

c 

c 


COUPLE < UPC ,A02 
COMPLt-X  f M,TP 
COMPLEX  SP,SN,3Y 
COMPLEX  3SF.NoF.Y3F 
COMDLtX  yf.nfm.xfp 
COMPLEX  8F.NF.MFM.HF0 

COMMON  /Oil /I  MX  ,t  HX 

OATA  PI/3.i*.i382b3/ 


CI-CM®LX ( C . j . b ) 

0M1 =ORO-l 
OPl=URO*i 

A02  = ARG/2  JTM  = A02**  ( -ORU  I T P = A 02  * * ( +0R J ) 

FO  =FF  (0*01  3Gl  =GF(  1) 

F0M1-FF ( 0 M 1 ) JG0P1=GF  (0P1 ) 

G0=G1FG0P1 

SP=9SF ( OPP, ARG) 

SN=NSF(0-<0,ARG,GO» 

S Y = Y SF (ORD.ARG) 

dF= ( TP/FO) »53 
YF=(2/PI)»CLOG<an2)*riF 
NFM=-(TM/PI)*FJM1*SY 
NFP  = - < TP/PI ) * (1/FCI *SN 

nf=yf»nfm+nfp 

HFM=BF-CI*NF 

MFP=BF+CI*NF 


rf  turn 

tNl) 


subroutine  sral(gro,acg,pf.nf,mfm,hfpi 
AsrMPrjrxc  expansion  < large  arguments) 


INTEGtR 

COMPLEX 

COMPLEX 

COMPLcX 

complex 

COMPLEX 
COM3LF  X 
complex 
C 0 M p i.  c X 
CONPLl X 
C.0MPl_  X 


d-0 
A-G 
Cl 
cc 
x .c 


C X . j X 
A MX , APX 
c M x , f < 
p ,o 

3F,  NF.HFM.HFP 

jcexp 


J A T A »I/3  .mi-3  1263/ 


CI=CMPLX (C.Otl.u) 

X = 6RC- <2*0RJ  ♦..)*(  = I/A  > TCx  = Cc03  f J1  * *_><=. o 
CC  = 2/(°I,,6R&>  iC  = CSOP  T f CC  I 

AMX  ~“CI*X  t_ MX=DC1 XP (AMY ) 

A PX  - *C  I * X i_P>  = CCZx°(APX  ) A*51 

C A L e SRFOCORO.AAC-.^  • " ) 


I N ( X ) 


JF  = C* ( = *CX-  I*Sx ) 

NF  = C*  < P*SX* I'LXI 
MFM=C* ( P-CI*1) *eMX 
HFP=C*  ( o«-Cl*Q»  *cp* 

RETURN 

END 

subroutine  srpq fo=o, a=g,p.q) 

HA  N<Ec  ASYMPTOTIC  EXPANSIONS 

INTEGER  CR,L=,pQ 
INTEGeR  0=0, UM 
I NT  E Ge  R sgn 
COMPLEX  A=G,AT3 
COMPLEX  T,s 
COMPLEX  P,Q 

COMMON/Aj.C/CR»LP»pQ 

COMMON/3i,./IMX,EMX 


>*««+«***  + *» 


iO  FORMAT  ("C'VeRROP  (SPPO 

° = 1 k3  = . 

Q=C  i£Q=i 

T = 1 

OM  = -,*ORO*GRO 
AT8=S* APG 


MAXIMUM  XT_-’ATION  LIMIT  E XCEE  OE  O’V  / » 


10  2 

1 C 1 
I Oil 


DO  101  IX=1,1MX,2  IIXPl=Ix»l 
DO  102  IY=IX,IXP1.1 
SGN=-(-il**IY 
IO=2*IY-i  $IS=IC*IO 
S=SGN* ( (OM-IS)  / (IY*AT#»)  ) 

T = T * S 

IF(SGN.EQ.-i)  P=P*T 
IF(SGN. EO.+l)  1=0+T 
CONTINUE 

I F ( CA  8S ( ° > . NE • 0 J e°=CAPS(T/P) 

IF(CAB3(a»  .NE.ii)  EQ  = CA0S(T/Q) 

FRR  = A MAXx ( t = » E Q • «IF(fRR.LT.EMX»  GO  TO  1 JO 

CONTINUE 

MRIT£<lp.1u) 

CONTINUE 


RETURN 

END 


C 

c 

c 

c 

c 


SUBROUTINE  MBNHF ( CRD*  APG* MBF, MNF, MHFM.MHFPI 

MOD  I F IEO  BESScL  FUNCTIONS 
MOOIFIEO  NEUMANN  FUNCTIONS 
MOD  I F IE  0 HANKeL  FUNCTIONS 

INTEGER  0=0 
COMPLEX  Cl 
COMPLEX  APG » C A 
COMPLEX  BF.NF.MFM.HFo 

COMPLEX  MBF.MmF, mhfm,mwff  a_52 


DATA  PI/s.l^lS-JEfc  :/ 


oooo  o oooo  o ooo 


c 


Cl =CMPlx  < o . o 


c 

c 

c 


CA  = CI*A\1  ICALL  BNHF (0RG,CA,8F,NF,HFM,HFP) 


Mi1F  = (Cl**  .’-OR J) ) *3F 
MNF=NF 

MHFM=(PI/2>*<CI**(0R0+1))*HFM 
M HF=>=  (3I/E)*(CI*M0RD*1>)*HFP 


RETURN 

END 


COM  °LE  X FUNCTION  3 3F ( OR 0 , A RG ) 


UcSScL  ->E^i£0 


I NT  i.Gc.  R CR,lP,RG 
INTEGER  3RD 
COMPLEX  A-G,ADE 
COM^l.X  TC»TR,r8»S3 


COMMON/AiL/CRf LP.c0 

C0MM0N/311/IMX,£MX 

i.0  F OR  MAT  (‘‘u',f"EFRO"  (HSF)  S ^ A X 1 MU  M ITERATION  LIMIT  EXCELOPD"//) 

ADt  =ARG/,_  *TC=AC1**E 

T B = l I38  = i IEB  = i 

30  « o i IX-iiiMXii 
T R = - TC  / (IX* ( IX  + ORC) ) 

TP=T8*  TR 

GB=S J*T  x 

IF(  CADS  (i°)  . Mi.  . 0 ) :.P  = GAPS  (TD/33)  J I F ( a 9 . L T . E MX  ) GO  TO  10o 

101  CONTIN'Jr. 

WRITE(L^»>.u) 

10  J CUNT  IN’Ja 
8SF=3n 


RETJRN 

END 

CUMPcc-X  FUNCTION  NSF  (ORD,  AkG,GO> 


NEUMA  4 iE  - I-G 


C 

c 

C 


INT EGlR  0”  tL 
INTEGER  ORO 
C0Mdll X AFG.AOB 
COMPLfcX  T i , T 3 
COMPLEX  TC,T R,TX , TN,3N 

C0MM0N/A.,.u/CR,Lo»ca 
COMMON /Dll /I Mx  , CMX 

10  FORMAT ("0 ", "ERROR  (NSF)  S MAXIMUM  ITERATION  LIMIT  EXCEtOEC  "//) 


ACE  = ARG/E  $rC  = ACE**C 

TS=GJ  SSN=TS  IEN=: 

TP=  1 
TN=T  J 


C 


A-53 


oooo  o ooooo  o oooo  o ooo 


I 


1 


c 

JO  • u x I A - i . J.  '-I A « I 

TX  = *F;_OAT ( . ) /IX ♦FLOAT (1) / f IX*0R0) 
rs  = r.>«-T* 

tr=-tc/ < i/*  ( ix  ♦os  cn 

IP-  TP*TR 
TN=TS*r- 
0N=oN+fN 

IFCCAFSISU  .Ur  .C)  E N - C A H S ( T S / 3 N ) ?I  F ( r.N . L T . E MX  I GO  TO  100 
lGi  CONTiNlI- 

i 0 0 CO  NT  IN')' 

NSF=3N 

C 

RETURN 

JNu 

COM^l^X  FUNCTION  YSF(C~C»ARG) 

NEUMANN  3_kI_3 

INTEGER  JrOrJM* 

CO  M 3i_  _ X A f G • A U E 
COMPl..  < TC.Ta.TX.SY 

YSF  = j $IF(ORO.£<J.G>  RETURN 
YoF  = i J IF ( OR  0 . r Q . 1 ) RtU'PN 

JM1 = OR J-l 

AOE=A rj/I  ITC=AQt*«I 

TV-j.  »5Y  = * 

00  1 3c  IX  = 1, 3Mi,i 
TR  = + Tc/ < IX* ( JR3-I X) ) 

TY=TY*T? 

SY=SY+TY 
103  CO  NT  I N JE 
ysf=oy 

RETURN 

END 

****»***»***•*»♦*¥♦»»*»*«»»*♦*•»*»***»»****»*»»*♦**»*»»*•»»»*****»»»»»» 

REAL  FUNCTION  ff  ( I rc » 

FACT0PI A^ 

FF=1  IIF(IR.LT.i)  RETURN 
F = i 

oo  m ix  = ifi"*fi 

F=F*IX 

111  CONTI N JE 
FF  = F 

RE  T JRN 
END 


REAL  FUNCTION  3F(IRI 

GANMA  a-5A 

JAM  .E77ElG'56‘.cLl537/ 


I 


oroo  oo  o oo  o oooo  o o ooo 


C 

GF=-GAMMA  fIP(IR.£3.i>  Pf  TURN 
C 

SUM=u 
IR  Mi.  - 1 R-i 

oo  u(j  xx  = i,  irmi»i 

TERM=FLOAT (1 ) /IX 
SUM=SUM+T;RM 
100  CONTI NU£ 

GP= -GAMMA ♦ SUM 


C 


RETURN 

END 


C* 


♦**¥*»*¥*****♦♦***, 


'++*++*+ * 


REAL  FUNCTION  CAN3(C> 


COrtPCcX  ANGLE 


R l A L R • I 

COMPLlX  C 


CANl=L 

R=  RlAC ( C ) SlMfi.EQ.GJ  RETURN 
I=AIMAG(C) 

CANG-AT  ANE (I *R» 


RETURN 

END 


»»**»*»**», 


COMPLEX  FUNCTION  DC£XP(C» 


DOU  Ul - =RrCIS10N  CO  M°LF  X EXPONENTIAL 


oouaut  precision  dr.di 

D0U3lE  PR- Cl j I ON  0E,0C.DS 

COMPllX  CC 
COMPllX  C 
COMPLEX  EXP 

CC  = CRPlX  J i ,ii 0 » 

CR-  RlAL(C)  iO^=nnLLICP) 

£ I=AIMAG<C)  S0I  = 0(JLf  (CI> 


Jf  = JExPfUO  i Gl  = SNGL ( DE I 
JC=OCOi(JI)  TSC=SNGL ( OO 
J5=OiIN<01)  Jc3=SFGL (OS) 

EXf>  = 3'-*  f3C+0C*SS)  FOCt  XP  = r XP 

RETURN 

ENO 


SU3R JUT  Inl  FTJRiniJ.r)ATA,Px.PYt 


FAG 

TRAN 

0~UM  R-Au 

INT 

EGc  R 

p>  »py 

I N T 

EG-R 

C-  ,Lj,30 

INT 

EGlR 

OU 

INT 

lGlR 

OMN 

INT 

lp^R 

- 

INT 

EG-  p 

no,  :of 

A-55 


oooooo  o ooo  o o o o o oooooo 


I N 7 _ 'j  t ? **  H N J • W A I T 

coh-»lc.x  oata 

c 

JIM^NSION  JA  T*»  ( aX  « aY  > 

c 

c 

jata 

jara  -o  / 1/ 

)ara  aut/ls/ 

c 

*0  r CRMAT  CV,  *R_  A3  T k A N3FE  F COMPLETE  / CKUM  : , 1 1 1 ) 
„1  F ORMAT  ( " j",  ,*->L  ao  r-ANSFff-  _F?On  / UNIT  »'*,riC/ 


♦ ••  / E.MD  i'Miu/ 

♦ / EOF  ,I1C) 

c 


t_Nr=0MN  T~CF=OMN 

C a c L NTR4  NO'JtRHN  C» 

CJiL  ntran ( )j  ,-?c  .eng,  data, eof> 
can.  ntranzdiwhait) 

rlaj<ou>  uflTa 

IF(-_OF._  J.E'JJI  Hi-  ITc  Ue.lt  ) CU 

IF  ( E OF . NE  . JO  I WRITE  ( L c » i 1 ) DU*END,EOF 

■it  TJR.N 
END 

**»*»*»****»«»»♦*»»»»**»*«*»**»***♦»»#»**»***»**»**♦*»*»**»*»»*»»»»***» 

sugROOTiuf  ftow<ou,data,°x,°y> 

FA3TRAN  OkUiH  WRIU 
INTEGER 

ZNUGiR  CF,LP,o0 
I N T c.  Gc.  R 01 
I NT  GE  R OMN 
INTtG-R  HFT 
INTEGER  -.NO.tOF 
INTEGER  -WNO.WAIT 

com jLt.x  oata 

dimension  oara  <fx ,rvi 

COMMON/aiL/CR,L°,crl 

OATA  RWNO/1U/ 

DATA  WRT  / 1/ 

DATA  mMT/LcL/ 

1C  FORMAT  ( "0  " * **HkI  T f TRANSFER  COMPLETE  / 0~UM  f‘*,riGI 


ii  FORMA  f ( "3 ",  “Hf- 1 Tt  TRANSFER  ERROR  / UNIT  I'MlC/ 

♦ “ * V / END  l” .110/ 

♦ '*  *V  / EOF  : " . 1 1 G » 


OMN=PX*PY  SENJ=CMN  tE  of  = DMN 


c 


CALL  NTRAN(OU.RWNC) 

CALL  NTRANI JJ.WRT.FNO, OATA, EOF) 
CAL»  NTRAMOU, HA  I T ) 

WRITEOU)  JATA 


A-56 


1 


**  y 


o o o 


UL-i.0 


APPENDIX  B 


0.0 


icj  : 
11.0 
l.G 

♦ 3 

4 


20.0 


Sample  Input  Data  Listing 


+ 1 


0. 11 J£-b  J 


NA 

PULSE  (At 
TIME  Co) 

i 

0.0 

FIGURE  2 
NA 

PULSc.  < A » 
ESC OUENCY 


l 

2l.Lw+'u3 
LIGHT  MING  P'JlSE 


( RAD) 


(FREQUENCY  DOMAIN)  VS  T R CQU-l  NC  Y 


- 1 • 42  ♦ • 3t> 

FIGURc.  3A  ! OUT  _ 
FREQUENCY 
EXCITATION  <V/M) 
POSITION  Ml 


R ^AGITATION  (FREQUENCY  DOMAIN)  VS  POSITION 

B-l 


IMX.EMX 
ML.RL, SL 
MIN, MAX 
NO. MO, MX 
NQ.NI 
MN.NM 


1 

i.  1 

1 1 

1 

M»S 

0 

0 0 

0 0 

0 

RPS 

0 

1 1 

1 1 

1 

WPS 

3 

MR 

1 

MP 

luu.bE-Co 

1. CE+06 

TM,  WM 

51 

51 

NK,  NC 

10  Ou  • (j  1 C 0 J • l l j 0 0 • 0 

CG 

1 . 0 z ♦ 0 j j . 

5 0.5 

a c 

PS 

20.0E+C30. 

5 u . 5 

G 1 

RX 

1 . 0 E ♦ u 30 « 

5 0.5 

o c 

RS 

0.0 

1.0E-J6  103.0E-0F 

l.CE-06  i.uE“06 

PS 

0.0 

l.lL-Jb  HC.Oi-Cfe 

1 • o £- 0 6 1 • b E-0  6 

RX 

c.o 

1 . 0E- j 6 iOC.CE-0  6 

40.CE-03  12Q.0E-23 

RS 

-2.J00E+C  3 + c 

.U00Et-02  41 

EI.EF.NE 

-0 . 0 OOt  *0 3+i 

• 0 CCt  + O 3 4l 

XI , XF, NX 

-0 . 0 0 U E ♦ u i*  ± 

. uO  vE»  0 3 3 

YI.YF.NY 

l.OE+tJ 

♦■1 . C t-  + J6  61 

WI.WF.NW 

0 • 0£-b  3 

l.CE-03  61 

TI .TF.NT 

-1 

-1  -1 

♦ 1 -1 

MEXYWT 

1 

•*1  U 

IEIFN 

1 

4 1 1C 

IXIFN 

1 

3 1 

IYIFN 

i 

bx  21 

IWIFN 

1 

o 1 1C 

ITIFN 

RW 


2.95 3£- 

C 3 3.953E-03 

RSM.RSP 

4.9532- 

03  5.953E-03 

RAM, RAP 

6 • 9532- 

03 

RC 

1.0 

OO 

0.0 

SO 

1.0 

5.0 

1.5E-G3 

RUE  OE 

i.  a 

2.26 

l.CE-03 

RUEOI 

1.0 

1.0 

5.72+P7 

RUEOA 

1.0 

2.26 

l.GE-03 

RIJE02 

1.0 

1.  J 

5.72407 

RUEOS 

1.0 

2.26 

1. 0E-C  3 

RUEOI 

1.0 

1.0 

5.72*07 

RUEOW 

50.u 

J . 0 

ZL1 

50.o 

b . u 

ZL2 

1 

MIO 

1 

ML  C 
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1 


2 
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FREQUENCY  (RAO) 

1 4 0 0 0 
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CURRENT  (A) 

POSITION  (Ml 

0 0 112 
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CURRENT  (A* 
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FIGURE  5 S TRANSFER  RATIO  (FEE  QUENC Y OOMAIN)  VS  FREQUENCY 
NA 

TRANSFER  RATIO 
FREQUENCY  (RAO) 


1 l 1 1 2 

-0.0290  *4.0358  -1 

FIGURt  6 A : INN-.R  EXCITATION  (FREQUENCY  OOMAIN)  VS  POSITION 

FREQUENCY 

EXCITATION  (V/M) 

POSITION  (M) 


j 


2 


-l.15t.-Qo  *1.22E-0i  -1 

FIGURc  6e  « INNER  _XCITATION  (FREQUENCY  OOMAIN)  VS  FREQUENCY 

POSITION 

EXCITATION  (V/M) 

FREQUtNCY  (RAO) 

1 o 0 0 0 

-1.15£-0o  *1.22E-0o  -1 

FIGURE  7 A « INNER  CURRENT  (FREQUENCY  DOMAIN)  VS  POSITION 
FREQUENCY 
CURRENT  (A) 

POSITION  ( M ( 


4 0 C 0 C 

-Q.Q0189  ♦,4,.Co?22  -1 

FIGURE  70  » INNER  CURRENT  (FREQUENCY  OOMAIN)  VS  FREQUtNCY 

POSITION 

CURRENT  (A) 

FREQUtNCY  (-.AO) 


2 


— 0.00189  *4.40522 

FIGURE  7 C I INNER  VOLTAGE 
FREQUENCY 
VOLTAGE  (V) 

POSITION  (M) 

4 4 

-Q.C05oC  *0.4^53 
FIGURt  70  ! INNER  VOLTAGt 
POSITION 
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FREQUtNCY  (*A9) 

1 4 

-0.J05o0  *4.Clo3 

FIGURE  8 A I INNcR  CURRENT 

TIME 


-1 
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0 4 0 

-1 
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-1 

(TIMt  OOMAIN)  VS  POSITION 


B-2 


MTPGW  E 


MTPGX  0 


MTPGW  0 


MTPGW  T 


MTPGX  E 


MTPGW  E 


MTPGX  I 


MTPGW  I 


MTPGX  I 


MTPGW  I 
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j o 

-o.u  **2. 5 
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POSITION 
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- J . u ♦! 2 • i 

FIGURE  dC  * INNER  VOLTAGE 

TIME 

VOLTAGE  (V) 

POSITION  M> 

J O 

-U.333  *39.3 

FIGURE  *0  * INNER  VOLTAGE 

POSITION 

VOLTaGE  CV) 

T I Me.  (3) 


C 0 U 

-1 

(TINS;  DOMAIN)  VS  TIME 


i 1 z 

-1 

(TIMt  DOMAIN)  VS  POSITION 


f U 0 

-1 
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MTPGT  T 


mtpgx  t 


MTPGT  T 
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C - 1 


Table  1 


List  of  Input  Cards  and  Formats 


Card 


Label 


format 


1 

1 MX , JMX , EMX 

21 10, CIO. 0 

2 

ML,RL,SL 

1 10, 2G 10.0 

3 

MIN,KAX,MP 

2C10.0,  1 10 

4 

ND.MD.MXX 

3 1 10 

3 

NQ , N 1 

2110 

6 

IFT.JFT 

2 110 

7 

NM.MN 

2 1 1C 

8 

MLS,MC0,MTT,MC 1 ,MTL,MFT 

6110 

9 

RLS  , RCO, RTT , RC 1 , RTL , RFT 

6110 

10 

WLS,WC0,WTT,WCI  ,WTL,WFT 

6H0 

1 1 

HR, MG 

2110 

12 

NR,  IR 

2110 

13 

TPK 

G10.0 

lit 

t'FK 

G10.0 

1 8 

HQ 

i'r 

16 

TH,WH 

2G10.0 

17 

NK.NC 

2110 

10 

PGPS.PGPX ,PGRS 

3G1C.0 

1? 

PSPS , DCPS , DKPS ,MPS ,NPS 

3G10.0.2I  10 

20 

PSRX , DCRX , DKRX ,MRX , NRX 

3G10.0.21 10 

21 

PSRS,DCRS,DKRS,HRS,NRS 

3G10.0.2I  10 

22 

TSPS,TRPS,TDPS,TPFS,TCPS 

5G10.0 

23 

TSRX,TRRX,TDRX,TPPX,TCRX 

5G10.0 

2 it 

TSRS ,TRRS ,TDRS ,TPRS ,TCRS 

5G10.0 

28 

SI  ,SF,NS 

2G10.0, 1 10 

26 

XI  ,XF,NX 

2G10.0, 1 10 

27 

Y 1 , YF , NY 

2G10.0, 110 

28 

Wl ,WF ,NW 

2G10.0, 1 10 

23 

T 1 ,TF,NT 

2G10.0, 1 10 

30 

MS, MX ,MX ,MW ,MT 

5110 

31 

IEI  , 1 EF.NEI 

3110 

32 

1 X 1 , 1 XF , MX  1 

3110 

33 

IYI , IYF.NYI 

3110 

34 

IWI , IWF.NWI 

3110 

35 

ITI , ITF , NT  1 

3110 

36 

RW 

G 1 Oo  0 

37 

RSM.RSP 

2G10..0 

30 

RAM, RAP 

2G10.0 

39 

RS 

G10.0 

i<0 

20 

G10.0 

it  1 

Y0 

G10.C 

42 

URE,ERF,0F 

3G10.0 

43 

UR  1 , ER 1 ,01 

3C10.0 

44 

URA , ERA ,0A 

3G10.0 

45 

UR2 , ER? ,02 

3G10.0 

46 

UR$,ERS,0S 

3G10.0 

47 

UR  1 ,ERI ,01 

3G10.0 

40 

URW,ERW,0W 

SG10.0 

49 

ZLM 

2G10.0 

50 

ZLP 

2G10.0 

C - l 


C-, 


UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

x 

UJ 

X 

X 

X 

X 

X 

UJ 

UJ 

UJ 

UJ 

UJ 

X 

X 

X 

X 

X 

— 

X 

— 

— 

— 

— 

— 

X 

X 

X 

X 







. 

or 

z 

— 

or 

or 

or 

or 

or 

Z 

z 

z 

z 

Z 

— 

— 

— 

— 

CL 

or 

or 

or 

or 

z 

Z 

z 

h- 

or 

or 

h- 

h- 

v- 

h- 

h- 

CL 

or 

or 

or 

X 

CL 

CL 

01 

CL 

or 

K 

h- 

*- 

CL 

or 

or 

CO 

ID 

>— 

CO 

CO 

co 

CO 

CO 

O 

ZD 

ID 

ID 

ID 

H* 

h- 

h- 

h~ 

h~ 

CO 

CO 

CO 

10 

CO 

Z> 

ZD 

ID 

ZD 

UJ 

H- 

CO 

UJ 

UJ 

UJ 

UJ 

UJ 

1 — 

1 — 

l — 

h- 

to 

CO 

CO 

CO 

CO 

UJ 

UJ 

UJ 

UJ 

UJ 

»— 

*— 

h- 

CL 

UJ 

UJ 

or 

or 

CL 

or 

or 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

or 

or 

or 

or 

or 

J UJ 

CL 

or 

or 

CL 

a_ 

CL 

Q- 

CL 

or 

or 

or 

CL 

or 

or 

or 

CL 

or 

or 

Q- 

a 

Cl 

o_ 

Q_ 

CL 

or 

CL 

or 

lD  L- 
£ O 


a:  

UJ  Q_ 


o o 
o o 
a:  or 


g -i 

g => 
3: 


D ID 

c.  a. 


— o o 

_i 

a.  a:  a; 

— UJ  UJ 
h-  CO  CO 

-J  XX 
ID  Z3  ZD 
X z z 


1/3  CO 
CO  CO 


ID  ID 
O O 


ID  ID 
cl  a. 


ai  CO 

< z 
uj  O 
a-  cj 

UJ  >- 
CO  < 
—I  CD 
ID  UJ 
CL  Q 


DOW  QJ  QJ 

cn  cn  cn  id  cri 

<D  03  03  • — • — — <d  , , , _ 

rccSSSS*J'^Jr0f0',,ratl, 

CU£_a»UJ(U(ycica)<Da)<DCJ 

u u J_  — — u u ^ u u 


H 

H* 

I- 

h- 

h- 

z 

Z 

z 

Z 

Z 

h- 

< 

< 

< 

<t 

Z 

z 

( — 

CL 

K 

f- 

CL 

H 

1 — 

or 

< 

< 

CO 

>- 

UJ 

X 

CO 

CO 

>- 

LU 

X 

co 

CO 

>- 

UJ 

> 

UJ 

H- 

>- 

UJ 

1— 

z 

< 

CO 

< 

z 

z 

<c 

QJ 

< 

z 

z 

< 

00 

or 

UJ 

<c 

CO 

to 

or 

UJ 

CO 

to 

0 

CD 

X 

UJ 

0 

CD 

CD 

X 

UJ 

0 

0 

CD 

X 

CO 

<D 

1 

■r** 

< 

CO 

-J 

z 

CD 

UJ 

ID 

CL 

<D 

LD 

UJ 

Z) 

CL 

CD 

0 

UJ 

ZD 

h- 

•— 

UJ 

ID 

b 

H- 

— 

UJ 

3 

0 

0 

Z 

0 

z 

O 

Z 

CO 

CL 

O 

CL 

CD 

CO 

or 

O 

CL 

CD 

>- 

UJ 
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